Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,0000TG35.2% 64.1 / 10.2 54Y557S (TAC→TCC) maeAmalate dehydrogenase (oxaloacetate‑decarboxylating)
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base G (19/0);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAGGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGAAGCAAT  >  NC_000913/1553753‑1554247
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
agcagGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTAtt                                                                                                                                                                                                                                                          <  2:228843/251‑1 (MQ=255)
  cagGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCg                                                                                                                                                                                                                                                        <  1:407018/251‑1 (MQ=255)
           gCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGcc                                                                                                                                                                                                                                                 >  1:551413/1‑249 (MQ=255)
            ccATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTGcccg                                                                                                                                                                                                                                               >  2:98331/1‑248 (MQ=255)
                      gTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTc                                                                                                                                                                                                                                   <  1:98331/251‑1 (MQ=255)
                                      aaggatcatcatgaCATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAATTATCGTCAATGGCCTGTTGCAgg                                                                                                                                                                                                                  >  1:295582/14‑251 (MQ=255)
                                       taCGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGGTTGCCAGAAATTATCGTCAATTGCCTGTTTCAggg                                                                                                                                                                                                                 >  1:427740/1‑251 (MQ=255)
                                          gTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTGCCAGAAATTATCGACAATGGCCTGTTGCAGGGCtt                                                                                                                                                                                                              >  1:513864/1‑251 (MQ=255)
                                               gTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGGTTGCCCGAAATTATCGGCAATGGCCTGTTGGAGGGGTTCggtg                                                                                                                                                                                                         >  2:21980/1‑249 (MQ=255)
                                                  aaCATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCATTTTCGGATGATGTTCTGCCTAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCCCGGGCTTTGGCTTGACAGAAATTATAGTCAATGGCCCGGTGCAGGGCtttgggaggg                                                                                                                                                                                                      >  1:276668/1‑243 (MQ=255)
                                                        atgatCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTc                                                                                                                                                                                                <  1:410278/251‑1 (MQ=255)
                                                                    aacaacTGCAGCCTTCATAGTTCCTCCTTTTCGGATGCTGTTCTGCATAGCAGGTGACGCAAATGAGATTTATGCGCCACTACCCAGTATGGATGAGATCTGAAACAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTAGCGGGCGCAGGCTTAGATGGAGGTACGGCGGCTGTCGCGGGATTCCGCGTGCCGGAAAAAATCGCCAATGGGCTTGTGGAGGGCTTCGGGAGAGGTTTTCAccgccccgccg                                                                                                                                                                                    >  1:141576/1‑250 (MQ=255)
                                                                    aacaacGGCAGCCTTCATAGTCCCCCCTTTCGGGATGATGTTCGGAAGAGAGGGTGAGGCAAAGGAGATTTATTCGCCACTACCCAGTATGGATGAGATCGGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTCCACCGCCACGCCt                                                                                                                                                                                    <  2:551413/251‑1 (MQ=255)
                                                                       aacTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGGGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCACGCCTtgc                                                                                                                                                                                 <  1:69676/251‑1 (MQ=255)
                                                                        acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCAAAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCGAAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGGGGTACGGGGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTGGCAGGGTTTCGGCAGAGTTTTTCACCGCCACGCCTtgct                                                                                                                                                                                <  1:208383/251‑1 (MQ=255)
                                                                                  tCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCCGCTTGCCAGAAATTATCGTCAAAGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTTTCCGCCCCGCCTTGCTGCTGCGccct                                                                                                                                                                      >  1:544064/1‑249 (MQ=255)
                                                                                     tAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTTCTGCTGCGCCAtttt                                                                                                                                                                   >  1:426995/1‑251 (MQ=255)
                                                                                                    ggATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGAAGTCGCGGGCATCGGCTTGGCAGAAAATATTCTCCAAAGGCCGTTGCAGGGCTTCGGGCGAGGTTTTTACCGCCCCCCCCTGCTGCTTCGCCACTTTGCCCACCGCaacacc                                                                                                                                                    >  2:128143/1‑247 (MQ=255)
                                                                                                            ttCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTTGGCTTGCCAGAAATTATCGTCAATGGCCCGGTGCAGGGGTTCGGCAGAGGTTTTCACCGCCCCGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACCCAATTgcgc                                                                                                                                          >  2:405883/1‑251 (MQ=255)
                                                                                                              cTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGGAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTTCTGCTGCGCCCTTTTGCCAACCCCCAACGCAATTGCGCgg                                                                                                                                        >  1:554991/1‑251 (MQ=255)
                                                                                                              cTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCCCCACGCCCTGCTGCTGCGCCATTTTTCCAACCGCAAACGCCATTGCGCgg                                                                                                                                        >  1:269748/1‑251 (MQ=255)
                                                                                                                gCATAGCAGGTGAGGAAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGcgcg                                                                                                                                         <  1:528238/248‑1 (MQ=255)
                                                                                                                          tGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCGTGCCAGAAATTTTTGGCAAATGGCTGTTTCAGGGGTTCGGCAGGGGGTTTTAACGCCACCCCTTGCTGCCGCGCCCTTTTTCCCACCGCAAAACGAAATTCGCGGGGGACTTTCttg                                                                                                                            >  2:175408/1‑249 (MQ=255)
                                                                                                                          tGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGTCGGGCTTTGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTTCCAACCGCCAACGCAATTGCGCGGGGGACTTTCTg                                                                                                                             >  1:528946/1‑250 (MQ=255)
                                                                                                                             ggaaaaTGAGATTTATTCGCCACTACCAAGTATGGATGAGATCTAAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAAAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAata                                                                                                                         <  1:175408/248‑1 (MQ=255)
                                                                                                                                  gtagttttgtttcccatcTACCCATTATGTATGAGATTTAAAAAAGGGGGGGGGAAATAGCCCGTTAGCCTTCCCTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGAACTGGCGTTATTCGGCTTGCCAAAAATTAACGTCAATGGCCTGTTGCAGGGCTCCGACAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCTCAAACGCAATTTCGCGGGAGACTTTCTGAATATCttt                                                                                                                    <  2:276668/234‑1 (MQ=255)
                                                                                                                                    gagaTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCCCGCCCTGCTGCTGCGCCCTTTTTCCAACCGCCAACGCAATTGCGCGGGGGACTTTCTGAATATCTTTc                                                                                                                   >  1:553823/1‑250 (MQ=255)
                                                                                                                                        tttATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCCATGGCCTGGTGCAGGGCTTCGGCAGAGGGTTTTACCGCCCCGCCTTGCTGCTGCGCCCTTTTGCCAACCGCAAACGCAATTGCGCGGGGGACTTTCTGAATATCttttcgttc                                                                                                              >  2:155030/1‑245 (MQ=255)
                                                                                                                                          tgttccccgctccccAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATGGCCCGGTAGCCTTCCCTACCGGGCGCAGGCTGAGAGGGAGGTAGGGGGGTGGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGTCTTCGGCAGAGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTcc                                                                                                             <  2:295582/238‑1 (MQ=255)
                                                                                                                                                  aCTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAATTATCGTCAATTGGCTGTTGCAGGGTTTCGGGCGAGGGTTTCACCGCCCCGCCTTTCTGCCGCGCCCTTTTTTCAACCCAAACAGCAAATTCGCGGGGGACTTTCTGGATATTTTTctgccatacggttctg                                                                                                    >  1:100230/1‑236 (MQ=255)
                                                                                                                                                      cccAGTAGGGATGAAATCTGAAAAAGGGGGGGGAAAATACCCTGGTACCCTCCACTACCGGGCCCAGGCTAGGAGGGGGTTACGGCGGTAGTCGCGGTTTCCGGCTTGCCAAAAATTATCGTCAATGCCCTTTTGCAGGCCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCAt                                                                                                  <  2:57775/249‑1 (MQ=255)
                                                                                                                                                          gTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAAATATTGTCCATGGCCTGTTGCCGGGCTTTGGCAGAGGGTTTCACCCCCACGCCCTGCTGGTGCGCCCTTTTTCCAACCGCCAACCCAAATGCGCGGGGGACTTTCTGGATATCTTTTAGTTCCGGCCGGACCATAccctc                                                                                            >  2:308736/1‑248 (MQ=255)
                                                                                                                                                            aTGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTTCAGGGCTTCGGCAGAGGGTTTTACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAAACGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCCGTACCATACCTTCGc                                                                                          >  1:334950/1‑251 (MQ=255)
                                                                                                                                                            aTGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCCATTGCGCGGGGGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGc                                                                                          >  1:261567/1‑251 (MQ=255)
                                                                                                                                                                aTGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGtt                                                                                      <  1:272084/251‑1 (MQ=255)
                                                                                                                                                                         gAAAAAGGGAGAGGGAAAAAGCCCGGTAGCCTTCCGTCCCGGGCGCAGGCTTAGATGGAGGTACGCCGGTAGCCGCGGTATTCGGCTTGCCAGAAATTATCGTAAATGGCCTGTTGCAGGGCTTCGGCATGGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCaa                                                                             <  1:249995/251‑1 (MQ=255)
                                                                                                                                                                              aGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATTGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTTCCAACCCCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCCCCGTTCCGCACCAATGGTg                                                                        >  1:132460/1‑251 (MQ=255)
                                                                                                                                                                              aGGGAGAGGGAAAAAGCCCGGTAGCCTTCACTACGGGGCGCAGGCTTAGATGGAGTTCGGGCGGTAGTCGCGGTATTCGGCTTGCAAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTg                                                                        <  1:155030/251‑1 (MQ=255)
                                                                                                                                                                                 gagaGGGAAATGCCCCGGTAGCCTTCCCTCCCGGGCGCAGTCTTAGATGGAGGTTCGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTTTTGCGGGGCTTCGGCAGAGGTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAAt                                                                     <  1:339704/251‑1 (MQ=255)
                                                                                                                                                                                     gggAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAAAGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCCTTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTg                                                                 >  1:360954/1‑251 (MQ=255)
                                                                                                                                                                                      ggAAATAGCCCGGTAGCCTTCACTACCGGGCGCGGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAACGGTGAATACTgc                                                                <  1:392164/251‑1 (MQ=255)
                                                                                                                                                                                                 ggTAGCCTTCACTCCCGGGCGCAGGCGTAGAGGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAAAAATTACCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGtttn                                                     <  2:261567/251‑2 (MQ=255)
                                                                                                                                                                                                    aGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAAACACGGTATTCGGCTTGCCAACAATTTACATAAAAGGCCTGTTTCAAGGCATCGGCCGAGGTTTTCACCCGCACGCCTTACTGCTTCGCCCTTTTTCCAACCCCGAAAACCATTTCGCCGGAGACGTTCTGAATATATTTCAGTTCCCGCAGTACCCTACCTTCCCCGCTCAGCACCAATTGTGTATACTGCGcccacgattccat                                                  >  1:249423/1‑240 (MQ=255)
                                                                                                                                                                                                       cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTTCTGCTGCGCCCCTTTTCCCACCCCCAACGCCAATGCGCGGGGGACTTTCTGAATATCTTTTCGTTCCCGCCGGACCACACCTTCCCCGTTCAGCACCCATGGTGAATAATGCGCCCGCGGTTCCCTTGc                                               >  2:48770/1‑251 (MQ=255)
                                                                                                                                                                                                               ccGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGGTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTTCGCCATTTTTCCAACCGCAAACGCAATTGCGCGGGAGGCTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCCGCACCAATGGTGAATACTGCGCCAGCGGTTCCATTGCCGACATc                                        >  1:46543/1‑250 (MQ=255)
                                                                                                                                                                                                                        ggCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAAAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGAAGGGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCAt                              <  1:78356/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCCCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCCTACCTTCGCCGTTTCGCCCCAATGGTGAAAACTGCGCCAGCGTTTCACTTGCCGACATCAGCCTCTCAt                              >  2:54211/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gAGGTACGGCGTTAGTCGCGGTATTCGGCTGGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACg                    <  1:405883/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gggcggTAGTGGGGGTATTCGGCTTGCCAGAATTTATCGTCATGGGCCTGTTGCAGGGCTTCGGCAGAGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGAcgc              <  2:554991/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ggagggaGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGTGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATTCCTTCGCCGTTCAGCACCAATGGTGCATACTGCGCCAGCGTTTCACTTGCCGACAGCAGCATCTCATCTGTGATACGTGAcg               >  2:152632/7‑249 (MQ=255)
                                                                                                                                                                                                                                            ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCg          >  2:52613/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gTCGCGGTTTTCGGCTACCCAGAAATTATCGTAAATGGCCTTTTGCAGGGTTTCGGCGGGGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGAAg      <  2:269748/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccgcgGTATTCGGCTTGCCAGAAATTGTCGTCAATGCCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGAAGc     <  1:21980/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACCTCAGCATCTCATCGGTGATACGTGACGCGCCGGAAGCAAt  >  2:383028/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
AGCAGGGCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGAAGCAAT  >  NC_000913/1553753‑1554247

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