Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 167,291 G→A 5.2% intergenic (+27/‑193) mrcB → / → fhuA fused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,2910GA 5.2% 274.2 / 2.8 95intergenic (+27/‑193)mrcB/fhuAfused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (52/38);  new base A (3/2);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATC  >  NC_000913/167145‑167439
                                                                                                                                                  |                                                                                                                                                    
aaTCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAg                                                                                                                                                    >  1:234431/1‑149 (MQ=255)
 aTCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGt                                                                                                                                                   >  1:522607/1‑149 (MQ=255)
          aTCAGTCTTCTCCGCCGCAGCAACAGCCGCAACAGCACCCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGATTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATcag                                                                                                                                          >  1:1182836/1‑147 (MQ=255)
          aTCAGGCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCAGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCg                                                                                                                                          <  2:519411/149‑1 (MQ=255)
             tgTCTTCTCAGCCGCAGCAACAGCCTCAACAGCAACCTGCTCATCAAGAGCAGAAAGACAGAGACGTTGTAGCAGGTTTGATCAAGGATATGTTTGTTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                       <  1:575709/148‑1 (MQ=255)
             aGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCGCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTa                                                                                                                                       <  2:37810/149‑1 (MQ=255)
               tcttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaa                                                                                                                                     >  1:1100408/1‑149 (MQ=255)
                cttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaaa                                                                                                                                    <  1:56146/149‑1 (MQ=255)
                 ttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAAt                                                                                                                                   <  2:1020210/149‑1 (MQ=255)
                   ctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCCGCCAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGAGATGTTGGGTAGGAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCTTCCGGTAATATaa                                                                                                                                 >  2:246125/1‑149 (MQ=255)
                   ctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATaa                                                                                                                                 >  1:932431/1‑149 (MQ=255)
                    tcAGCCGCAGCAACAGCCGCCACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGGGTAGCCGGTTGGATCATGGATATGTTTGGTAGTAATTAACATCGAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                >  1:664125/1‑149 (MQ=255)
                    tcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                >  2:1139176/1‑149 (MQ=255)
                      aGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACAt                                                                                                                              <  2:1139048/149‑1 (MQ=255)
                       cccGCAGCAACAGCCGCAACAGCAACCTGCTCATCAAGAGCAGAAAGACAGCGACGGTGGAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATc                                                                                                                             <  2:799755/148‑1 (MQ=255)
                        ccGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATcc                                                                                                                            >  1:303917/1‑149 (MQ=255)
                         cGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATccc                                                                                                                           >  2:427924/1‑149 (MQ=255)
                          gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCa                                                                                                                          >  2:388254/1‑149 (MQ=255)
                              caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                      >  1:1122801/1‑149 (MQ=255)
                              caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                      >  2:345953/1‑149 (MQ=255)
                              caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                      <  2:289256/149‑1 (MQ=255)
                                aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                    <  2:207314/149‑1 (MQ=255)
                                aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGACTAGCCGGTAGGATCAAGGATATGTTTGGTAGAAATTAACATCTAAGCGTGAAATACCGGAAGGCGAGTTGCCAACCGGAAAAATAACATCCCAACTAAg                                                                                                                    >  1:209562/1‑149 (MQ=255)
                                    ccGCAACAGCAAGCTGCTCAGCAAGATCAGAAAGACTGCGACGGTGTAGCCGGTTCGATCAAGGATTTGTTTGGTAGTTATTATCATCTAAGCGTGAATTACGGTATGGCGAGTTGCCATCCGGTAAAATATCATCCCATCTAAGatat                                                                                                                <  2:702490/149‑1 (MQ=255)
                                    ccGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGatat                                                                                                                >  1:8247/1‑149 (MQ=255)
                                     cGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGTTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATAtt                                                                                                               >  1:972805/1‑149 (MQ=255)
                                      gCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTa                                                                                                              >  1:1106539/1‑149 (MQ=255)
                                      gCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTa                                                                                                              >  1:327116/1‑149 (MQ=255)
                                           aGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGACGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACAct                                                                                                         >  1:1021659/1‑149 (MQ=255)
                                                ccTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttctt                                                                                                    >  2:291176/1‑149 (MQ=255)
                                                ccTGCGCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttctt                                                                                                    <  1:610379/149‑1 (MQ=255)
                                                  tGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCtttt                                                                                                  <  2:554550/149‑1 (MQ=255)
                                                   tcTCAGCAAGAGCAGAAAGACATCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGAAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTc                                                                                                 <  1:1036379/148‑1 (MQ=255)
                                                         cAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTgg                                                                                           >  1:552380/1‑149 (MQ=255)
                                                           agagCAGAACGACAGCGACGGTGTCGCCGGGTGGATCAAGGTTATGTGTGGTTGTATTTAACATCTAAGCGTGAAATACCGGATGGCGAGGTGCCATCCGGTAAAATAACATCCCAGCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                         <  1:696255/149‑1 (MQ=255)
                                                            gagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttg                                                                                        <  2:399064/149‑1 (MQ=255)
                                                              gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCTGgttgtt                                                                                      >  2:646927/1‑149 (MQ=255)
                                                              gCAGAAAGACAGAGACGGTTTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCTGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                      <  2:891729/149‑1 (MQ=255)
                                                              ccAGAAAGACAGCGACGGGGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                      <  2:970657/148‑1 (MQ=255)
                                                               cAGAAAGACAGCTACGGTGTAGCCGGTTGGATGAAGTATATTTATGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                     <  1:661571/149‑1 (MQ=255)
                                                               cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                     >  1:1012914/1‑149 (MQ=255)
                                                                aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGTTAGTAATTAACATCTAAGCGTTAAATACCTGATATCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACACCTTTCTTTTCATCTGGTTGTtta                                                                                    >  2:789274/1‑149 (MQ=255)
                                                                aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                    >  1:464092/1‑149 (MQ=255)
                                                                 gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgt                                                                                          >  1:432990/1‑142 (MQ=255)
                                                                 gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgt                                                                                          <  2:432990/142‑1 (MQ=255)
                                                                 gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACACCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttat                                                                                   <  2:403662/149‑1 (MQ=255)
                                                                  aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                  >  2:593873/1‑149 (MQ=255)
                                                                    aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTACCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                                >  1:216562/1‑149 (MQ=255)
                                                                    aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAAAACCGGATGGCGCGTTGCCATCCCGTAAAATAACATCCCATCAAAGATATTAACCCTTTCTTTTCATCTGGTTCTtaattac                                                                                >  1:707572/1‑148 (MQ=255)
                                                                          ctaCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTCAGCGTTAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                          <  2:540420/147‑1 (MQ=255)
                                                                          cGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTATCTTTTCATCTGGTTGTTTCTTAACCCTTc                                                                          <  1:688291/149‑1 (MQ=255)
                                                                           gACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                       >  2:560941/1‑87 (MQ=38)
                                                                           gACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                       <  1:560941/87‑1 (MQ=38)
                                                                            aCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAg                                                                        <  1:755086/149‑1 (MQ=255)
                                                                               gtgtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCATTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGaa                                                                     >  2:1127734/1‑149 (MQ=255)
                                                                                 gtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATACCATCTAAGATATTAACCCTTTCTTTTAATCTGGTTGTTTAATAACCCTTCAGGAACg                                                                   >  1:1049610/1‑149 (MQ=255)
                                                                                    gCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTc                                                                >  2:421031/1‑149 (MQ=255)
                                                                                      cGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAg                                                              <  2:1106539/149‑1 (MQ=255)
                                                                                        gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGAt                                                            >  1:16207/1‑149 (MQ=255)
                                                                                        gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGAt                                                            <  1:617216/149‑1 (MQ=255)
                                                                                         ttGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGAtt                                                           <  2:303917/149‑1 (MQ=255)
                                                                                           ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGc                                                         <  1:287228/149‑1 (MQ=255)
                                                                                           ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGc                                                         <  2:361085/149‑1 (MQ=255)
                                                                                           ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGc                                                         <  1:297771/149‑1 (MQ=255)
                                                                                            gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCg                                                        >  2:618311/1‑149 (MQ=255)
                                                                                             aTCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGt                                                       <  1:355708/149‑1 (MQ=255)
                                                                                              tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTa                                                      >  1:564494/1‑149 (MQ=255)
                                                                                               cAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTAc                                                     >  2:1090884/1‑149 (MQ=255)
                                                                                                aaGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGATGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACTCTCAGATTGCGTAcc                                                    >  1:93632/1‑149 (MQ=255)
                                                                                                  ggATATGTTTTGTAGTAATTAACATATAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGc                                                  <  2:216562/149‑1 (MQ=255)
                                                                                                  ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGc                                                  >  2:401420/1‑149 (MQ=255)
                                                                                                     tatGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTg                                               <  2:1100408/149‑1 (MQ=255)
                                                                                                      atGATTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAAGCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGc                                              <  2:1182836/149‑1 (MQ=255)
                                                                                                            ggTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAccc                                        <  1:354480/149‑1 (MQ=255)
                                                                                                             gtagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCg                                       >  2:1125999/1‑149 (MQ=255)
                                                                                                             gtagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCACATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCg                                       >  2:614499/1‑149 (MQ=255)
                                                                                                                gtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACACTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGcaa                                    >  1:1158977/1‑147 (MQ=255)
                                                                                                                   aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCg                                 >  2:571302/1‑149 (MQ=255)
                                                                                                                      aaCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                              >  2:57136/1‑149 (MQ=255)
                                                                                                                       aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                             <  2:689722/149‑1 (MQ=255)
                                                                                                                       aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTACAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                             >  2:477421/1‑149 (MQ=255)
                                                                                                                       aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTGTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACACGCCAGCGTTTc                             >  2:863932/1‑149 (MQ=255)
                                                                                                                         aTCTAAGAGTGAAATACCGGAGGGCGAGTTGCCATCCGGTAAAATAACATCCCATGTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTAATAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGa                           <  2:217868/149‑1 (MQ=255)
                                                                                                                          tCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGaa                          <  2:988302/149‑1 (MQ=255)
                                                                                                                            tAAGCGTGAAATACCGGATGGCGAGTTGCCCTCCGGTAAAATAACATCCCAACTAAGATATTAACCCTTTCGTTTCATCTGGTTGTTTATTAACCCTTCCGGAACGCTCAGCTTGCGTACCGCTTGCGATCCCGCCAGCGTTTCGaaaa                        >  2:1074121/1‑147 (MQ=255)
                                                                                                                            tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAata                        >  1:405820/1‑149 (MQ=255)
                                                                                                                                 gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                   >  2:79877/1‑149 (MQ=255)
                                                                                                                                 gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                   <  2:1072315/149‑1 (MQ=255)
                                                                                                                                  tGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatct                  >  1:935528/1‑149 (MQ=255)
                                                                                                                                      aTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTTACCCTTTCTTTTCATCTGGTTGTTTATTTACCCTTCAGGAACGCTCAGATTGAGTACCTCTTGCGAACCCGAAAGCGATTCGAATAttatcttata              >  2:431790/1‑148 (MQ=255)
                                                                                                                                      aTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatcttatc              >  2:942315/1‑149 (MQ=255)
                                                                                                                                         ccGGATGGGGTTTTGCCATCCGGTAAAATAACATCCCCTTTTAGTTATTAACCCTTTCGTTTGATCTGGTTGTTTATTAACCCTTCAGGAACGCGCAGTTTGCGTACCGCTTGCGTACCCGCCAGCGTTTGGAATATTATCTTATCttt           <  2:972805/149‑1 (MQ=255)
                                                                                                                                         ccGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTAAGATTGCGTACCGCTTGCGAACCAGCCAGCGTTTCGAATATTATCTTATCttt           >  1:603514/1‑149 (MQ=255)
                                                                                                                                          cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGGTTGCGAAGCCGCCAGCGTTTCGAATTTTATCTTATCtttt          >  2:374253/1‑148 (MQ=255)
                                                                                                                                          cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTGTTTTACCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATATTATCTTta          >  1:667583/1‑149 (MQ=255)
                                                                                                                                                  gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATc  <  1:620000/149‑1 (MQ=255)
                                                                                                                                                  gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATc  <  1:538936/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATC  >  NC_000913/167145‑167439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑