Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG37.5% 89.7 / 13.8 56G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/32);  new base G (21/0);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.50e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGTAATTAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACA  >  NC_000913/2326202‑2326699
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
aaaGTAATTAAAGCGGCCATTGAACGGGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGAGTTACACGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGGGCAGAGCATTGTGTCGCGGGGTAt                                                                                                                                                                                                                                                             <  2:42555/249‑1 (MQ=255)
  aGTAATTAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                                                                                                                                           >  1:63176/1‑251 (MQ=255)
       ttAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaa                                                                                                                                                                                                                                                       >  1:54007/1‑250 (MQ=255)
        tAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                                                                                                                                                      <  1:394268/250‑1 (MQ=255)
        tAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGaaaaa                                                                                                                                                                                                                                                      >  2:203369/1‑249 (MQ=255)
         aaaGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAata                                                                                                                                                                                                                                                     <  1:258101/250‑1 (MQ=255)
               gccATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                                                                                                                                                                                              <  1:191054/251‑1 (MQ=255)
                 cATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGttt                                                                                                                                                                                                                                              >  1:164523/1‑249 (MQ=255)
                     gAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                                                                                                                                                         <  1:18472/250‑1 (MQ=255)
                         gTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAGGGGAAATATGGGTTTTGACCCCtac                                                                                                                                                                                                                                    >  2:203371/1‑251 (MQ=255)
                           gCAAAAATCGATTCACAACACGGTGATGAATTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGACATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtact                                                                                                                                                                                                                                   <  1:242715/250‑1 (MQ=255)
                             aaaaaTCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAGccccccacctac                                                                                                                                                                                                                                >  1:146321/1‑244 (MQ=255)
                                     aTTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATCTGAGTTTAGCCCCCtcccttctcggtgcc                                                                                                                                                                                                                        >  1:363350/1‑237 (MQ=255)
                                       tCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTTGcccccacctttttcgattcca                                                                                                                                                                                                                       >  2:433459/1‑234 (MQ=255)
                                               cGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGct                                                                                                                                                                                                             <  2:17980/251‑1 (MQ=255)
                                                    atgaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGt                                                                                                                                                                                                        <  1:98901/251‑1 (MQ=255)
                                                     tgaAGGGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAAACCGGCGCGTCAGGCACTGTTAAAGAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTAGGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCATGTCAGGCGCAGAGCATTGGGGTGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCGCTGGt                                                                                                                                                                                                        <  1:168632/250‑1 (MQ=255)
                                                      gaAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGAAAATATGGGGTTAGCCCCCCACTTACTCCGTTCAAAAACACGCTCTgggtt                                                                                                                                                                                                      >  2:112614/1‑248 (MQ=255)
                                                        aGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                                                                                                                                               <  2:29368/225‑1 (MQ=255)
                                                        aGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAAATGGGTTTAGTCCCCTACTTACTCGATCTAAAAACACGCTCTGGttttc                                                                                                                                                                                                    >  1:129806/1‑248 (MQ=255)
                                                        aGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAGGTTTGCCCCCtacttact                                                                                                                                                                                                                               >  1:29368/1‑225 (MQ=255)
                                                          ggATTATGGGTAACGTGTTACAAGACGGGCTGGGGCAAAAGCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAACCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCCGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGt                                                                                                                                                                                                  <  1:338044/250‑1 (MQ=255)
                                                            aTTATGGGTAACGTGGTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGCGGTTCGGGTCGTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGTGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                                                                                                                                                <  1:340799/251‑1 (MQ=255)
                                                                   gTAACATGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGAc                                                                                                                                                                                         <  1:314386/251‑1 (MQ=255)
                                                                     aaCGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACgg                                                                                                                                                                                       <  1:96300/251‑1 (MQ=255)
                                                                                 gCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGACAGGGTTATg                                                                                                                                                                           >  1:214586/1‑250 (MQ=255)
                                                                                          gggCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAGATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCt                                                                                                                                                                 <  1:61304/251‑1 (MQ=255)
                                                                                           ggCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTg                                                                                                                                                                <  2:358305/251‑1 (MQ=255)
                                                                                              aaaaTCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGGAAATATGGGGTTATCCACCTATTTACTTCTATCAAAAAAACGCCCTGGTTATGCTCTTGGGGATGGACAGGTTTAttacagtatgctgcgc                                                                                                                                                             >  1:37771/1‑236 (MQ=255)
                                                                                               aaaTCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                                                                                                                                           <  1:392204/251‑1 (MQ=255)
                                                                                                              ggCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAAGGAAAATATGAGTTTTGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATTGTCTTTGAGACGGGCAGGTTTATGACGTAATCCTGCGCGATGGCCTGCtgtggg                                                                                                                                           >  1:367119/1‑249 (MQ=255)
                                                                                                                    gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTGGCCCCTACTTACTTGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGAAAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAccc                                                                                                                                     >  1:443169/1‑251 (MQ=255)
                                                                                                                       aaaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGCGGGGATGGGAAATATGAGTTTAGACCCCTACTTACTTGATGCAAAAACACGCTCTGGTTATCGTCTTGGGGGCGGGCAGGGTTATGACGTAATCCCTCGCGGTGGCCTGAATGGCGCCACCCa                                                                                                                                    >  1:120418/1‑249 (MQ=255)
                                                                                                                         aaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTGCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGt                                                                                                                                <  1:203369/251‑1 (MQ=255)
                                                                                                                                           cGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAcc                                                                                                              <  1:82473/251‑1 (MQ=255)
                                                                                                                                               gtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCCTGGTTATCATATGGGGGTTAccgccg                                                                                                          >  1:296344/1‑251 (MQ=255)
                                                                                                                                                   ggATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGt                                                                                                                                                                       <  1:421898/188‑1 (MQ=255)
                                                                                                                                                       tCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAACCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGACGGACAGGGTTATGACGTAATCCTGCGCGGTGGCCTGCTGGGGGCCACCCACGGTTATTATATGGGGGGTAACGCCGCAAACggg                                                                                                  >  2:454664/1‑249 (MQ=255)
                                                                                                                                                       tCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGGCTTGGAGACGGAAAGGGTTATGACGGAAACCTGCGCGATGGCCCGGTGTGCGCCCCCCAAGGGTTTCATATGGGGATTACCGCCCAAAACggg                                                                                                  >  1:236989/1‑249 (MQ=255)
                                                                                                                                                                    aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGCGGGGAAAAATGAGTTTAGCCCCCCACTTAATCGATTCAAAAGCACGCCCTGGTTTTCCGCTTTGAGACGGGCAAGGTTTAGACGTAATACTTCGCGATGGGCTGTTGGGTGCCACCCATGGGTTTCATATGGGGATTACCCCCCAAAAcaggggccaaaaggtcg                                                                                     >  2:320660/1‑235 (MQ=255)
                                                                                                                                                                      tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGAAAATATGAGTTTTGCCCCCCACTTACTCGATGCAAAAAGACGCGCTGGTTATCGTCTTTGGGAAGGAAAGGTTTATGACGTAATCCTCCGCGGTGGCCCGAAGTGCGCCACCCCTGGTTATCATATGGGGGTTACCGCCGAGAAAGTGGGTAAAGAGTACgg                                                                                    >  1:246899/1‑250 (MQ=255)
                                                                                                                                                                         gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAAt                                                                                 <  1:186684/250‑1 (MQ=255)
                                                                                                                                                                                gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATGCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGGCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTg                                                                         <  2:329330/251‑1 (MQ=255)
                                                                                                                                                                                 ggTCTTAACAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTTCGCGGTGGCCTGATGTGCGCCACCCATGGCTTTTATATGGGGATTAACGCCGAAAAAGGGGGTTAAGAGTACGGAATTACCCGTg                                                                         >  1:108560/1‑250 (MQ=255)
                                                                                                                                                                                       aaaaGTGGGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGTTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGAATGACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGc                                                                   <  2:129806/250‑1 (MQ=255)
                                                                                                                                                                                                     tGCCGCCCAGGCGATTCAGGCAGGTCAGGCGCAGAGCATTTTGGCGGGGGGTAGGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACCCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGTATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGgcg                                                     <  2:435357/250‑1 (MQ=255)
                                                                                                                                                                                                           ccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACCCTCTGGTTATCGGCTTGGAGACGGGCAGGTTTATGACGTAATCCTGCGCGATGGGCTGGTGTGCGCCCCCCATGGGTTTCATATGGGGGTTACCCCCCGAAACGGGGCTAAAGAGTACGGAAATACCCGTGAAATGCAGGATGAAATGGCGCTACAtt                                              >  1:89945/1‑251 (MQ=255)
                                                                                                                                                                                                                  aTTCAGGCAGGTCAGGCGCAGAGCATTGGGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCg                                       <  2:63176/251‑1 (MQ=255)
                                                                                                                                                                                                                     caggcaggTGAGGCGGAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCAACTTACTCGATCCAAAAGCACGCGCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGGAACCGGGGTCGAGGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGGTTACCGCCGAAAACGTGGCTATAGAGTACGGAATTACCCGTGATATGCAGGATGAACTGGCTCTACCTTCACAGCGtc                                     <  2:139374/250‑2 (MQ=255)
                                                                                                                                                                                                                     caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTaa                                    <  1:112614/251‑1 (MQ=255)
                                                                                                                                                                                                                      aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTaaa                                   <  1:172483/251‑1 (MQ=255)
                                                                                                                                                                                                                        gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGGTATCATATGGGGGTTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGAAGAAATGGCGCTACATTTACAGCGTaaaac                                 >  1:292536/1‑249 (MQ=255)
                                                                                                                                                                                                                              cAGGCGAAGAGCATGGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGc                           <  2:246899/251‑1 (MQ=255)
                                                                                                                                                                                                                                       agCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGa                  <  2:363350/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTc               <  1:23419/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gTATGGAAAATAGGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCGCTGGTTATCGTCGTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGTTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATGCACAGCGGAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACa  <  2:146321/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
ACAGTAATTAAAGCCGCCATTGAACGTGCAAAAATCGATTCACAACACGTTGATGAAGTGATTATGGGTAACGTGTTACAAGCCGGGCTGGGGCAAAATCCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAACGTGGCTAAAGAGTACGGAATTACCCGTGAAATGCAGGATGAACTGGCGCTACATTCACAGCGTAAAGCGGCAGCCGCAATTGAGTCCGGTGCTTTTACA  >  NC_000913/2326202‑2326699

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: