Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,378,9680CT34.2% 36.1 / 16.8 41E98K (GAA→AAA) menD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base C (8/19);  new base T (12/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 9.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGCGTTAATGCTT  >  NC_000913/2378733‑2379213
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ccATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGc                                                                                                                                                                                                                                        <  1:367975/251‑1 (MQ=255)
 catccAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCt                                                                                                                                                                                                                                       <  1:373366/247‑1 (MQ=255)
 cATACAGCGGTTCAGCAAACGGGCAGTTGATATGGAATCACCCCGCATGAAGCGTACCGCGAGCGTGGTCGATGGTAGAAACCAGCCAACGGGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGTCGGGTTAGAGAATAAGTATTTTTcccggtccaccgcgt                                                                                                                                                                                                                                        >  2:347834/1‑238 (MQ=255)
    aCAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCTCATGAAGCGTACCGAGCGCGTGGTCGATGGTAGAAACCAGCCATCGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CACGGCtt                                                                                                                                                                                                                                      <  2:8207/249‑1 (MQ=255)
    aCAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCa                                                                                                                                                                                                                                    <  1:442997/251‑1 (MQ=255)
          gTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCTTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCCTCAATCAGt                                                                                                                                                                                                                              <  2:131/251‑1 (MQ=255)
             cAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGc                                                                                                                                                                                                                            <  2:310172/250‑1 (MQ=255)
               gCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCAGTGTTTCTCCCGGTAA‑CCCGGCTTCAAACAGTGCCgg                                                                                                                                                                                                                         >  1:213924/1‑251 (MQ=255)
                     gggCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGGCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATagag                                                                                                                                                                                                                   <  2:19664/251‑1 (MQ=255)
                       gCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGGAACCGATTTTTTCCGGGTAA‑CCCCGCTTCaaaccctgcccgggaagaat                                                                                                                                                                                                                 >  2:519917/1‑233 (MQ=255)
                         aGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCATGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGAGAATTCGCTTTTTTCCGGGTAA‑CACGGCGTCAAACAGAGCCGGATAgagagtt                                                                                                                                                                                                               >  1:25882/1‑248 (MQ=255)
                                     cTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGGGAAATAGCGTTTTTCCGGTTAAACCGGGCTTCA‑TTAGAGCCGGAAAGAGGATTGGCAccgccttg                                                                                                                                                                                                   >  2:427446/1‑248 (MQ=255)
                                            cGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAAACAGTTTTTCTCCCGTTAA‑CCCGGCCTCCATTAGTGCCGGATAGAGAATTTACAccgccggggcggggg                                                                                                                                                                                            >  2:408214/1‑242 (MQ=255)
                                               aTGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAAAGGAACAGTTTTTTTCCGGTTTA‑CCCGGCGTCAATCAGTGCAGGGTAATGAATTTCCACCGCCGTGCCGGGGGTcc                                                                                                                                                                                         >  1:298898/1‑250 (MQ=255)
                                               aTGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGGGAATCAGGGTTTTTTCGGTTTA‑CCCAGCCTCAATCCGTGCCGGGAAGGGATTTTTCAccgccggtgcggggggcc                                                                                                                                                                                         >  2:287134/1‑239 (MQ=255)
                                                          cGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGGTTTTTTCCGGGTAA‑ACCGGCTTCAATCAGTGCCGGGTAGAGATTTGCCACCGCCCGGCCGGGGGTCCCAATCACCGcc                                                                                                                                                                              >  1:87486/1‑251 (MQ=255)
                                                                     gggagggTAGAAACCAGCCAACGTGCGGGGATAGCCTGGGT‑CGGGCGCGGCAATGAAATACTGGGCGTGGGGTGAGAGGCGAACATCCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTAAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCGCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCtgctt                                                                                                                                                                   <  2:458223/246‑1 (MQ=255)
                                                                       gATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGCTTAGAGGATCAGGTTTTTTCCGGTTAA‑CCCGGCTTCAATCAATGCCGGATTGAGATTTGCCACCGCCGGGCCGGGGGTCACAATCAcccccccccgccgcttgg                                                                                                                                                                 >  2:411082/1‑235 (MQ=255)
                                                                              gAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGGATCAGTTTTTTTCCCGTTTA‑CCCCGCTTCAATCAGGGCCCGATAGGGATTTGCCACCCCCCTGGCGGAGGGCACAATCAcccccccccgccgccttctgccttt                                                                                                                                                          >  1:444730/1‑228 (MQ=255)
                                                                                agcaagccacagTGCGGGGATATCCTGGGC‑CGGGCGCGGCAAGGAAAAACTGTGCGGGGGGGGAGAGGCGAACATGCCCGGGGGGAGAATTGCCGGATGCGCGCTGCAGCCAATGAGCTCGGGGGGGCGATCGGCGGTTAAGATAATCAGTTTTGCTCCGGTGAA‑CCCGGCTTCAATCATGGCGGGATAGAAATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTGGCTGACTTTCg                                                                                                                                                        <  2:50780/240‑1 (MQ=255)
                                                                                     gCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGGGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGGTTTTTTTCGGTTTA‑CCCCGCTTCTAACAGGGCCGGGGAGAGATTTTCCACCGCCGTGCCGGGGGGCACAATCACCGCCAACCGCCGCTTGCTGACTTTTGcccgc                                                                                                                                                   >  2:78724/1‑248 (MQ=255)
                                                                                       cAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCCGTTTTTCTCCGGTTTA‑ACCGGCTTCAATCAGTGCCGGATAGAGATTTTCCACCGCCGTGCCGGGGGGCACAATCACCGCCACCGGCCGCCTGCTGACTTTTCCCAGccc                                                                                                                                                 >  1:68325/1‑251 (MQ=255)
                                                                                       cAACGTGCGGGGATATCC‑GGGTGCGGGCGCGGCAAAGAAATACTGTGCGTGGGTTGATAGGCAAACCTGCCCGGCTGGCAAATTGCCTGGTTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAAAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGGGCCGGATTGAGATTTTCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTCCGCCAGccc                                                                                                                                                 <  2:223233/251‑1 (MQ=255)
                                                                                        aaCGTGCGGGGATATCCTGGGC‑CGGGCGGGGCAAAGAAATACTGGGCGTGGGTTGAGAGCCGAACTTCCCCGGCTGGCGAATTGCCGAATTCGCCCCGCAGGCAATTAGCTCCGGCGGGCGACCGGGGGTGAAGAGAATCAGTTTTTCGCCGGTTAA‑CCCGGTTTCAATCGGGGCCGGATAGAGATTTGCCGCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGcccc                                                                                                                                                <  2:165582/251‑1 (MQ=255)
                                                                                            ggCGGGGATATCCGGGGT‑GGGGCGAGGCAAGGAAATACTGTGAGTGGGGAGAGAGGCGAAAATGCCCGGCTGGCGAATTGCGGGATTCGCGCCGCAGTCAATTAGCTCCGGGGGGCGATCGGGGGTTAAGAGAATAAGTTTTTTTCCGGTTAC‑CCCGGCTACAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAACCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAgcg                                                                                                                                            <  1:196503/250‑1 (MQ=255)
                                                                                                                      ggCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATc                                                                                                                   <  1:345153/251‑1 (MQ=255)
                                                                                                                         aaTGAAATACTGTGCGTGGGGGGAGAGGCGAACATCCCCGGCGGGCGAATTGCCTGATGCGCCCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTCCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATGGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCTCCAGCCCCAGCGCCAGATTCCCCAACCCACGCTCATCGaa                                                                                                                <  2:257149/251‑1 (MQ=255)
                                                                                                                               aTACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAGAGAAACCGATTTTTTCCGGGTGA‑CCCGGCTTCAATCAATGCCGGGTAGAGGTTTGCCACCGCCGTGCCGGGGGTCACAATCACCGCCCCCCGCTGCTTCCTGGCTTTTTCCACCCCCACCGCCCGATGCCCCCACCCCCCCTCATCGAAATgggg                                                                                                          >  1:281360/1‑250 (MQ=255)
                                                                                                                                                              ccggtgggggAATTGCCTGATTCGCGCCGCAGTCAATTAGCCCCGGGGGGCGATCGGCGGTTAAGAGAAGCAGTTTTTCTCCGGTTAC‑CCCGGTTTCAATCAGTGCCGGATAGAGATTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTccgccgcc                                                                           <  2:403910/242‑1 (MQ=255)
                                                                                                                                                               cGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGGATCAGTTTTTTTCCGGTTAA‑ACCGGCTTCAATCAGTGCCGGAGAAAGATTTTGCACCCCCGGGCCGGGGGGCACAATCACCCCCACCCGCTGCTTGGTGGCTTTCCCCCGCCCCCGCGCCCGATGCCCCCACCCCCCCTCATCGAAATGGGGGTGGGGGATGAATGGGGGAGTCTcccccgccc                                                                          >  2:384351/1‑250 (MQ=255)
                                                                                                                                                                ggCTGGTGAATTGCATGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGATAACAGAAACCGGCTTTCTCCGGTTCA‑ACCGGCCTCAAAAAAGGGCGGATAGAGAATTGCCACCCCCACCCCGGGGGGCACAAACAACCCCAACGGAAGCTTGCTGGCTTTCGCCAGACCCACCCGCAGACGCCCCAACACCCCCTCATCGCACTAGGTGTGGCGAATGAATGCGGGATTCTccacctccc                                                                          >  2:347273/1‑243 (MQ=255)
                                                                                                                                                                                       gccgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGGATAAAATTTTCTCCGGGTAA‑ACAGGCGTCAATCAGAGCCGGATAGAGATTTTTCAACGCCGTGCCGGAGGGCCAAATCACCGCCACCGGCTGCTCGGTGGCTTTCGCCAACACCACCGCCAAATGCCCCAACCCACGCTCAATGAAATTGGGGGGGTGTATGGATGAGGGATTCTCAGACGCCGCCAACGTTTACGGTGTGGaaaagg                                                  >  1:289327/1‑247 (MQ=255)
                                                                                                                                                                                        ccgcAGTCAATTAGCTCCGGCGGGCGATTGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGa                                                 <  1:138932/251‑1 (MQ=255)
                                                                                                                                                                                         cgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGGATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACTTTCACCGCCACCCGCTGCTTGCTGACTTTCTCCCGCCCCAGCGCCAGATGCCCCAACCCCCCCCCATCGAAATggg                                                                                                           >  1:410454/1‑192 (MQ=255)
                                                                                                                                                                                                 aTTAGCCCCGGCGGCCGATCGGCTGTTAAGAGAACCAGTTTTTCCCCGGTTAC‑CCCGGCTGCATTCAGTGCCGGATAGAGATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCATCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGTGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGc                                        <  1:259337/251‑1 (MQ=255)
                                                                                                                                                                                                           gcgggcgACCGGCGGTTAAGAGAAGTAGTTTTTCTCCGGTTAC‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCTCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGTCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGAtg                              <  1:361795/251‑1 (MQ=255)
                                                                                                                                                                                                              ggcgatgggggtttAAGAGAATCATTTTTTCCCGGGTTAC‑CCCGCCTTCAATCAGTGCCGGAAAGAGATTGGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCCCCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATgtgt                           <  2:7389/239‑1 (MQ=255)
                                                                                                                                                                                                                            aaGAGAATCAGTATTTTTCCGTTTAA‑CCGGGCCTCAATCATGGCCGTATAGATATTTCCCACCGCCGTGCCGGAGTTCACAATCACCGCCACCGGCGGCTGGCTGCCTTGCGCCAGCCCCAGCGCCAGATGCCCCACCCCGCGCTCAGCGAAATGTGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCTGTGAcg             <  1:87733/251‑1 (MQ=255)
                                                                                                                                                                                                                                gaATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGTATTCTCCGCCGCCGCTaa                                                                      >  1:480551/1‑190 (MQ=255)
                                                                                                                                                                                                                                 aaTCGGTTTTTCCCCGGTAAC‑CCCGGCTTCAAACAGTGCGGGATAGAGTTTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCGGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGGGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGCGTTa        <  2:68325/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tttttCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGCCGCCGTGACGCGTTACTGCtt  >  2:164142/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
CCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGT‑CGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTCCGGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGCGTTAATGCTT  >  NC_000913/2378733‑2379213

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