Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG43.5% 150.6 / 13.4 86G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (6/42);  new base G (37/0);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.21e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326304‑2326598
                                                                                                                                                  |                                                                                                                                                     
gcgcGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                                                                                   <  2:1044151/151‑1 (MQ=255)
gcgcGTCAGGCACTGTTAAAAAGCGGCGTTGCAGAAACGGTGGGCGGATTCACGGTCAATAAAGTAGGTGGTTCGGGTCTTAAAACTGTGGGGGTTGCCGCCCAGACCATTCAGGCAGGTCAGGCTCAGAGCATTGTGGTGGGGAGTAggg                                                                                                                                                   <  2:319020/151‑1 (MQ=255)
    gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                               <  1:404521/151‑1 (MQ=255)
    gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                               <  1:298912/151‑1 (MQ=255)
      cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                                              <  1:1008022/150‑1 (MQ=255)
       aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAatat                                                                                                                                            <  2:1411355/151‑1 (MQ=255)
       aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAatat                                                                                                                                            <  2:619024/151‑1 (MQ=255)
       aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAatat                                                                                                                                            <  1:203592/151‑1 (MQ=255)
        ggCACTGTTAAAAAGCGTGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATGAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTGGCCGCCCAGGCCATTCAGTCAGGTCAGGTGCAGAGCATGGTGACGGGGGGTTTGGaaa                                                                                                                                                <  2:206366/146‑1 (MQ=255)
        ggCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCAATCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATg                                                                                                                                           <  2:344662/151‑1 (MQ=255)
         gCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAata                                                                                                                                             >  1:684987/1‑148 (MQ=255)
          cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTATTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGaaaa                                                                                                                                               >  2:1209132/1‑145 (MQ=255)
            cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATggggt                                                                                                                                       >  2:917082/1‑147 (MQ=255)
              gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                                                                                      <  1:559366/150‑1 (MQ=255)
              gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                     <  2:161152/151‑1 (MQ=255)
              gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                     <  2:1333047/151‑1 (MQ=255)
                 aaaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATAAGAGGTTAGc                                                                                                                                   >  2:1354550/1‑150 (MQ=255)
                    aaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGccccc                                                                                                                               >  1:1110686/1‑151 (MQ=255)
                       cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGTTTAGcccccca                                                                                                                             >  2:1260609/1‑148 (MQ=255)
                        gggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtact                                                                                                                           <  1:574635/151‑1 (MQ=255)
                            tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGTGTTTATGCCCCgacttac                                                                                                                        >  2:1134291/1‑150 (MQ=255)
                             ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                                       <  1:44811/150‑1 (MQ=255)
                                tgAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATTAGGTTAGCCCCCCACTTTCTCGAt                                                                                                                   >  2:257369/2‑151 (MQ=255)
                                aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGAt                                                                                                                   <  1:994888/151‑1 (MQ=255)
                                  aaaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCtacttac                                                                                                                        >  2:720773/1‑144 (MQ=255)
                                  aaaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGGGGTTTGCCCCCCACTTTCTCGagg                                                                                                                  >  2:1396705/1‑148 (MQ=255)
                                           gCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                                        <  1:1268003/151‑1 (MQ=255)
                                           gCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTCCTTACTCGATGGaaaa                                                                                                             >  1:40934/1‑146 (MQ=255)
                                             ggATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGct                                                                                                      <  1:959393/151‑1 (MQ=255)
                                              gATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTTCTCGATGCAAAAGCACg                                                                                                        >  1:1223546/1‑148 (MQ=255)
                                               aTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                                        >  1:453656/1‑147 (MQ=255)
                                                   aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGtt                                                                                                <  2:922475/151‑1 (MQ=255)
                                                   aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGtt                                                                                                >  1:861296/1‑151 (MQ=255)
                                                    cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGtt                                                                                                <  2:1504901/150‑1 (MQ=255)
                                                    cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAACAAATGAGTTTATCCCCCCACTTACCCGATGCAAAAACACGCCCTGGtt                                                                                                >  1:1141526/1‑150 (MQ=255)
                                                     ggTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTAt                                                                                              >  1:291432/1‑151 (MQ=255)
                                                       tCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCCCTgggttt                                                                                              >  2:791478/1‑145 (MQ=255)
                                                         aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTc                                                                                         <  1:1232429/151‑1 (MQ=255)
                                                         aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGt                                                                                          >  1:1235845/1‑150 (MQ=255)
                                                         aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGGGGTTAGCCCCCTACTTAATCGATGCAAAAGCACGCTCTGGTTATCGt                                                                                          >  1:1438645/1‑150 (MQ=255)
                                                            aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGGAAAAATTAGTTTAGCCCCCAACTTTCTAGGTGTAAAAGGACGGTCTGGTGATCgcatt                                                                                        >  2:83526/1‑146 (MQ=255)
                                                            aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                       >  1:1409404/1‑150 (MQ=255)
                                                            aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCtt                                                                                       >  1:920919/1‑150 (MQ=255)
                                                              aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGga                                                                                    <  1:1261510/151‑1 (MQ=255)
                                                                tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATTGtttt                                                                                       >  2:778489/1‑146 (MQ=255)
                                                                tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTgg                                                                                     >  1:409026/1‑148 (MQ=255)
                                                                  tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTgaga                                                                                  >  1:1320518/1‑149 (MQ=255)
                                                                   gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACg                                                                                <  2:1377799/150‑1 (MQ=255)
                                                                    tgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTTTCGTCTTGGAGACgg                                                                               >  1:1162753/1‑150 (MQ=255)
                                                                     ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGa                                                                              <  2:1430421/150‑1 (MQ=255)
                                                                     ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTGCCCCCTACTTACTCGATGCAAAAGCACGCTCTTGTTATCGTCTTTGGGACgg                                                                               >  1:1411194/1‑149 (MQ=255)
                                                                      gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGAGTTTTGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGGGACgg                                                                               >  1:1295054/1‑148 (MQ=255)
                                                                       ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGCATGGAAAAAATGAGCTTAGACCCCCCCTTTTTCGAGGCAGAAACACGGCCTGGGTTTCGTTCTTgggacg                                                                                 >  2:203777/1‑142 (MQ=255)
                                                                        tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg                                                                           <  1:176949/150‑1 (MQ=255)
                                                                         cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGCGTTTAGCCCCCTACTTACTCGATGCAAAAACACCGCCTGGGTATCGTGTTGGAGACgggacag                                                                          >  2:1069260/1‑145 (MQ=255)
                                                                               ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTTTCGTCTTGGAGACGGACAGGGTTAt                                                                     >  1:1406848/1‑149 (MQ=255)
                                                                                tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                                     <  1:149963/148‑1 (MQ=255)
                                                                                tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGa                                                                   <  1:998404/150‑1 (MQ=255)
                                                                                 aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACg                                                                 <  1:561493/151‑1 (MQ=255)
                                                                                 aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCCCTGGGTATCTTCTTGGAGACGGACAGGGTTTAGACg                                                                 >  2:1093807/1‑151 (MQ=255)
                                                                                    aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTaa                                                              <  1:1355513/151‑1 (MQ=255)
                                                                                    aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTa                                                               >  2:1032836/1‑150 (MQ=255)
                                                                                      tgtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACAAGCCCTGGTTATCGTCTTGGGGACGGACAGGTTTATGACGTaa                                                              >  1:780992/1‑149 (MQ=255)
                                                                                       gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATcc                                                           <  1:501401/151‑1 (MQ=255)
                                                                                         ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTg                                                         <  1:990545/151‑1 (MQ=255)
                                                                                           cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAAGGGAAATATGAGTTTAGCCCCCTTCTTACTCGATGCCAAAGCACGCTCTGGTTATCGTTTTGGAGACGGACAGGGTTATGACGGAATTCTgcg                                                       >  1:1200865/1‑151 (MQ=255)
                                                                                            gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                      <  1:63460/151‑1 (MQ=255)
                                                                                               tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTTTCGTCTTGGAGACGGAAAGGGTTATTACGTAATTCTGcgcg                                                     >  1:617378/1‑149 (MQ=255)
                                                                                                  cgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGc                                                <  1:1164357/151‑1 (MQ=255)
                                                                                                  cgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGAGGGGTTGGGAAATTTTAGTTTAGCCCCCCACTTAATCGATGCAAAAACAGGATCTCGTTTTAGCGTTGGAGAAAGACAGGGTTATTACGTAAtaccacg                                                       >  2:1432060/1‑138 (MQ=255)
                                                                                                    cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGcc                                               >  1:991217/1‑150 (MQ=255)
                                                                                                        ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGGGATGGCCTGAtg                                          <  1:920721/151‑1 (MQ=255)
                                                                                                               caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGcc                                    >  1:1292156/1‑150 (MQ=255)
                                                                                                               caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGcc                                    >  1:333237/1‑150 (MQ=255)
                                                                                                                    aggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGATAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                              <  1:94130/151‑1 (MQ=255)
                                                                                                                    aggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                              <  1:1050176/151‑1 (MQ=255)
                                                                                                                      gTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATgg                            <  2:610389/151‑1 (MQ=255)
                                                                                                                        cAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGATTTGCCCCCCACTTTTTCGATGCAAAAAAAAGCTCTGGGTTTCGTCTTGGAGACGGACCGGTTTATTACGTAAATCTGCGCGATGGCCTGATGTGGGCCACCCCTGGtt                          >  1:429023/1‑151 (MQ=255)
                                                                                                                          ggCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTAt                        <  1:223608/151‑1 (MQ=255)
                                                                                                                                       gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTa           <  1:927709/151‑1 (MQ=255)
                                                                                                                                       gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTa           <  1:1126514/151‑1 (MQ=255)
                                                                                                                                        tgGCGGGGGGTATGTAAAATATGAGTTTAGCCCGCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGATACTGACAGGTTTATGACGTAATCCTGCGCGATGGGCTGATGTGCGCGACCCATGGTTAGCAGATGGGGAtt            <  2:146220/149‑1 (MQ=255)
                                                                                                                                        tgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAc          >  2:584355/1‑151 (MQ=255)
                                                                                                                                          gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccg        >  1:692400/1‑151 (MQ=255)
                                                                                                                                           cGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgc       >  1:1233082/1‑151 (MQ=255)
                                                                                                                                                ggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaa  <  1:154146/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326304‑2326598

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