Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 92,843 A→G P477P (CCA→CCG)  ftsI → peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091392,8430AG93.2% 164.8 / ‑1.7 59P477P (CCA→CCG) ftsIpeptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base G (30/25);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

AACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCAT  >  NC_000913/92599‑93073
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
aaCGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTAAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGgg‑cggc                                                                                                                                                                                                                                    >  2:414171/1‑251 (MQ=255)
   gaagtagCTTATGTCCTCAACAACAACGGTGGACTGACATAGAGAGGGCCACATTCTCTTTCGGCTACGAGCTAATGGTCACACCGTTACAGTTAGAGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCCAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGACCGCACCGAGGTGCACATGATGGCAAACGAGGCGCGACCAgcgggggac                                                                                                                                                                                                                                  >  2:452830/7‑243 (MQ=255)
            tataTCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGATACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGa                                                                                                                                                                                                                          >  1:113239/1‑249 (MQ=255)
                 ccTCAAAAACAACGGTGGTCTGAAATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAggcgg                                                                                                                                                                                                                    <  1:96474/250‑1 (MQ=255)
                  cacAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCGTCCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCGCCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGgcggcg                                                                                                                                                                                                                  <  2:158866/249‑1 (MQ=255)
                      aaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTa                                                                                                                                                                                                              <  1:41040/251‑1 (MQ=255)
                       aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCCACTATCGGCAGCTACGGCCTTTATCGCCCCCTGTCGCTTACCAAAGTTCACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGGGCATTTGATGGAAAGCGGGGCGCTACCGGG‑ATGCGGCGGCGTGAcggcggggctta                                                                                                                                                                                                              >  2:90969/1‑244 (MQ=255)
                            acGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGACCGCACCGTGGTGCCTATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCt                                                                                                                                                                                                        >  2:446989/1‑251 (MQ=255)
                                 ggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgt                                                                                                                                                                                                   >  1:364789/1‑251 (MQ=255)
                                             atagGGCCACCCTCTCTTTCGGCTACGGGCTAATGGTAACCCCGTTACAGTTAGCGCGAGTCTACGCAACTATCGTCAGCTACGGCATTTAGCGCCCACTGGCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa                                                                                                                                                                                       <  2:352400/249‑1 (MQ=255)
                                              gagGGCCACCTTCTCTTTCGGCTACGCGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTGAACGTGTCTGCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaaa                                                                                                                                                                                      <  1:19411/251‑1 (MQ=255)
                                                  gCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTATACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGCCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCgg                                                                                                                                                                                  >  2:406496/1‑251 (MQ=255)
                                                   ccACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGAGAGTCTAAGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCACGGTGAAAGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGAGCCTATGATGGAAAGCCTGGCGCT‑CCAGGGAGGCGGCGGCGTGAAGGCGGCGATTAAAGGAGACCGTATCACCCTTAAAAcagg                                                                                                                                                                                  >  2:133663/1‑247 (MQ=255)
                                                     cccTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAc                                                                                                                                                                               <  1:414171/250‑1 (MQ=255)
                                                                  gCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTc                                                                                                                                                                <  2:362128/251‑1 (MQ=255)
                                                                   cTACGGGCTACTGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAGGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAGAGGCTATCGTATCGCCCTTAAAACCGGTAC‑C‑GCGAAACAGGTCg                                                                                                                                                               >  2:4484/1‑251 (MQ=255)
                                                                   cTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCg                                                                                                                                                               >  1:314107/1‑251 (MQ=255)
                                                                       gggCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGcc                                                                                                                                                           <  1:363148/251‑1 (MQ=255)
                                                                       gggCTAAGTGTAACACCGGTACTGTTAGCGCGAGGCTACGCAACTATCGTCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGcc                                                                                                                                                           <  2:20532/251‑1 (MQ=255)
                                                                         gCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTCTCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCACGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGAATAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑A‑GCGAAAAAGGTCGGGCCgg                                                                                                                                                         >  1:62683/1‑251 (MQ=255)
                                                                           tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGCTTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGAGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTAGCGTATCGCCATTAAAACCGG‑AC‑CAGCGAAAAAGGGCGGGCCGGAc                                                                                                                                                       >  2:1803/1‑251 (MQ=255)
                                                                           tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGCATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTCAAACCGGTAC‑C‑GCGGAACAGGGCGGGCCgtcc                                                                                                                                                       >  2:206291/1‑248 (MQ=255)
                                                                              tGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCACACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGGGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTGCCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGCAC‑C‑GCGAAAAAgggcgggcagggaggg                                                                                                                                                    >  2:175186/1‑243 (MQ=255)
                                                                                  aaCACCGTTTCAGTTAGCGCGAGTCTACGCAACTACCGGCAGCTACGGCATTTATCGCCCACTGTCGATTTCCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCt                                                                                                                                                <  2:354930/251‑1 (MQ=255)
                                                                                  aaCACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGc                                                                                                                                                 <  1:160286/250‑1 (MQ=255)
                                                                                   acacCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAATTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCTGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTAC‑C‑GCGAAAAAGGTAGGGCCCGACgggccctc                                                                                                                                               >  2:179159/1‑245 (MQ=255)
                                                                                    caccgtttcaggtAGCGCGAGTCTACGCACCTATCGGCAGCTACGGCAGTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGGGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGTTCGGGCCGGACGGGCGCTAc                                                                                                                                              <  2:416168/239‑1 (MQ=255)
                                                                                      ccgttACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACAt                                                                                                                                            <  1:203751/248‑1 (MQ=255)
                                                                                              gTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGAATTTATCGCACACTGGCGATTACCACAGTTGGCCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAATGCGTGGCGCTACAGGG‑CGGAGACGGCGTGGAGGCGGCGATTAAGGGCTATCGTATCGCCATTAGAACCGGTACAC‑G‑GAAAAAGGTCGGAACGGAAGGTCGGTACCTCaataaat                                                                                                                                    >  2:406084/1‑251 (MQ=255)
                                                                                                tAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGGCGGGCCGGACGGTCGCTACATCAATAAatat                                                                                                                                  >  2:98664/1‑251 (MQ=255)
                                                                                                                  cTATCGGCAGCTACGGCATTTATCGACCACTGTCGATTACCACAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCTCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGc                                                                                                                >  1:133950/1‑251 (MQ=255)
                                                                                                                            cTACGGCATTTATCGCCCACTGTCGATTCCCAAAGTTGACCCCCCGGTTCCCGGTGCACGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTGCGTATGATGGAAAGCGTGGCGCTACCGGA‑CGGCGCCGGCGTGAAGGCGGCGATTAACGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGACCACGGTCGAGCCGGACGGTCGTTACATCATTATAAACATTGATTATACCGCAGGCGTTGCGCCTg                                                                                                      >  2:107860/1‑251 (MQ=255)
                                                                                                                             tACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCg                                                                                                                                                               >  2:38795/1‑193 (MQ=255)
                                                                                                                             tACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCg                                                                                                                                                               <  1:38795/193‑1 (MQ=255)
                                                                                                                                  cATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCTCTACCGGT‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCAGTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTc                                                                                                <  1:305758/251‑1 (MQ=255)
                                                                                                                                              aCTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAAACCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTcg                                                                                    >  2:468299/1‑251 (MQ=255)
                                                                                                                                               atGTCGATTACCACAGTTGACCCCCCGGATCCCGCGGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCACATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTACAGGCTATCGAATCGCCATTAGAACCGGAAC‑C‑GCGTGACGGGTAGGGCCCGACGGGGGCTACATCCGTAAATATATAGGTTAAAACGCAGGCCATGCGccgtggggtccacagcacgtcac                                                                                   >  2:259890/2‑230 (MQ=255)
                                                                                                                                                 gTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCgcgc                                                                                 <  1:243068/251‑1 (MQ=255)
                                                                                                                                                           aaaGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgt                                                                       >  1:436574/1‑251 (MQ=255)
                                                                                                                                                                      cccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                                            <  1:361780/251‑1 (MQ=255)
                                                                                                                                                                       ccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                                           >  1:144835/1‑251 (MQ=255)
                                                                                                                                                                        cGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTTTCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTTCTTATCCCGCAGGCGTTGCGCCTGCGCGTCCGCCGCGCATCGCGCTGGTTGTTGTTATCAACGATCCg                                                          >  2:111370/1‑251 (MQ=255)
                                                                                                                                                                            tCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAgg                                                      >  1:283990/1‑251 (MQ=255)
                                                                                                                                                                                ggTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCggg                                                  <  1:4484/251‑1 (MQ=255)
                                                                                                                                                                                         gtCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                                           <  2:24739/185‑1 (MQ=255)
                                                                                                                                                                                         gtCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                                           >  1:24739/1‑185 (MQ=255)
                                                                                                                                                                                             tCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGTCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                                        <  1:198396/248‑1 (MQ=255)
                                                                                                                                                                                                cGGAAGCCAGTGGCCGCACCGTGGGGAATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGTCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTAGACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc                                  <  2:174636/251‑1 (MQ=255)
                                                                                                                                                                                                  gAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTAGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAATGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGACAGCCGCGCTTCGCGCTGGTTGTTGTTTTCAACGATCCGCAGGCggggacataccaacgccgggg                                >  2:338162/1‑233 (MQ=255)
                                                                                                                                                                                                   aaTCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGcc                               <  1:206291/251‑1 (MQ=255)
                                                                                                                                                                                                       cATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGACTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaaatcttcggaggagcccgtt                           >  1:19009/1‑232 (MQ=255)
                                                                                                                                                                                                           gTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAttataaggagggggcgttttccg                       >  1:405883/1‑229 (MQ=255)
                                                                                                                                                                                                           gTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcgc                       <  1:394048/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggcaaaggtgggAAATCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTAACGCCATTAAAACCGGTCC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGTAGGCGGGTAAATACTACGGCGGCGTCGTTTCCGCGCGGGTCTTTGGt          <  2:58602/241‑1 (MQ=255)
                                                                                                                                                                                                                           atatGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAAAACTACGGCGGGGCCCGTTCCACGCCCGTCttttgggcc       >  1:347024/1‑245 (MQ=255)
                                                                                                                                                                                                                            tatGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACCGCGGCGCCGGTTTCTCGCCgggcgtttgggacc      >  1:458412/1‑239 (MQ=255)
                                                                                                                                                                                                                             atgatgGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCca      <  1:98664/250‑1 (MQ=255)
                                                                                                                                                                                                                               gatgGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAATACCACGCCGGCGCCGGTTCCGCGCCgggcgtttggggcatcc   >  1:194240/1‑236 (MQ=255)
                                                                                                                                                                                                                                atgGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGGTTTCGCGCCGGTCttttggtcgctcat  >  1:288161/1‑240 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
AACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCAT  >  NC_000913/92599‑93073

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