Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,845 A→C 19.2% V197G (GTG→GGG)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,8450AC19.2% 120.5 / 3.6 62V197G (GTG→GGG) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (28/17);  new base C (3/12);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 6.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCG  >  NC_000913/2967611‑2968080
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
cctgTTTTTCCCCCCCACCCTGCCCCCG‑GGCTTTAAACACCTCCGAAAAAAATTCCCCGACGGAGCGGGGTTTCCGAAGGCGGGTGGGAGAAAAGACCAAAAAGAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATacac                                                                                                                                                                                                                                              <  2:98722/233‑1 (MQ=255)
aCCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTGGCGTTCCAGG‑GCCGGATCCAGCCT‑TGATGCCTGTTAcccccg                                                                                                                                                                                                                                           >  1:1293154/1‑237 (MQ=255)
            gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTTGTTTGCCGCTTCTTCCACGGCCCCGGTCAA‑TCGTTCGCGGTCCTAG‑GCCGGCTCCACGCT‑TGTGCCCTGATAAAACC‑GTCCCATTaaaag                                                                                                                                                                                                                              >  2:1281508/1‑237 (MQ=255)
            gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCCTTTTCCGCTTCTTCCACCGCCCCGGGCAG‑TCGTTCGCGTTCCAGC‑GCCGGATCCCACGT‑GGATGCCCGATACAACT‑GGTCCATTAAAgg                                                                                                                                                                                                                              >  2:1062448/1‑239 (MQ=255)
                            tctttttATCAACCCGCGAAAAAAATT‑CCGGGCGGCGCCGGGTGTCCGAAAGGAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGTTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGc                                                                                                                                                                                                              <  2:1293154/235‑1 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTTTTCCGCTTCTTCCAGCGCCCCGGGCGG‑TCTCTCCCGTTCCCGA‑CCCGGCGCCA‑CGTGTGAAGCCTGGCACCTCT‑GCTCCCTTAAAAGGAAAGGGGcctcctt                                                                                                                                                                                                               >  2:549934/1‑237 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCCTTAAAGGGAACGTGGCATCCTGc                                                                                                                                                                                                              >  1:294342/1‑239 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTg                                                                                                                                                                                                               >  1:436198/1‑238 (MQ=255)
                            aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAACGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATAAACCT‑GTTCCATTAAAGGTAACGTGGCACCCTg                                                                                                                                                                                                               >  1:56631/1‑238 (MQ=255)
                                   aTCAACCTGGGCAAACAATC‑CCGGACCAAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCtt                                                                                                                                                                                                       <  1:683396/239‑1 (MQ=255)
                                         cTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGGAACGTGGCATCCTGCCAGCCTTCGGCa                                                                                                                                                                                                  >  1:464027/1‑238 (MQ=255)
                                           cGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTTTGCCGCTTCTTCCATCGCCCCGGTCAA‑GCGGTCACGTTCCAGA‑GCCGGGTCCAGCCT‑TTATGCCTGATAAACCT‑GGTCCATTAAAGACACCGAGGCAACCGGCCACCCTTCGACCATc                                                                                                                                                                                               >  2:406947/1‑239 (MQ=255)
                                            ggCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAACCTCTTTTTCCGCTTCTTCCCGCGCCCCGGTCAG‑GCGTTCGCGTTCCAGA‑TCCGGATCCACCGT‑TGATGCCTGATACCCCC‑TTTCCATTCAAGGTACCGCGCCACCCTCCCAccccctggccacg                                                                                                                                                                                               >  2:85209/1‑228 (MQ=255)
                                            ggCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCCGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATc                                                                                                                                                                                               >  1:53215/1‑238 (MQ=255)
                                            ggCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATc                                                                                                                                                                                               >  1:217390/1‑238 (MQ=255)
                                                                        gtCCGAAAGCGGGGGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCCCCGGCGGCAAAGCGTTTGCTGTAGCTGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATgcg                                                                                                                                                                    <  1:510801/238‑1 (MQ=255)
                                                                                  gggggCGAGTAAGAATCGAAAATGGCCGCCGTTTCTTTTTGGGCCCCGGCGGCAAAGCCTTTGCTGGAGCGGCGAAACTGGTTTGCCGGTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                                                                                                                          <  2:536349/234‑1 (MQ=255)
                                                                                   ggTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCC‑GTTCCCTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGCCAATCCCCACACCCGGTGCTCCCGGTAATTGCCGGATTctcg                                                                                                                                                          >  1:407001/1‑237 (MQ=255)
                                                                                    gTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                                                                                                                          >  1:177079/1‑236 (MQ=255)
                                                                                           tAAAGTCCGAAAATCGCCCCCGTTTTTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                                                                                                                                 <  2:294342/238‑1 (MQ=255)
                                                                                                 ccGAAAATAGCCGCCGTTTCTTGTTGGGCACCGGCGGCAAAGCGTTGGCTGTAGCGGCGAAACTCGGTTGCCGCTGCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGGAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                                                                                                                            <  1:85209/236‑1 (MQ=255)
                                                                                                         aGCCCCCGTTTCTTTTTGGGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                                                                                                                  <  1:1062448/239‑1 (MQ=255)
                                                                                                                            gcgcgggggggaaaGCGTTTGCTGTAGCGGCGAACCCCGTTTCCCTCTTCTCCCAGCGCCCCGTCCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCGGCGT‑TGAGGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCGGCCAGCCTTCGGCAATCGCCACGCCAGGTGCGGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTgaga                                                                                                              <  2:56631/228‑1 (MQ=255)
                                                                                                                            gccccGGGGGCAAAGCGTTTGTTGTGGGGGCCAAATCCTTTTTCCGTTTCTTCCAGCGCCCCGGTGAG‑TCGGTCGGGTGCAAGA‑GCCGGATCCAGGGT‑TGATGCCTGATACACCA‑GTTCCATTAAAGGTACCGTGGCATCCTGCCAGCCTTCGGCAATTGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTgaga                                                                                                              <  2:1139700/236‑1 (MQ=255)
                                                                                                                                       aaaGCGTTTGCTGGAGCGGCCAATCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTGCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATCCCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCGTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCt                                                                                                    <  1:654057/238‑1 (MQ=255)
                                                                                                                                            gtttgtgttcgcgacgAACATCTTTTGCCGCTTTTTCCGGCCCCCCGGTCAG‑CCCTTCGTGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACGCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTGCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATc                                                                                              <  2:562472/229‑1 (MQ=255)
                                                                                                                                             tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCCGCGCCCCGGTCAG‑TCGTTCGAGGTCCCGA‑GCCGCATCCAGCGT‑TGCTGCCTGATACCCCT‑GTTCCATTAAAGGCAACGTGGCATCCTGCCAGCCTTCCGCCATCGCCACACCAAGAGCTGACGGTCATGCCGGGAATCTCGGCTGCCGATATTGCCCCACCAACGCAGACAACTGCT‑ACGGagaaagaaaaacagcat                                                                                                >  2:5396/1‑227 (MQ=255)
                                                                                                                                             tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCCG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGTGATTGCCCAAACAAGCGAGTCAACTGCG‑ACTGAGAAAGAATAACTGCCa                                                                                                >  1:382798/1‑236 (MQ=255)
                                                                                                                                             tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGGAAAGGGGCATCCTCCCAGCCTTCCGCAATCACCACACCCGGTGTTGCAGGCAATGAGGGGATTCGGGGCTGGCGATATTGCCCAACCCAGGCAGTC‑ACCGCGCACTGagaaagaaaaacggccc                                                                                                >  2:939845/1‑227 (MQ=255)
                                                                                                                                                     aGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGGGGCa                                                                                     >  1:1219787/1‑239 (MQ=255)
                                                                                                                                                          cgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCGGCAACCTGGGTGGCAAGGGt                                                                                >  1:134949/1‑239 (MQ=255)
                                                                                                                                                                 cGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GGCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCagg                                                                         <  1:380559/239‑1 (MQ=255)
                                                                                                                                                                   ttttCCGTTTCTCCCGGCGCCCCGGTCCG‑TCTGTCGCGTTCCGGA‑GCCGGATCCGGCGT‑TGAGGCCTGATCCCCCT‑GTTCAATTAAAGGTAACGTGGCAGCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                       <  2:924538/239‑1 (MQ=255)
                                                                                                                                                                   tttGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                       <  1:1209132/239‑1 (MQ=255)
                                                                                                                                                                    ttGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGCGAAAGACAAGCTGCCACCTGGGGGGCAGGCGCAACCaggaag                                                                      >  1:1104803/1‑239 (MQ=255)
                                                                                                                                                                    ttGCCGCTTCTTCCAGCGCCCCGGTCAA‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaag                                                                      >  1:728491/1‑239 (MQ=255)
                                                                                                                                                                      gCCCCTTCTTCCCCCACCCCGGTGAC‑CCCTTCGCGGTCCAGG‑CGCGGACCCCGAGT‑TTGTGCCTGA‑ACCCCTGGTTCGTTTAAAGGTAACGGGGCACTCTCCAAGCTTGCGGAAATCGGCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCCACGGGCGATATTGCCCAAACAAGGCATTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGTTGGCAAGCGTCACCaggaagga                                                                    <  2:1161628/239‑1 (MQ=255)
                                                                                                                                                                            tcttcCAGCGCCCCGGTCAG‑TCGTTTGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGACGTCAACTCCG‑ACTGAGAAAGAGTAGCTTCCATCTGGGTGGCAACGGTGACCaggaaggtgtctc                                                               >  1:533069/1‑232 (MQ=255)
                                                                                                                                                                              tcccAGCCCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCAc                                                            <  1:1112422/237‑1 (MQ=255)
                                                                                                                                                                                ccAGCGCCCCGGGCAG‑TCGTTCGCGTTCAAGG‑GCGGGATCCCAGGT‑TGAGGCCTGATACACCT‑GTCCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACtt                                                          <  2:464027/239‑1 (MQ=255)
                                                                                                                                                                                 cAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACttt                                                         <  1:159466/239‑1 (MQ=255)
                                                                                                                                                                                   gcgcCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACGCCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                       <  1:130456/239‑1 (MQ=255)
                                                                                                                                                                                    ggcCCCCGGCAGTTCGGTCGGGTCCAAGA‑GCCGGGTCCGGTGT‑TGATGCCTGAAACACCC‑TGTCCATTAAAGGGACCGTGGCATCCTGCCGGCCTGCGGCAAGCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                       <  2:53215/238‑1 (MQ=255)
                                                                                                                                                                                       cccGGTACG‑TCGTTGGCTTTCCAGG‑GCGGGATCCAGCGT‑TGAGTCCTGTGACCCCG‑GGTGCATTAAGGGTAACGTGGCATCCTGCCGGCCTCCGGCATACGCCACCCCAGGTGCTGGCGGTAATGCGCGGATTCGCGACTGGCGATATTGCCCAATCGAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTGCACTTTCGt                                                      <  2:382798/236‑1 (MQ=255)
                                                                                                                                                                                                acggttgccgTTTCAGAGGCAGGAACCCGCGT‑TAATGCCTGATCCCCCT‑GTTCCATTAAAGGGAACGGGGCACCCTGCCAGCCTTCTGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCGCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCTCCAGGAAGGATTCTTCACTTTCGTCATACTGgcg                                            <  2:709286/231‑1 (MQ=255)
                                                                                                                                                                                                  cGTTCGGGTTCCAGA‑GCCGGATCCAGCGT‑GGATGCCTGATCCCCCT‑GTTCCATTAAGGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCa                                          <  2:882410/238‑1 (MQ=255)
                                                                                                                                                                                                      cgcgTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGTCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCACCTGCG‑ACTGAGAAAGAATAGCTGCCATCTTGGTGGCAAGCGTCACCAGGAAGGCTTCTTCCCTTTCGTCATACTGGCCCACctc                                      >  1:1129020/1‑238 (MQ=255)
                                                                                                                                                                                                      cgcgTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGATTAGCTGCCATCTGGGTGGCAAGCGCCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                     >  1:43456/1‑239 (MQ=255)
                                                                                                                                                                                                         gTTCCAGA‑CCCGGCTCCGGCTT‑TGATGCTTGATACCCCT‑GTTCCATTAAAGGTAACGGGGCATCCGGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGt                                  <  2:436198/239‑1 (MQ=255)
                                                                                                                                                                                                            ccAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTAACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTg                                >  1:882410/1‑238 (MQ=255)
                                                                                                                                                                                                             cgga‑gCCGGCTCCCGGTT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCt                              <  1:509362/237‑1 (MQ=255)
                                                                                                                                                                                                                   ccGGATCCAGCGT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                          <  1:455317/238‑1 (MQ=255)
                                                                                                                                                                                                                    cGGATCCAGCGG‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                          <  1:5263/237‑1 (MQ=255)
                                                                                                                                                                                                                        tCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGGATAACTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGACATACTGGCGCAACCCCCGTTGCGGTacaaacac                    >  2:646037/1‑233 (MQ=255)
                                                                                                                                                                                                                        tCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTAccacca                     <  1:742631/238‑1 (MQ=255)
                                                                                                                                                                                                                         ccAGCGT‑TGATTCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAg                    <  1:406947/238‑1 (MQ=255)
                                                                                                                                                                                                                          cAGCGT‑TGATGCCTTATACGCCG‑GTTCCATTAAAGTTAACGTGGCATCCTGCCAGCCTTCGGAAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                  <  1:1308562/239‑1 (MQ=255)
                                                                                                                                                                                                                            gCGT‑TGGTGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaca                <  1:435556/239‑1 (MQ=255)
                                                                                                                                                                                                                             cGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaccc               >  1:85319/1‑237 (MQ=255)
                                                                                                                                                                                                                             cGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAGTTTTCCCTTTCGTCATACTGGCGCCACTCCCGTTGCGGGACCCCCACAacac               >  2:998392/1‑239 (MQ=255)
                                                                                                                                                                                                                                   aTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATCTTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAAAAGCTGCCATCTGGGAGGCAGGCGTCACCAGGAAGGATTCTTCACTTTCGTCATAAAGGGGCAACTCTCGTTGCCGTACCACCAGTA‑ACCaaacc         >  2:1172694/1‑236 (MQ=255)
                                                                                                                                                                                                                                   aTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAct                                       >  1:946905/1‑210 (MQ=255)
                                                                                                                                                                                                                                   aTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTCCTGTACCACCAGACCACCAAGc          >  1:546545/1‑239 (MQ=255)
                                                                                                                                                                                                                                        tGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGGCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGACCCACCAGAAAACCaacaaacgg     >  2:348963/1‑232 (MQ=255)
                                                                                                                                                                                                                                            ccacCT‑GTTCCATTAAGGGTAACGTGGCCTCTGGCCAGGCTTCGGGAAGCGCCACAACAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCATACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCg  <  1:796196/237‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
ACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCG  >  NC_000913/2967611‑2968080

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