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→CCGftsIpeptidoglycan 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

AACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCAT  >  NC_000913/92599‑93073
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
aaCGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTAAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGggcggc                                                                                                                                                                                                                                    >  2:414171/1‑251 (MQ=255)
   gaagtagCTTATGTCCTCAACAACAACGGTGGACTGACATAGAGAGGGCCACATTCTCTTTCGGCTACGAGCTAATGGTCACACCGTTACAGTTAGAGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCCAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGACCGCACCGAGGTGCACATGATGGCAAACGAGGCGCGACCAgcgggggac                                                                                                                                                                                                                                  >  2:452830/7‑243 (MQ=255)
            tataTCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGATACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGa                                                                                                                                                                                                                          >  1:113239/1‑249 (MQ=255)
                 ccTCAAAAACAACGGTGGTCTGAAATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAggcgg                                                                                                                                                                                                                    <  1:96474/250‑1 (MQ=255)
                  cacAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCGTCCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCGCCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGgcggcg                                                                                                                                                                                                                  <  2:158866/249‑1 (MQ=255)
                      aaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTa                                                                                                                                                                                                              <  1:41040/251‑1 (MQ=255)
                       aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCCACTATCGGCAGCTACGGCCTTTATCGCCCCCTGTCGCTTACCAAAGTTCACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGGGCATTTGATGGAAAGCGGGGCGCTACCGGGATGCGGCGGCGTGAcggcggggctta                                                                                                                                                                                                              >  2:90969/1‑244 (MQ=255)
                            acGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGACCGCACCGTGGTGCCTATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCt                                                                                                                                                                                                        >  2:446989/1‑251 (MQ=255)
                                 ggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgt                                                                                                                                                                                                   >  1:364789/1‑251 (MQ=255)
                                             atagGGCCACCCTCTCTTTCGGCTACGGGCTAATGGTAACCCCGTTACAGTTAGCGCGAGTCTACGCAACTATCGTCAGCTACGGCATTTAGCGCCCACTGGCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa                                                                                                                                                                                       <  2:352400/249‑1 (MQ=255)
                                              gagGGCCACCTTCTCTTTCGGCTACGCGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTGAACGTGTCTGCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaaa                                                                                                                                                                                      <  1:19411/251‑1 (MQ=255)
                                                  gCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTATACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGCCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCgg                                                                                                                                                                                  >  2:406496/1‑251 (MQ=255)
                                                   ccACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGAGAGTCTAAGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCACGGTGAAAGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGAGCCTATGATGGAAAGCCTGGCGCTCCAGGGAGGCGGCGGCGTGAAGGCGGCGATTAAAGGAGACCGTATCACCCTTAAAAcagg                                                                                                                                                                                  >  2:133663/1‑247 (MQ=255)
                                                     cccTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAc                                                                                                                                                                               <  1:414171/250‑1 (MQ=255)
                                                                  gCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTc                                                                                                                                                                <  2:362128/251‑1 (MQ=255)
                                                                   cTACGGGCTACTGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAGGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAGAGGCTATCGTATCGCCCTTAAAACCGGTACCGCGAAACAGGTCg                                                                                                                                                               >  2:4484/1‑251 (MQ=255)
                                                                   cTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCg                                                                                                                                                               >  1:314107/1‑251 (MQ=255)
                                                                       gggCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGcc                                                                                                                                                           <  1:363148/251‑1 (MQ=255)
                                                                       gggCTAAGTGTAACACCGGTACTGTTAGCGCGAGGCTACGCAACTATCGTCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGcc                                                                                                                                                           <  2:20532/251‑1 (MQ=255)
                                                                         gCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTCTCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCACGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGAATAAAGGCTATCGTATCGCCATTAAAACCGGTACAGCGAAAAAGGTCGGGCCgg                                                                                                                                                         >  1:62683/1‑251 (MQ=255)
                                                                           tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGCTTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGAGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTAGCGTATCGCCATTAAAACCGGACCAGCGAAAAAGGGCGGGCCGGAc                                                                                                                                                       >  2:1803/1‑251 (MQ=255)
                                                                           tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGCATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTCAAACCGGTACCGCGGAACAGGGCGGGCCgtcc                                                                                                                                                       >  2:206291/1‑248 (MQ=255)
                                                                              tGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCACACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGGGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTGCCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGCACCGCGAAAAAgggcgggcagggaggg                                                                                                                                                    >  2:175186/1‑243 (MQ=255)
                                                                                  aaCACCGTTTCAGTTAGCGCGAGTCTACGCAACTACCGGCAGCTACGGCATTTATCGCCCACTGTCGATTTCCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCt                                                                                                                                                <  2:354930/251‑1 (MQ=255)
                                                                                  aaCACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGc                                                                                                                                                 <  1:160286/250‑1 (MQ=255)
                                                                                   acacCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAATTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTACCGGGCGGCGGCGGCGTGAAGGCTGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTACCGCGAAAAAGGTAGGGCCCGACgggccctc                                                                                                                                               >  2:179159/1‑245 (MQ=255)
                                                                                    caccgtttcaggtAGCGCGAGTCTACGCACCTATCGGCAGCTACGGCAGTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGGGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGTTCGGGCCGGACGGGCGCTAc                                                                                                                                              <  2:416168/239‑1 (MQ=255)
                                                                                      ccgttACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACAt                                                                                                                                            <  1:203751/248‑1 (MQ=255)
                                                                                              gTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGAATTTATCGCACACTGGCGATTACCACAGTTGGCCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAATGCGTGGCGCTACAGGGCGGAGACGGCGTGGAGGCGGCGATTAAGGGCTATCGTATCGCCATTAGAACCGGTACACGGAAAAAGGTCGGAACGGAAGGTCGGTACCTCaataaat                                                                                                                                    >  2:406084/1‑251 (MQ=255)
                                                                                                tAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGGCGGGCCGGACGGTCGCTACATCAATAAatat                                                                                                                                  >  2:98664/1‑251 (MQ=255)
                                                                                                                  cTATCGGCAGCTACGGCATTTATCGACCACTGTCGATTACCACAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCTCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGc                                                                                                                >  1:133950/1‑251 (MQ=255)
                                                                                                                            cTACGGCATTTATCGCCCACTGTCGATTCCCAAAGTTGACCCCCCGGTTCCCGGTGCACGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTGCGTATGATGGAAAGCGTGGCGCTACCGGACGGCGCCGGCGTGAAGGCGGCGATTAACGGCTATCGTATCGCCATTAAAACCGGTACCGCGACCACGGTCGAGCCGGACGGTCGTTACATCATTATAAACATTGATTATACCGCAGGCGTTGCGCCTg                                                                                                      >  2:107860/1‑251 (MQ=255)
                                                                                                                             tACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCg                                                                                                                                                               >  2:38795/1‑193 (MQ=255)
                                                                                                                             tACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCg                                                                                                                                                               <  1:38795/193‑1 (MQ=255)
                                                                                                                                  cATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCTCTACCGGTCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCAGTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTc                                                                                                <  1:305758/251‑1 (MQ=255)
                                                                                                                                              aCTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAAACCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTcg                                                                                    >  2:468299/1‑251 (MQ=255)
                                                                                                                                               atGTCGATTACCACAGTTGACCCCCCGGATCCCGCGGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCACATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTACAGGCTATCGAATCGCCATTAGAACCGGAACCGCGTGACGGGTAGGGCCCGACGGGGGCTACATCCGTAAATATATAGGTTAAAACGCAGGCCATGCGccgtggggtccacagcacgtcac                                                                                   >  2:259890/2‑230 (MQ=255)
                                                                                                                                                 gTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCgcgc                                                                                 <  1:243068/251‑1 (MQ=255)
                                                                                                                                                           aaaGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgt                                                                       >  1:436574/1‑251 (MQ=255)
                                                                                                                                                                      cccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                                            <  1:361780/251‑1 (MQ=255)
                                                                                                                                                                       ccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                                           >  1:144835/1‑251 (MQ=255)
                                                                                                                                                                        cGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTTTCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTTCTTATCCCGCAGGCGTTGCGCCTGCGCGTCCGCCGCGCATCGCGCTGGTTGTTGTTATCAACGATCCg                                                          >  2:111370/1‑251 (MQ=255)
                                                                                                                                                                            tCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAgg                                                      >  1:283990/1‑251 (MQ=255)
                                                                                                                                                                                ggTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCggg                                                  <  1:4484/251‑1 (MQ=255)
                                                                                                                                                                                         gtCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                                           <  2:24739/185‑1 (MQ=255)
                                                                                                                                                                                         gtCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                                           >  1:24739/1‑185 (MQ=255)
                                                                                                                                                                                             tCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGTCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                                        <  1:198396/248‑1 (MQ=255)
                                                                                                                                                                                                cGGAAGCCAGTGGCCGCACCGTGGGGAATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGTCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTAGACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc                                  <  2:174636/251‑1 (MQ=255)
                                                                                                                                                                                                  gAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTAGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAATGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGACAGCCGCGCTTCGCGCTGGTTGTTGTTTTCAACGATCCGCAGGCggggacataccaacgccgggg                                >  2:338162/1‑233 (MQ=255)
                                                                                                                                                                                                   aaTCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGcc                               <  1:206291/251‑1 (MQ=255)
                                                                                                                                                                                                       cATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGACTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaaatcttcggaggagcccgtt                           >  1:19009/1‑232 (MQ=255)
                                                                                                                                                                                                           gTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAttataaggagggggcgttttccg                       >  1:405883/1‑229 (MQ=255)
                                                                                                                                                                                                           gTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcgc                       <  1:394048/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggcaaaggtgggAAATCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTAACGCCATTAAAACCGGTCCCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGTAGGCGGGTAAATACTACGGCGGCGTCGTTTCCGCGCGGGTCTTTGGt          <  2:58602/241‑1 (MQ=255)
                                                                                                                                                                                                                           atatGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAAAACTACGGCGGGGCCCGTTCCACGCCCGTCttttgggcc       >  1:347024/1‑245 (MQ=255)
                                                                                                                                                                                                                            tatGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACCGCGGCGCCGGTTTCTCGCCgggcgtttgggacc      >  1:458412/1‑239 (MQ=255)
                                                                                                                                                                                                                             atgatgGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCca      <  1:98664/250‑1 (MQ=255)
                                                                                                                                                                                                                               gatgGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAATACCACGCCGGCGCCGGTTCCGCGCCgggcgtttggggcatcc   >  1:194240/1‑236 (MQ=255)
                                                                                                                                                                                                                                atgGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGGTTTCGCGCCGGTCttttggtcgctcat  >  1:288161/1‑240 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
AACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCAT  >  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: