Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,8550AC32.9% 119.8 / 18.8 70G247G (GGT→GGG) gpmA2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (36/11);  new base C (0/23);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.60e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786608‑787105
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           
gATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                                                                                                                                                                                                                                          >  1:186695/1‑251 (MQ=255)
       gCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGtt                                                                                                                                                                                                                                                    <  1:238373/249‑1 (MQ=255)
         aTCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGAGCGTCGCGGGGAAGGTTTAATTGTCGGATGCGCCGGCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTT‑TCATTTTAAACGAATGACGTTTACTTTGCTTTACCCT‑GGTTTGc                                                                                                                                                                                                                                                 >  1:452604/1‑250 (MQ=255)
           cGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCCC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAc                                                                                                                                                                                                                                              >  1:107230/1‑251 (MQ=255)
            ggatGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAcc                                                                                                                                                                                                                                             <  2:242341/248‑1 (MQ=255)
               tGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAccgcc                                                                                                                                                                                                                                          >  1:357884/1‑251 (MQ=255)
                   aaGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCgctg                                                                                                                                                                                                                                      >  1:158622/1‑251 (MQ=255)
                    aGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCgctgc                                                                                                                                                                                                                                     >  1:144454/1‑251 (MQ=255)
                      aTACCGATGCTCTTTTACGTCAACTGGCGAAAGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑CCAGCGTCGCCTCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACCTTTACTTCGCTTTACCCT‑GGTTTTCCACCGCCGCTGCtt                                                                                                                                                                                                                                   >  2:124383/1‑251 (MQ=255)
                             tGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                                                                                                                            >  2:54345/1‑251 (MQ=255)
                             ggcgttgttgcgaaaaCTGGCGAAGGCAGTAAAAGCGGGTCGGGGGTAAGGTTTAATGGGCGGAGGCGCCGTAAGAGGGGCG‑TATCCGAGAAATCCCCACAGGCCTGATAAGTCGC‑GCGGCGTCGCATCAGGCAATGTGCTCCATGGTTAGCAAAAAAAAAGCCGAC‑TCCCTTGCGGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTCCTTCGCTTACCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                                                                                                                            <  2:104714/237‑1 (MQ=255)
                                   gTTACGTCAACTGGCGAATGCAGAAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCT‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGCCGCATCAGGCAATGTGCTCCATTGTTAGAAACAAAAAAGCCGAC‑TCACTTGCAGTCGGTTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                                                                                                                      <  1:73031/251‑1 (MQ=255)
                                      aCGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                                                                                                                                   >  2:339781/1‑251 (MQ=255)
                                        gacAACGGGCGAATGCAGAAAAAGCGGGGGGGGGGTAAGGTTAAATGGGCGGATGCGCCGTCAGAGGGGCG‑TATCCGAGGAATCACCACAGGCCTGATAAGGCGC‑GCAGGGTCGCATCAGGAAATGTGCTCCATGGTTAGAAAAAAAAAAGCCGAC‑CCACTTGCAGTCGGTTTTC‑TCATTTTAAACGAATGCCTTTTACTTCTCTTTCCCCT‑GGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                                                                                                                                 <  2:175138/249‑1 (MQ=255)
                                          cAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGTCAATGTGCTCCCTTGTTAGCAACCAATAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCCTTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCCCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                                                                                                                                >  2:212239/1‑250 (MQ=255)
                                             cTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTTGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGCTCCCGTCAGCAtt                                                                                                                                                                                                             >  2:388561/1‑250 (MQ=255)
                                              gggggAATGCAGTAAAGGGGCGTCGCGGGTAGGGTTTAATTGTGGGATGCGCCGCGAGAGGGGCG‑TATCCGATAAACCACCACAGGCCTGATAAGGCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTTTAAGCAAAAAAAAAGCCGAC‑TCTCTTGCAGTGGGCTTTC‑TCATTTTAAACGAATGACGTTTCCTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                                                                                                                           <  2:24887/247‑1 (MQ=255)
                                                  gAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGAGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAACAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  1:374216/1‑249 (MQ=255)
                                                  gAATGCAGCAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  1:99649/1‑249 (MQ=255)
                                                   aaTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  1:320788/1‑248 (MQ=255)
                                                   aaTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  2:74195/1‑248 (MQ=255)
                                                    aTGCAGTAAAAGCGCGTCGGGGGTAAGGTTTAATTGGGGGAGGCGCCGTCAGAGGGGCT‑TCCCCGAGGAACCACCACAGGCCTGATAAGTCGC‑GCAGCGCCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                                                                                                                                    <  1:188672/236‑1 (MQ=255)
                                                    aTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                                                                                                                                    >  2:188672/1‑236 (MQ=255)
                                                      gCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                                                                                                                    >  1:158057/1‑250 (MQ=255)
                                                      gCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                                                                                                                    >  2:72714/1‑250 (MQ=255)
                                                          aaaaaGCGCGGGGGGGGTAAGGTTTAAGTGTGGGAGGGGCCGTCAGAGGGGGG‑TAACCGAAGAACCCCCACAGCCCGGATAAGGCGC‑GCAGCGTCGCATCAGGCAGTGGGCTCCATTTTTAGCAACAAAAAAGCCGCC‑TCACTGGCAGGCGGTTTTC‑TCATTTAAAACGATGGACGTTTACTTCGTTTTACCCG‑TTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                                                                                                                               <  2:387798/250‑1 (MQ=255)
                                                           aaaaGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                                                                                                                              >  1:189407/1‑251 (MQ=255)
                                                               gcgcGTCGGGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGGGGCG‑TATCCGATGAATCCCCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGTTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                                                                                                                          <  1:458685/251‑1 (MQ=255)
                                                                   ggggcgGTTAAGTTTTAATTGTCGGTGGCGCCGACAGAGGGGCG‑TACCCGATAAATCACCACAGCCCTAAGAAGTCGC‑GCAGCCTCGCAACAGGCAATGTGCCCCTTTGTTAGCAACAAAAAAGCCGAC‑TCCCTTGCAGGCGGCTTTC‑TCATTTTAACCGAATGACGTTTACTCCCCTTTCCCCT‑GGTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                                                                                                                                      <  2:399147/248‑1 (MQ=255)
                                                                        gggTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                                                                                                                                 >  1:377872/1‑251 (MQ=255)
                                                                                              gcgcCGTCAGAG‑GGCGGTATCCGATGAATCCCCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCCCTGGCGGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                                                                                                                                           <  1:39240/251‑1 (MQ=255)
                                                                                                    ccAGAGTGGCG‑TACCCGATAAATCACCACAGGCCGGATAAGTCGCCGC‑GCGTCGAAACGGGCAATGTGCCCCATTGTAAGCAAAAAAAAATCCGAC‑TCTCTTGCTGTCGGTTTTC‑TCATTTAAAACGAATGACTTTTGCTTCGCTTTACCCG‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                                                                                                                                     <  1:240797/250‑1 (MQ=255)
                                                                                                      agagTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGc                                                                                                                                                   >  2:241702/1‑251 (MQ=255)
                                                                                                         gTGGCG‑TACCCGTTGACTCACCACGGGCCGTATAAGCCGC‑CCACCGTCACTTCAGGCTGTGTGCTCAATGATGAGAAAAAAAAACGCCGAC‑TCACTTCCATGCGGTTTTC‑CCTTTTTAAACTAAGGTCTTTTACCTCGCTTTCCCCG‑TGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                                                                                                                                <  1:362206/251‑1 (MQ=255)
                                                                                                                          atcacaccgggCCTGATAAGCCGC‑GCAGCGTCGCATAAGGCAATGTGCCCCATTGTAAACAACAAAAAAGCCGAC‑TCCCTGGCGGTGGGTTTTC‑TCATTTTAAACGAATGCCTTTTACTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTaa                                                                                                                                <  2:44020/242‑1 (MQ=255)
                                                                                                                                     ccTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                                                                                                                     >  1:182127/1‑251 (MQ=255)
                                                                                                                                          aaaagccgcgcggcgccccATCGGGCAATGTGCCCCATTTTAAGCAAAAAAAAAGCCGCC‑TCTCTTGCAGTGGGTTTTC‑TCATTTTAAACGAAGGACGTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt                                                                                                                 <  2:186695/233‑1 (MQ=255)
                                                                                                                                           aaGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct                                                                                                               >  1:143692/1‑251 (MQ=255)
                                                                                                                                            aGTCGC‑GCAGCCTCGCATCAGCCAATGTGCCCCATTGTTAGCAAAAAAAAACCCGCC‑CCCCTTGCAGTCGGCTTTC‑TCATTTAAAACAAATGATTTTTACTTCGCTTTCCCTT‑GGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt                                                                                                              <  1:339781/251‑1 (MQ=255)
                                                                                                                                             gTCGC‑GCGGCGCCCCATCAGGCAATGTGCTCCATTTTGAGAAAAAAAAAACCCGCC‑TCCTTTGCATTCGGTTTTC‑TCATTTTAAACGAATGACTTTTCCTTCCCTTTCCCCG‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                                                                                             <  1:74095/251‑1 (MQ=255)
                                                                                                                                             gTCGC‑GCAGGGCCGCATCAGGCAATTTGCTCTATTGTTAACAAAAAAAAAGCCGCC‑TCACTTGCAGTCGGCTTCC‑CAATTTAAAACGAATGACGTTTCCTTCGCTTTCCCCG‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                                                                                             <  1:75274/251‑1 (MQ=255)
                                                                                                                                                gc‑gcAGCGTCGCATCAGGAAATGTGCTCCATTGTTAGAAAAAAAAAAGCCGCC‑TCACTTGCAGTCGGCTTTC‑CCATTTAAAACGAATGAGTTTTACTTCGCTTTCCCTT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt                                                                                                          <  1:206644/251‑1 (MQ=255)
                                                                                                                                                             aTCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCa                                                                                                >  2:192565/1‑249 (MQ=255)
                                                                                                                                                               cagcaaagtggCTCCTTTGTTAACAACAAAAAACCCGCC‑TCACTTGCTTTCTTTTTCC‑TCATTTAAACCAAATGACTTTTCCTTCGCTTTCCCCG‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGat                                                                                            <  1:421376/241‑1 (MQ=255)
                                                                                                                                                                   aaaTGTGCCCCTTTGTAAGCAACAAAAAAGCCGAC‑TCCTTTGCAGTCG‑CTTTCCTCATTTAAAACGAATGCCGTTTCCTTCGTTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATAttt                                                                                        <  1:345527/250‑1 (MQ=255)
                                                                                                                                                                       ttccTCCATTGTAAGAAAAAAAAAAGCCCAC‑TCATTTGCGGTCGGTTTTC‑TCATTTTAAACGAATGACTTTTACTTCGCTTTACCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCAc                                                                                     <  1:54345/247‑1 (MQ=255)
                                                                                                                                                                        tgCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCAcc                                                                                    >  2:355276/1‑250 (MQ=255)
                                                                                                                                                                            ccATTGTTAGCAACAAAAAAACCGCC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTCCCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCCCGCCAGTCGGGATATTAAGCTCAAGACTCTCTTCTTCGCTCATTTTATCAAGATATTTCACCAgcgc                                                                               >  2:381769/1‑251 (MQ=255)
                                                                                                                                                                               ttgtaaacaacaaAAAAGCCGAC‑TCCTTTGCGGTCGGTTTTC‑TCATTTTAAAGAAATGACTTTTACTCCGTTTTCCCCT‑GTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCAc                                                                             <  1:432724/243‑1 (MQ=255)
                                                                                                                                                                                   gaGCAACAAAAA‑GCCCACCTCTCTTGCAGTCGGTTTTC‑TCTTTTTAAACAATTGACTTTTCCTTCGCTGTCCCTT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaaa                                                                        <  2:384868/250‑1 (MQ=255)
                                                                                                                                                                                      aaaaaaaaaaGCCTCC‑TCATTTGCAGTGGGTTTTC‑TAATTTAAAACAAAGGACTTTTACTCCG‑TTTACCCTTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAagag                                                                     <  2:247782/247‑1 (MQ=255)
                                                                                                                                                                                            aaaaGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCg                                                               >  2:204231/1‑251 (MQ=255)
                                                                                                                                                                                              aaGCCGAC‑TCACTTGCAGTCGTTTTTC‑TCATTTAAAACGAATGACGTTTACTCCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGATTTACCgtg                                                             <  1:72714/251‑1 (MQ=255)
                                                                                                                                                                                              aaGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCgtg                                                             >  1:403479/1‑251 (MQ=255)
                                                                                                                                                                                               aacccgtctcattttcagtcttttttC‑TCTTTTAAAACGAATGACTTTTTCTTCTCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCgtg                                                             <  1:245621/228‑1 (MQ=255)
                                                                                                                                                                                               aGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAACGCCGCTGCTTTCGCTGCGATCTCGTCAGCCTTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAACAGCGGCACGCCAGTCGGGATATTAAGCTCAAGACTCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTACAGAGTTACCgtgg                                                            >  2:380656/1‑250 (MQ=255)
                                                                                                                                                                                                                ggcGGCTTTC‑TCATTTTAAACGAATGCCTTTTCCTTCGCTTTCCCCG‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcgc                                            <  1:81848/249‑1 (MQ=255)
                                                                                                                                                                                                                   ggCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCt                                         >  1:415122/1‑251 (MQ=255)
                                                                                                                                                                                                                               tttAAACGAATGACTTTTCCTTCGCTTTCCCCG‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgat                                                    <  1:94117/229‑1 (MQ=255)
                                                                                                                                                                                                                               tttAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgat                                                    >  2:94117/1‑229 (MQ=255)
                                                                                                                                                                                                                                 tAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCAt                            >  1:370119/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaaacgattgttgTTTACTTCGTTTTCCCTT‑GTTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCa                             <  2:418414/238‑1 (MQ=255)
                                                                                                                                                                                                                                  aaccaattGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCAt                            <  1:207388/243‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGACGTTTCCTTCGTTTTCCCCT‑GTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                     <  1:236608/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                     >  2:278325/1‑249 (MQ=255)
                                                                                                                                                                                                                                           cggTTTCCTTCGCTTTCCCCT‑GGTTCGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGa                  <  1:38744/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    tcgttttACCCT‑GGTTTGCCCCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTc          <  2:18359/246‑1 (MQ=255)
                                                                                                                                                                                                                                                    tCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTc          >  2:191741/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     cGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCAt        >  2:19987/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     cGCTTTAACCT‑GGTTTGCAACCGCCGCTGGTTTCGCTGCGATCTCGTCAGCATTACCCAGATACTAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCCGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCCGAATCGTTTCAt        >  2:153158/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           aCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGCATAGTTTCATTCCAGt  >  2:318329/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           
GATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCG‑TATCCGATGAATCACCACAGGCCTGATAAGTCGC‑GCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGAC‑TCACTTGCAGTCGGCTTTC‑TCATTTTAAACGAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786608‑787105

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