Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG30.8% 112.1 / 14.4 65G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/41);  new base G (20/0);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 3.75e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATTAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTAC  >  NC_000913/2326207‑2326698
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                             
aaTTAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATggga                                                                                                                                                                                                                                                         >  1:295829/1‑248 (MQ=255)
          gccATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAattacaggtt                                                                                                                                                                                                                                              >  1:6289/1‑243 (MQ=255)
             aTTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAgt                                                                                                                                                                                                                                          >  1:31787/1‑250 (MQ=255)
               tGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAt                                                                                                                                                                                                                                                             <  2:25110/231‑1 (MQ=255)
                      gCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtact                                                                                                                                                                                                                                  <  1:559751/250‑1 (MQ=255)
                        aaaaaTCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttac                                                                                                                                                                                                                               <  1:626346/251‑1 (MQ=255)
                          aaaTCGATTCACAACACGTTGATGAATTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTGAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCTGGTCTTAAAAGTGTGGCGCTTGCGGCCCAGGCCATTCAGGCAGGGCAGGCGCATAGGGTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTc                                                                                                                                                                                                                             <  1:50948/251‑1 (MQ=255)
                                aTTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATTAGGTTAGccccccacctatttgatgc                                                                                                                                                                                                                        >  1:543124/1‑236 (MQ=255)
                                       acgctgtgatgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGGATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCGACTTACTCGATGCAAAAGCAc                                                                                                                                                                                                                <  2:427587/245‑1 (MQ=255)
                                       acacGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAACAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGTGGGGATGGAAAATATTGGTTTTGcccccccctttcaccctgccaaagaac                                                                                                                                                                                                                >  2:707795/1‑229 (MQ=255)
                                        cacGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                                                                                                                                               <  1:549931/251‑1 (MQ=255)
                                          cGTTGATGAAGTGATTATGGGTAACGTGTTACACGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGCTGGGAAAAATGGGTTTTTCCCCCCACTTTCTCGATGCAAAAACACGcc                                                                                                                                                                                                             >  2:293892/1‑250 (MQ=255)
                                              gatgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTg                                                                                                                                                                                                          <  1:4694/250‑1 (MQ=255)
                                               atgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttac                                                                                                                                                                                                                               <  2:557816/228‑1 (MQ=255)
                                               atgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaa                                                                                                                                                                                                                     <  2:22650/238‑1 (MQ=255)
                                               atgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGATGccaa                                                                                                                                                                                                                     >  1:22650/1‑235 (MQ=255)
                                               atgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCtccttacagatcgggagagcacacggctgg                                                                                                                                                                                                        >  1:557816/1‑228 (MQ=255)
                                                   aGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATc                                                                                                                                                                                                    <  1:92449/251‑1 (MQ=255)
                                                       aTTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATc                                                                                                                                                                                                    <  1:478821/247‑1 (MQ=255)
                                                         tATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAGGGGAAAAATGGGTTTAGCCCCCTCCTTTCTCGCTGCAAAAGCAAGCTCTGGTTATCGTTTTg                                                                                                                                                                                             >  2:623465/1‑251 (MQ=255)
                                                            gggTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTAGGCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTgg                                                                                                                                                                                            >  1:500077/1‑249 (MQ=255)
                                                             ggTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTAACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTATCCCCCCACTTAATCGATGCAAAAGCAAGCTCTGGGTTTCGTCCTGgagg                                                                                                                                                                                         >  1:142382/1‑250 (MQ=255)
                                                                        aCAAGCCGGGCTGGGGCAAAATCCGTCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGTAAATATGTGCTTCAGCCCCCTCTTACTCGATTTGCAAGCACGCACTGGTTATTGGCTTGGGGACggcacggctt                                                                                                                                                                             >  1:156722/1‑243 (MQ=255)
                                                                        aCAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGttt                                                                                                                                                                             <  1:589110/251‑1 (MQ=255)
                                                                          aaGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                                                                                                                                           <  1:473065/251‑1 (MQ=255)
                                                                            gCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGa                                                                                                                                                                         <  1:308745/251‑1 (MQ=255)
                                                                             ccGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGAc                                                                                                                                                                        <  2:165183/251‑1 (MQ=255)
                                                                               gggCTGGGGCAAAATCCGGCGCGTCAGGCTCTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATGGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGAGGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGt                                                                                                                                                                      <  2:269343/251‑1 (MQ=255)
                                                                               gggCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGGTGGCATATACGGTGTGCGGATTCACGGGCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGGCTTGGAGACGGACAGGTTTATGACGt                                                                                                                                                                      <  2:128055/251‑1 (MQ=255)
                                                                                      ggCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGGCTTGGAGACGGACAGGTTTATGACGTAATCCTg                                                                                                                                                               <  2:353709/251‑1 (MQ=255)
                                                                                       gCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgc                                                                                                                                                              <  1:122030/251‑1 (MQ=255)
                                                                                       gCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAGCCTgc                                                                                                                                                              <  1:550319/251‑1 (MQ=255)
                                                                                         aaaaTCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGACACGGTGTGCGGCTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGGAATATTAGTTTAGCCCCCTACTTAATCCCTGCACAACGACGCCCTGGGTATTCTCTTTGAGGCGGGCAGGGTTATGACGTAATCCttcgg                                                                                                                                                            >  2:630125/1‑247 (MQ=255)
                                                                                         aaaaTCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATTAGTTTAGGCCCCTCCTTTCTCGATGCAAAAAAACGCCCTGGGTATCGTCTTGGAGACGGGCAGGGTTTTTACGGAATCCTgcgc                                                                                                                                                            >  2:274247/1‑251 (MQ=255)
                                                                                                     gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGCAAAAATGAGGTTTGCCCCCTCCTTTCTCGATTGAAAAACCCGACCTGGTTATCCTCTTTGAGCAGGGAAAGGTTATTACGTAAACCTTCGCGATGGCCTggtg                                                                                                                                               >  1:317231/1‑248 (MQ=255)
                                                                                                                 aaaaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                                                                                                                                  <  1:415862/250‑1 (MQ=255)
                                                                                                                    aaaGAGGGCTGGCAGAAACGGTGTGCGGATTAACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGt                                                                                                                               <  2:53876/251‑1 (MQ=255)
                                                                                                                          ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCTGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCat                                                                                                                         <  1:445941/251‑1 (MQ=255)
                                                                                                                          ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCat                                                                                                                         <  1:485881/251‑1 (MQ=255)
                                                                                                                           gCTGGCAGAAACGGTGTGCGGATTCCCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCata                                                                                                                        <  2:551728/251‑1 (MQ=255)
                                                                                                                            cTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGGCAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCCCCCAGGGTTCTCatat                                                                                                                       >  1:9838/1‑251 (MQ=255)
                                                                                                                                 aGAAACGGTGTGCGGATTCCCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACCCTCTGGTTATCGTCTTGGGGACGGACAGGGTTTTGACGGAATCCTGCGCGATGGGCTGATGTGCGCCACCCATGGGTATAATATggggg                                                                                                                  >  1:161281/1‑250 (MQ=255)
                                                                                                                                        gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgc                                                                                                           <  2:543124/251‑1 (MQ=255)
                                                                                                                                               gATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGAGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATTGGTATTGCCCCCCACTTACTCGCTGGGAAAGAACGCTCCGGGTATCGTCTTTGAGACGGACAAGGTTATTACGGAAACCCGCCCGGAGGCCCGATGGGCGCCCCCCATTGTTATCATTTTGGGGTTAACGCCGaaaa                                                                                                     >  2:48154/1‑250 (MQ=255)
                                                                                                                                                     cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTTGCCCTTGCCGCGCAGGCCAGTCAGGCGGGTCCGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGAAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTGTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATCACCGCCGAAAACGTGGCt                                                                                              <  2:142382/251‑1 (MQ=255)
                                                                                                                                                          aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAaga                                                                                         <  1:623465/251‑1 (MQ=255)
                                                                                                                                                           aTAAAGGATGGGGCTCGGGTCTTCAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGCGCATTGTGGAGGGGGGGAGGGCAAAAATTGGCTTATCCCACTACTTACTCGATGCACAAGCACGCCCCTGTTGTCGTACTTAGGCCGGACAGGTTTCTGACGTAATCCTCCGCGACTGCGCGCTGTGCGCCACCCCTCGTTATAATATTGGGGTTACCGCCGAAAACgggggccaag                                                                                          >  2:732590/1‑240 (MQ=255)
                                                                                                                                                               aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACg                                                                                    <  1:274247/251‑1 (MQ=255)
                                                                                                                                                                 tATGTGGTTCGGGTCTTCAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAGGGAAAAAATGAAGTTAGTCCCCTACTTACCCGATGCAAAAGCACGCTCGGGTTATAGCCTTGGGGACGGGCAGGGTTATGACGTGATCCTGGCCGATGTGCCGATTGGCGCCACCCATGGGTATTATATGGGGGATGCCGCCGAAAACGGGGGTAAagaggaaggg                                                                                  >  1:550549/1‑245 (MQ=255)
                                                                                                                                                                    gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAAAATGGGTTTAGCCCCCTACTTACTCGGTGCAAAAGCAAGCTCTGGTTATCGTCTTGGAGAAGGGCAGGGTTATGAAGTAATCCTGCGCGATGGCCTTGTGTGGGCCACCCATGGGTATCATATTGGGGTTACCGCCGAAAACGTGGCTAAagaggatgggaat                                                                               >  1:287654/1‑242 (MQ=255)
                                                                                                                                                                                 gaaaaGGGTGGCGCTTGCCGCCCAGGCAAGTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTATCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAGGGTTATCATAGGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATg                                                                   <  2:275125/249‑1 (MQ=255)
                                                                                                                                                                                  aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGTGCTTAGCCCCCCACTTTCCCGATGCAAAAGCACGCCCTGGTTATCGTCTTTGAGACCGGCAGGGTTATTACGGAATCCTGCGCGCTGGCCTGATGTGGGCCCCCCCATGGTTTCATATGTGGATTACCCCCCCAAACGTGGCTAAAAAGTACGGGAtttaccgttaaaaggc                                                                 >  2:692782/1‑237 (MQ=255)
                                                                                                                                                                                        gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGTTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATGTGGGGATTACCGCCGAAACCGTGGCTAAAGAGTACGGAATGACCCGTGAAATGCAGGATGa                                                           <  1:147178/251‑1 (MQ=255)
                                                                                                                                                                                          ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCTTATTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTGGGAGACGGACAGGTTTATGACGTAAGCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTTGATAAAGAGTACGGAATTACCCGTGAAATGCAGGATGaa                                                          <  1:736380/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGATTAGGCCCCTACTTACCCGCTGCACAAGCACGCCCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGGGCCCCCCCTGGTTCTCATATGGGGGTTACCGCCGAAAACGTGGCTAAAGAGGACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTaaa                                  >  1:694721/1‑250 (MQ=255)
                                                                                                                                                                                                                              gcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                                                                                                                                           <  1:147391/103‑1 (MQ=255)
                                                                                                                                                                                                                               cAGAGCATTGTGGCGGGGGGTATGGAAAATATGATTTTAGCCCCCTACTTACTCAATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAAt                    <  2:68822/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAAt                    <  1:293892/251‑1 (MQ=255)
                                                                                                                                                                                                                                  agCATTGTGGCGGGGGGTAGGGAAAATATGAGTTTAGCCCCCTTCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGa                 <  2:287654/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAg                <  1:11216/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tgtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCTAACTTCCTTGATGCAAAATCACGCTCTGGTTATCGTCTTGGAGACGGGCAGGTTTATGACGGAATCCTGCGGGATGGCCTGATGGGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTGCATGCACAGCGTAAAGCGGCAGCCGCAATTGAGTcc             <  2:37373/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCAATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCg            <  1:211157/250‑1 (MQ=255)
                                                                                                                                                                                                                                         tgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGAAAAAGCACGCTCTGGTTATCGTCTTGGCGACGGACAGGTTTATGACGTAATCCTGCGCGCTGGCCTGATGTGCGCCACCCATTGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGGGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTAAATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCggg          >  2:620831/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ggCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCCCCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGt          <  1:45005/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCGTGGAGACGGACAGGTTTATGACGTACTCCTGCGCGATGGCCTGATGTGCGCCACACATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTAc  <  2:550549/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                             
AATTAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTAC  >  NC_000913/2326207‑2326698

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