Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,982,5230AG43.8% 222.6 / 43.6 137intergenic (‑136/+31)otsB/araHtrehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (15/62);  new base G (60/0);  total (75/62)
Fisher's exact test for biased strand distribution p-value = 5.68e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

AACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTT  >  NC_000913/1982374‑1982672
                                                                                                                                                       |                                                                                                                                                     
aaCGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaa                                                                                                                                                      <  2:1108903/151‑1 (MQ=255)
 aCGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGaaa                                                                                                                                                     <  2:1345143/151‑1 (MQ=255)
  cGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGgaaa                                                                                                                                                    >  1:148093/1‑151 (MQ=255)
      tCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGgaaaaa                                                                                                                                                  >  2:714927/1‑149 (MQ=255)
         ggcccccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAAAACGCGCTTTACAGACATGTAAGTAAAGGCAGTTTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTa                                                                                                                                             <  2:879294/147‑1 (MQ=255)
              ccGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAGCTaaaatt                                                                                                                                        >  4:1446609/1‑148 (MQ=255)
               cGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGGAAAAAAACTCAAtttt                                                                                                                                       >  2:523350/1‑149 (MQ=255)
                 gTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGGAAAAAAGCTAAAattttg                                                                                                                                     >  1:75921/1‑151 (MQ=255)
                   tctcCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGc                                                                                                                                    <  2:589827/150‑1 (MQ=255)
                   tctcCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCAGTTGGCTATGGTGGGAAAAAACGCTAAACTGTTGc                                                                                                                                    <  2:761800/150‑1 (MQ=255)
                       ctcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaa                                                                                                                               <  1:308681/151‑1 (MQ=255)
                        tcctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaa                                                                                                                              <  3:229110/151‑1 (MQ=255)
                         cctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa                                                                                                                             <  4:605367/151‑1 (MQ=255)
                         cctTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaa                                                                                                                               <  2:961621/149‑1 (MQ=255)
                           ttCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTTTTGGGGGGAAAAAAGCTAAATTTTTGCAGaaaaaa                                                                                                                            >  2:1187258/1‑150 (MQ=255)
                            tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaaa                                                                                                                            <  2:924996/149‑1 (MQ=255)
                            tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGGATTGGGGGGAAAAAACCCAAATTTTTTGGGAaaaaaa                                                                                                                            >  4:96937/1‑150 (MQ=255)
                             ctctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAGCTAAATTGTTGCAGAAAAAAGc                                                                                                                          >  2:421076/1‑150 (MQ=255)
                                 tttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTTTTGGGGGGAAAAAAGCCAAAATTTTTGAaaaaaaaaaatca                                                                                                                       >  2:477761/1‑145 (MQ=255)
                                    cTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCCTTGGGGGGGAAAAAAGCCAAATTTTTGTAGAAAAAAAAATCAGAc                                                                                                                   >  2:520653/1‑150 (MQ=255)
                                      tATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTTGGGATGGGGGGGAAAAAAGCGAAAATTTTGCGGAAAAAAACAACAGAAAgg                                                                                                                >  4:1250204/1‑150 (MQ=255)
                                        ttGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTATGGGGGGGAAAAAAACCAAAATTTTTCAGaaaaaaaaaacaaacaatgg                                                                                                              >  2:932914/1‑137 (MQ=255)
                                                    gAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAATTGTTTGAGGAAAAAACATCACACAAGGGGGTTTGCttttt                                                                                                  >  3:386904/1‑151 (MQ=255)
                                                     aaCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTg                                                                                                 <  2:866126/151‑1 (MQ=255)
                                                     aaCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTg                                                                                                 <  2:430319/151‑1 (MQ=255)
                                                     aaCGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTg                                                                                                 <  2:299427/151‑1 (MQ=255)
                                                         cgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTg                                                                                             <  2:1057056/151‑1 (MQ=255)
                                                           cTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAAATGTTGCAGAAAAAAACATCAGACAGTGCGGTTCGCTTTTTTCTTGtt                                                                                           >  4:251056/1‑150 (MQ=255)
                                                           cTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAAATGTTGCAGAAAAAAACAACAGACAGTGCGGTTTGCTTTTTTCTTGTa                                                                                           >  1:419965/1‑151 (MQ=255)
                                                           cTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAACATAAAATGTTTCAGAAAAAAAAAACAAACAATGAGTTTTGCTTTTTTTTTGtt                                                                                           >  3:40086/1‑150 (MQ=255)
                                                           cTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTTTGGTGGGGAAAAAAGCTAAATTTTTTCAGAAAAAAACATCAGACAGTGCGGTTCGCTTTTTTCttttt                                                                                           >  1:363264/1‑148 (MQ=255)
                                                            tttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACAACAGAAAGTGCGGTTTGCttttttgtttta                                                                                           >  2:1148898/1‑143 (MQ=255)
                                                             ttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTATGGGGGGGAAAAACGCTAAAATTTTGCAGAAAAAAAAAACAGACAGGGGGGTTCGCTTTTTTCTTGTaaa                                                                                         >  1:237544/1‑150 (MQ=255)
                                                             ttACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAACGCTAAAATGTTGCAGAAAAAAACATCAGACAGGGCGGTTTGCTTTTTTTTTGtta                                                                                          >  2:1354250/1‑148 (MQ=255)
                                                              tACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAACGCTAAAATTTTGTAGAAAAAAAAATCAAACAAGGCGGTTCGCTTTTTTTTTGTaa                                                                                          >  4:924574/1‑149 (MQ=255)
                                                               aCAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACg                                                                                        <  4:1082036/150‑1 (MQ=255)
                                                                 aGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAACTAAAATTTTTTAGAAAAAAAAATTAAACAGTGCGGTTTGCTTTTTTCTTGTAAAGGtg                                                                                     >  4:901633/1‑150 (MQ=255)
                                                                    cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTATGGGGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACATCAGAAAATGCGTTTTGCTTTTTTCTTTTAAAGGGCGaa                                                                                  >  1:167470/1‑151 (MQ=255)
                                                                    cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATTGGGGGGAAAAAAGCGAAATTGTTTCAGAAAAAAACATCAGACAGTGCGGTTTGCTTTTTGCTTGTTACGGGCgga                                                                                  >  1:203350/1‑149 (MQ=255)
                                                                     aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAACTAAAATGTTGCAGAAAAAAACAACAAACAGTGCGGTTTGCTTTTTTCTTTTAAAGGGCGAAg                                                                                 >  1:17108/1‑151 (MQ=255)
                                                                     aTGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTTTGGTGGGGAAAAAAGCTAAATTGTTGTAGGAAAAAAAATCAGAAAGGGCGGTTCGCCTTTTTTTTGTTAAGGTCGa                                                                                   >  4:1385291/1‑149 (MQ=255)
                                                                      tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACATCAGACAGTGCGGTTTGCTTTTTTTTTGTAACGGTCGAAGa                                                                                >  1:436332/1‑151 (MQ=255)
                                                                        tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATc                                                                              <  1:554361/151‑1 (MQ=255)
                                                                         aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCa                                                                             <  4:858986/151‑1 (MQ=255)
                                                                         aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGGATTGTGGGGAAAAAAACTAAAATGTTGCAGAAAAAAAAAACAAACAAAGCGGTTTGCTTTTTGTTTTTTAAGGGCGaaaaaca                                                                             >  3:226815/1‑146 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTTGGGGGGGAAAAAAACCAAATTTTTGCAGAAAAAAACATCAGAAAGTGGGGTTTGCTTTTTGCTTTTAACGGTCGAAGAACa                                                                             >  3:291333/1‑150 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAACACTAAAATGTTGGAGAAAAAAAAATCAGACAGGGGGGTTTGCTTTTTTCTTTTTAAGGGCGAAGAtaa                                                                             >  3:198995/1‑148 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACATCAGAAAGTGCGGTTTGCTTTTTTCTTTTAACGGGCGAAGATCa                                                                             >  1:516956/1‑150 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAGCTAAATTTTTGCAGAAAAAAACATCAGACAGTGCGGTTCGCTTTTTTCTTGTAAAGGTCGAAGATCa                                                                             >  1:444474/1‑150 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGGATGGGGGGGAAAAAAAATAAATTGTTGCAGAAAAAAACATCAGACAAGGGGGTTGGCTTTTTGTTTTTTACGGGCGaaaaacaa                                                                            >  3:423274/1‑145 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                                            <  2:146294/151‑1 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGGGGGGAAAAAAGCTAAATTGTTGCAGAAAAAAAAATCAGACAGTGCGTTTTGCTTTTTTTTTGTAAAGGGCGAAGATCa                                                                             >  1:78964/1‑150 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTCGCTATGGGGGGGAAAAAAGCTAAATTTTTGGAGAAAAAAAAATCAGACAGAGCGTTTTGCTTTTTTTTTTTAAAGGGCGAAGATCa                                                                             >  1:196756/1‑150 (MQ=255)
                                                                          aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTGTGGTGGGGAAAAAAACAAAAATTGTGGAGAAAAAAACACCAAAAAGGGGGGTTTGCTTTTTTTTTGTTAAGGTCGAAGAACa                                                                             >  3:29769/1‑150 (MQ=255)
                                                                           gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAACCAAAATTTTGCAGAAAAAAACAACAAACAGTGCGGTTTGCTTTTTTCTTTTAAAGGGCGAAGAacac                                                                            >  1:496517/1‑150 (MQ=255)
                                                                             ataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGTGGGGAAAAACGCTAAATTGTTTCAGAAAAAAACATCAGACAATGCGGTTCGCTTTTTTCTTGTTACGGTCGAAGATCACCGc                                                                         >  2:701647/1‑151 (MQ=255)
                                                                              taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGTGGGGAAAAAAGCTAAATTGTTTCAGAAAAAAACATCAGACAGTGCGTTTCGCTTTTTTCTTGTAACGGGCGAAGATAACCGCt                                                                        >  2:155757/1‑151 (MQ=255)
                                                                                 aCAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATTGGGGGGAAAAACGCGAAAATGTTGTAGAAAAAAACAACAGACACAGCGCTTTGCCTTTTTTTTGTTAAGGGCGAAGATTACCGCTGcc                                                                     >  2:227819/1‑151 (MQ=255)
                                                                                 aCAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGcc                                                                     <  4:1127050/151‑1 (MQ=255)
                                                                                    gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTTTTGGGGGGAAAAACGCTCAAATTTTTCAGAAAAAAACAACAGACAATGGGTTTTGCTTTTTTCTTTTAAACGGCGAAAATCACCCccgcccgc                                                                  >  2:939801/1‑151 (MQ=255)
                                                                                    gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGGGGGGAAAAAAGCTAAATTGTTGCAGAAAAAAAAAACAGACAGTGCGGTTTGGTTTTTTCTTTTAAAGGTCGGAGATCACCGCTGcccg                                                                   >  1:71907/1‑148 (MQ=255)
                                                                                        gACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAACCAAAATTTTGCAGAAAAAAAAAACAAAAAGTGCGGTTTGCTTTTTTCTTTTTACGGGCGGAAATAACCCCTGCCCGCAggg                                                              >  4:186478/1‑150 (MQ=255)
                                                                                                 cGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcg                                                      <  2:96123/150‑1 (MQ=255)
                                                                                                     aTTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATTGTGGGGAAAAAAAATAAAAATTTGGCGAAAAAAACATCAGAAAGGGCGTTTTTCCTTTTTGTTTTAAAGGTCCAAAAGCCACCCTGCCCGGAGGATTAAGcctccga                                                  >  4:1279798/1‑145 (MQ=255)
                                                                                                        tGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCa                                               <  2:1446606/150‑1 (MQ=255)
                                                                                                          cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGt                                            <  4:900165/151‑1 (MQ=255)
                                                                                                          cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGt                                            <  4:88380/151‑1 (MQ=255)
                                                                                                          cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGt                                            <  3:148092/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGGTTTGGTGGGAAAAAAAACTAAATTTTTGGCACAAAAAGCAATAAACAGGGCCATTCCCTTTTTGCTTGTCATTGTTGTAAATCACCGATGCCCGCTGGTTTTAGAcgtgcacaaaatcc                                          >  2:1464573/1‑139 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  1:42180/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  4:589028/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  1:244909/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  3:553141/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  3:501044/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  2:1414187/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  3:248216/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  3:320506/151‑1 (MQ=255)
                                                                                                            aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc                                          <  1:563204/151‑1 (MQ=255)
                                                                                                              gaagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTg                                        <  4:1063049/151‑1 (MQ=255)
                                                                                                              gaagaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTg                                        <  2:1268067/151‑1 (MQ=255)
                                                                                                                 gaaCAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAACGCTAAATTGTTGCAGAAAAAAAAATCAGACAGTGGGTTTTGCTTTTTTCTTGTAACGGTCGAAGAACACCGCTGCCAGCAGGGTTAAACCGCGAACCACGGACTgcg                                      >  1:362436/1‑150 (MQ=255)
                                                                                                                 gaaCAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTgcgc                                     <  4:534302/151‑1 (MQ=255)
                                                                                                                      gCCAACCGTTATGGCGGTTGGCTATGGGGGGGAAAAACACCAAAATTTTGCAGAAAAAAAAATCAGACAGTGCGGTTTGCTTTTTTCTTTTAACGGTCGAAGAACACCGCCGCCCGCAGGGTTAAACCGCGGAACACGGACTGCGCGAAAg                                >  3:559325/1‑151 (MQ=255)
                                                                                                                         aaCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                             <  4:1187266/151‑1 (MQ=255)
                                                                                                                         aaCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag                             <  2:86729/151‑1 (MQ=255)
                                                                                                                           ccGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaa                           <  4:421082/151‑1 (MQ=255)
                                                                                                                                  ggCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTaa                    <  4:520657/151‑1 (MQ=255)
                                                                                                                                  ggCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTa                     <  2:387233/150‑1 (MQ=255)
                                                                                                                                    cGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGc                  <  3:363259/151‑1 (MQ=255)
                                                                                                                                      gTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGc                  >  4:1444834/1‑149 (MQ=255)
                                                                                                                                          ggCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt             >  2:300426/1‑151 (MQ=255)
                                                                                                                                             tATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCAt          >  4:1113734/1‑151 (MQ=255)
                                                                                                                                                tggtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATgg        <  3:172608/151‑1 (MQ=255)
                                                                                                                                                 ggtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATgg        <  4:1352677/150‑1 (MQ=255)
                                                                                                                                                 ggtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGc       <  2:392013/151‑1 (MQ=255)
                                                                                                                                                  gtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCg      >  3:424292/1‑151 (MQ=255)
                                                                                                                                                   tggGAAAAAACGCTAAATTGTTGCAGAAAAAGGCATCAGACAGGGCGTTTTGCTTTTTGTTTGAAACGGCCGAAGATCACCGCTGCCAGCGGGATTAAGCGGCGAACCCCGAACGGGGCGAAAGGAGAATTATTAAGCAGGTTCAGGGCGt     >  2:248352/1‑151 (MQ=255)
                                                                                                                                                   tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGt     >  2:598024/1‑151 (MQ=255)
                                                                                                                                                   tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGt     >  2:1438035/1‑151 (MQ=255)
                                                                                                                                                      gAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGtttt  >  4:571212/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AACGGTTCTGTCACCCGGTTCTCCTCCTTCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTT  >  NC_000913/1982374‑1982672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: