Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,8550AC24.1% 117.9 / 12.5 54G247G (GGT→GGGgpmA2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (35/6);  new base C (0/13);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 2.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCT  >  NC_000913/786718‑786991
                                                                                                                                         |                                                                                                                                        
gTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAcc                                                                                                                                        <  2:1302377/140‑1 (MQ=255)
gTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAcc                                                                                                                                        >  2:1415611/1‑140 (MQ=255)
 tATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                                                                                                                       >  1:675550/1‑140 (MQ=255)
 tATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                                                                                                                       <  2:742337/140‑1 (MQ=255)
   tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                     >  2:503969/1‑140 (MQ=255)
    ccGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTgg                                                                                                                                    >  1:348041/1‑140 (MQ=255)
        tGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                >  1:572447/1‑140 (MQ=255)
        tGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                >  1:269533/1‑140 (MQ=255)
        tGAATAACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                                                                                                                >  1:306998/1‑140 (MQ=255)
                cacaGGCCTGAAAAGTCGCGCAGCGCCGCATCAGGCAATGTGCTCCATTGTTAGAAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAAccgcc                                                                                                                        <  2:247574/140‑1 (MQ=255)
                 acaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgccg                                                                                                                       >  2:748381/1‑140 (MQ=255)
                  caGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACcgccgc                                                                                                                      >  2:877330/1‑140 (MQ=255)
                      ccTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgctgct                                                                                                                  >  2:510434/1‑140 (MQ=255)
                      ccTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgctgct                                                                                                                  >  2:551806/1‑140 (MQ=255)
                       cTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCtt                                                                                                                 <  1:938880/140‑1 (MQ=255)
                        tGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCTTTTTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCttt                                                                                                                <  2:1401146/140‑1 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                              >  1:27278/1‑140 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                              >  1:291155/1‑140 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                              >  1:34676/1‑140 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                              >  1:480122/1‑140 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                              >  2:1149768/1‑140 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                              >  2:1008285/1‑140 (MQ=255)
                          aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGc                                                                                                             >  2:232130/1‑140 (MQ=255)
                             aGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                           >  2:1130928/1‑140 (MQ=255)
                               tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                         >  1:140293/1‑140 (MQ=255)
                                      gCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                  <  1:242455/140‑1 (MQ=255)
                                      gCGTCGCATCAGCCATTGTCCTCCTTTGTCAAAAACAAAAAAGCCTCCTCTCTTGCAGTGGGTTTCCTCATTTTAAGCTAATGTGTTTTCTTTCGCTTTCCCCTGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                  <  2:291155/140‑1 (MQ=255)
                                      gCGTCGCACCAGGCAATGTGCTCCATTGTTAGAAAAAAAAAACCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                  <  2:140293/140‑1 (MQ=255)
                                      gCGTCCCATCAGGAAATGTGCTCCTTTGTTACAAAAAAAAAAGCCGACTCCTTTGCAGTCGGCTTTCTCTTTTAAAACGAATGACTTTTCCTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                  <  2:367754/140‑1 (MQ=255)
                                       cGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                 >  2:530237/1‑140 (MQ=255)
                                           gcctCAGGAAATGTGCTCCATTTTTAGAAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAACCGAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCa                                                                                             <  2:572447/137‑1 (MQ=255)
                                           gCATCAGGCAATGTGCTCCATTTTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCa                                                                                             <  2:846347/140‑1 (MQ=255)
                                             aTCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAtt                                                                                           <  2:27278/140‑1 (MQ=255)
                                               cAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                         <  2:399313/140‑1 (MQ=255)
                                                 ggCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                       <  2:1433347/140‑1 (MQ=255)
                                                    agtgtgCTCTATTTTTAGAAAAAAAAAAGCCGACTCGCTTGCAGTCGGTTTGCTATTTTTAAACGAATGACGTTTACTTCGTTTTACCTTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGa                                                                                    <  2:1461143/138‑1 (MQ=255)
                                                      tgtgCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                  >  2:275649/1‑140 (MQ=255)
                                                           tCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                             >  1:521960/1‑140 (MQ=255)
                                                             cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                           >  2:682301/1‑140 (MQ=255)
                                                                  ttAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAg                                                                      >  1:1222197/1‑140 (MQ=255)
                                                                    aGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                    >  2:38126/1‑140 (MQ=255)
                                                                    aGCAACAAAAAACCCGCCTCCCTTGCGGTCGGTTTTCTCATTTTAAACGAATGACTTTTTTTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                    <  2:1222197/140‑1 (MQ=255)
                                                                               aGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtc                                                         <  2:34676/140‑1 (MQ=255)
                                                                                   gccTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                     <  2:262582/138‑1 (MQ=255)
                                                                                    aCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                    >  1:1172928/1‑140 (MQ=255)
                                                                                    aCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                    >  1:676379/1‑140 (MQ=255)
                                                                                     cTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                   <  1:1008285/140‑1 (MQ=255)
                                                                                      tCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGa                                                  <  1:551806/140‑1 (MQ=255)
                                                                                               gTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                         >  1:130847/1‑140 (MQ=255)
                                                                                                tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                        >  1:465110/1‑140 (MQ=255)
                                                                                                   gCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATa                                           <  1:143272/134‑1 (MQ=255)
                                                                                                   gCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATa                                           >  2:143272/1‑134 (MQ=255)
                                                                                                    cTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                    >  2:1224167/1‑140 (MQ=255)
                                                                                                    cTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                    >  1:1220417/1‑140 (MQ=255)
                                                                                                           tttttAAACGAAGGTGTTTTCCCTCGCTTTCCCGTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGc                             <  2:465110/139‑1 (MQ=255)
                                                                                                            tttaAAACGAATGACGTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                            <  2:306998/140‑1 (MQ=255)
                                                                                                              ttAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                          <  1:275649/140‑1 (MQ=255)
                                                                                                               tAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                         <  1:1329859/140‑1 (MQ=255)
                                                                                                                aaaCGAATGACTTTTCCTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                        <  2:675550/140‑1 (MQ=255)
                                                                                                                 agcaaaTGTCCTTTCCTTCTCTTTCCCCTGGTTTAAAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                       <  2:521960/136‑1 (MQ=255)
                                                                                                                   cGAATGACTTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                     <  2:167584/140‑1 (MQ=255)
                                                                                                                     aaTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                   >  2:51820/1‑140 (MQ=255)
                                                                                                                        gccttttCCTTCCTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                <  2:676379/136‑1 (MQ=255)
                                                                                                                        gACTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                <  2:269533/140‑1 (MQ=255)
                                                                                                                        gACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                <  1:232130/140‑1 (MQ=255)
                                                                                                                           ttttCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg             <  1:1415611/139‑1 (MQ=255)
                                                                                                                                    gCTTTACCTTGTTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt    <  2:57111/140‑1 (MQ=255)
                                                                                                                                    gCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt    >  2:233446/1‑140 (MQ=255)
                                                                                                                                      tttCCCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct  <  1:1437261/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCT  >  NC_000913/786718‑786991

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