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,8430AG98.6% 230.0 / ‑4.1 71P477P (CCA→CCGftsIpeptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (39/31);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

GCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTAT  >  NC_000913/92608‑93084
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
gctttaaaccTCAAAAACAACGGTGGTATGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTAGCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTg                                                                                                                                                                                                                                     <  2:584314/243‑1 (MQ=255)
gCTTATATCCTCAAAAACAACCGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCACAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGACTCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTg                                                                                                                                                                                                                                     >  1:851345/1‑251 (MQ=255)
   tatttCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTAGCGGGCGGCGGCGGCGTGaa                                                                                                                                                                                                                                   <  2:127651/246‑1 (MQ=255)
        ccTCAAAAACCAAGGTGGCCTGCCATAGAGAGTGCCACCTTCTCTTCCGACTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTAGATTACCAAAGTTAAACCCACGGGTCCCGGTGACCGAGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCCTATGTTGGCACGCGTGACGCTACCGGGCGGCGCCGGAGTGCAggcggc                                                                                                                                                                                                                             >  2:73225/1‑251 (MQ=255)
           cAAAAACAACGGTGGTCGGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGGTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGAt                                                                                                                                                                                                                          <  2:149414/251‑1 (MQ=255)
            aaaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGAtt                                                                                                                                                                                                                         <  1:503274/251‑1 (MQ=255)
              aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGGAGGCGGCGATTa                                                                                                                                                                                                                        >  2:364626/1‑250 (MQ=255)
              aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTaa                                                                                                                                                                                                                       <  2:649071/251‑1 (MQ=255)
              aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGGGAAGGCGGCGATTa                                                                                                                                                                                                                        >  1:351947/1‑250 (MQ=255)
                 caacGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAgg                                                                                                                                                                                                                    >  1:657004/1‑251 (MQ=255)
                  aacGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGc                                                                                                                                                                                                                   >  1:618188/1‑251 (MQ=255)
                   acGGGGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGAAACACAGTTACAGTTAGCGAGAGTCTCCGCAACTATCGGCAGCTACGGCATTTAACGCCCACTGTCGATTACAAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCGTTGTCCGACCCGTGGAGCACATGGTGGAAAGAGTGGGGCTACCGAGCGGCGGCGGCGGGAAGTGGGGGGTTAAAGGCt                                                                                                                                                                                                                  >  2:153829/1‑251 (MQ=255)
                       tggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcg                                                                                                                                                                                                              >  1:729156/1‑251 (MQ=255)
                       tggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcg                                                                                                                                                                                                              >  1:499693/1‑251 (MQ=255)
                        ggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcg                                                                                                                                                                                                              <  1:461323/250‑1 (MQ=255)
                            tGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgtatcg                                                                                                                                                                                                        >  1:533047/1‑251 (MQ=255)
                             gACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgtat                                                                                                                                                                                                           <  1:448122/248‑1 (MQ=255)
                                 aagagagGGCCCCCTTTCCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACTGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCAtt                                                                                                                                                                                                   <  2:618188/250‑1 (MQ=255)
                                   gagagGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTa                                                                                                                                                                                                  >  1:167279/1‑250 (MQ=255)
                                    agagGGCCAACTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTGCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa                                                                                                                                                                                                <  1:79128/251‑1 (MQ=255)
                                     atgggccgacttatctgtcgtcgaCGGGCTAAGGGTGACGCCGTTCCAGTTAGCGTGAGTCTACGCAACTTTCGGCTGCTACTGCATTTATCGCCCACTGTCGATTACCAAAGTTGTCCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCTGTGTCCACGCCGTGGTACATATGATGGAAATCGTTGCGCTACCGGGCGGCGGAAGCGTGAAGGCGGCGATTACAGGCTAGCGTAGCTCCATTaaaa                                                                                                                                                                                               <  2:675989/228‑1 (MQ=255)
                                      agggccgcctgttctTTCGGCTACGTGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTTCGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAAc                                                                                                                                                                                              <  2:147392/239‑1 (MQ=255)
                                       gggCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGCACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGGGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTCAAAcc                                                                                                                                                                                             >  2:447424/1‑251 (MQ=255)
                                        ggCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGTCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCg                                                                                                                                                                                            >  1:235996/1‑251 (MQ=255)
                                            aCCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAc                                                                                                                                                                                        >  1:815671/1‑251 (MQ=255)
                                              cTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTAACGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTAACGTATCGCCATTAAAACCGGTAc                                                                                                                                                                                        >  1:421229/1‑249 (MQ=255)
                                                  tctTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGa                                                                                                                                                                                   >  1:778154/1‑250 (MQ=255)
                                                       cGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAg                                                                                                                                                                              <  1:364626/250‑1 (MQ=255)
                                                         gCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTACCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTc                                                                                                                                                                           >  1:537830/1‑251 (MQ=255)
                                                           tACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGAAGCTACGGCATTTATCGCCCACTGTCGATTACCACAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGATACCGCGAAAAAGGTCgg                                                                                                                                                                         >  1:831218/1‑251 (MQ=255)
                                                             cGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTaaa                                                                                                                                                                                                                      <  2:714135/205‑1 (MQ=255)
                                                             cGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTaaa                                                                                                                                                                                                                      >  1:714135/1‑205 (MQ=255)
                                                                   aaTGGTAACACCGTTACAGTTAGCGCGAGTCTACGCGACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGAc                                                                                                                                                                  <  1:648208/250‑1 (MQ=255)
                                                                      ggTAACACCGTTACAGTTAGCGCGAGGCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCTCTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGt                                                                                                                                                               <  2:868767/250‑1 (MQ=255)
                                                                             ccgttACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACAt                                                                                                                                                       <  1:378726/248‑1 (MQ=255)
                                                                             ccgttACAGTTAGCGCGAGTCGACGCAACTTGCGGCATCTACGTCATTGATCTCCCACTGTCGATGACCAAAGTTGACCCCCCGGTTCCCGGTGGACGTGGCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGGGGCGGCGGCGGGAATGCGGCGATTAAAGGCGATCGGATCGCCATTAAGACCGGGGCCGCGAAAAATGGCGGTCCGGACGGTCGCTACAt                                                                                                                                                       <  2:808420/248‑1 (MQ=255)
                                                                                       tAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCTCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTTAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAATGCTATCGTAACGCCATTAAAACCGGTACCGCTAAAAAGTTCGGGCCGGACGGTCGCTACATCaataaata                                                                                                                                              <  1:410947/250‑1 (MQ=255)
                                                                                          cgcgAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAAAGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCCTTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGCAGGTCGGTACATCACTAAATATAtt                                                                                                                                          >  2:454704/1‑251 (MQ=255)
                                                                                           gcgAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATCTGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAAGAGGCGATTAAAGGCTGTCGTATCGGCATTAACACCGGTGCCGCGTAAAAGACCGGGCCGGACGGGCGCTCCATCAATAAATATATTg                                                                                                                                         >  2:659144/1‑251 (MQ=255)
                                                                                           gcgAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTg                                                                                                                                         <  1:695580/251‑1 (MQ=255)
                                                                                             taGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCCAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGACGGCGGCGAATAAAGGCTATCGTATCCGCCTTAAAACCGTTACCGCGAAACAAGTCGGGCCGGACGGTCGCTAACTCCATAAATGTAttcct                                                                                                                                       >  2:208561/2‑248 (MQ=255)
                                                                                             gAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCACACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGGTGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTCACGGCTATCGTATCGCCATTAAACCCGGTACCGCGAAAAAGGTAGGGCCGGCAGGTCGCTACATCAATAATTATATTGCt                                                                                                                                       >  2:531624/1‑251 (MQ=255)
                                                                                               gTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCtt                                                                                                                                      >  1:827022/1‑250 (MQ=255)
                                                                                                             cGGCAGCTACACCATCTATCGCACACGCTCGCTTACCAAAGTAGAACCCCCGGTTCCCCGTGACCGTGTATTCCCGGAACCCATTGTCCGCGCCGTGATGCATATGATGGAAAGCGTGGCGCTACCGAGCGGCGGCGGCGGGAAGGCGGCGATTAAAGGCTATCGTATCACCATACAAACCGGTACCGCGAAAAAGGTCGGGCCGGAAGGTAGCTACACCAATAACAATATTGCTTATACCGCAAGCGTTg                                                                                                                       >  2:119025/1‑251 (MQ=255)
                                                                                                                cAGCTACGGCCTTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgcgc                                                                                                                    >  1:164683/1‑251 (MQ=255)
                                                                                                                cAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgcgc                                                                                                                    >  1:578216/1‑251 (MQ=255)
                                                                                                                     tcGGCATTTATCGGCCCCTGTCGATTACCAAAGTGGTCCCCCCGGTTCCCGGGTAACGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTGCATATGATGGAAAGCTTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGATAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGc                                                                                                                <  2:783143/249‑1 (MQ=255)
                                                                                                                      cGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTTAACGTGTCTTCCCGGAATCCATTGGCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                              <  1:615461/251‑1 (MQ=255)
                                                                                                                         cATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGAATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTc                                                                                                           >  2:13906/1‑251 (MQ=255)
                                                                                                                              aTCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGTCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcg                                                                                                      <  2:235996/251‑1 (MQ=255)
                                                                                                                              aTCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGAGTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcg                                                                                                      <  1:545479/251‑1 (MQ=255)
                                                                                                                                          cGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGGCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCGACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCGTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCt                                                                                           <  2:698686/250‑1 (MQ=255)
                                                                                                                                           gATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTgg                                                                                         >  1:10902/1‑251 (MQ=255)
                                                                                                                                             ttACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGCAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgtt                                                                                       >  2:262588/1‑251 (MQ=255)
                                                                                                                                             ttACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGAGCTGgtt                                                                                       >  2:568413/1‑251 (MQ=255)
                                                                                                                                                ccAAAGGTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCAGTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                    <  2:793502/251‑1 (MQ=255)
                                                                                                                                                        gACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                            <  1:170692/251‑1 (MQ=255)
                                                                                                                                                        gACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCa                                                                             <  1:865073/250‑1 (MQ=255)
                                                                                                                                                           cccccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACg                                                                          >  1:616630/1‑250 (MQ=255)
                                                                                                                                                                    cccGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAgtc                                                                >  1:749850/1‑249 (MQ=255)
                                                                                                                                                                               tggcatCCCGGAATACATTGTACGCACCGTGGTGCATATGATGGAACGTGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGCATCGCCATTAACACCGGTACCGCGAAAAAGGTCGGGACGGCCGGTCGCTACTTAAATAAATATATTGCTTATACCGCCGGCGCTGCGCATGAGAGTCAGCCGCGCTTCGCGCGCGTTTATTATATCAACGTTCCGCTGGCGGGTacacact                                                     >  2:332791/6‑245 (MQ=255)
                                                                                                                                                                                        ggAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAATACTAaggcgggg                                            >  1:29612/1‑249 (MQ=255)
                                                                                                                                                                                           aTCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGCCGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAAAACTAcgggcgcaccc                                         >  2:690117/1‑243 (MQ=255)
                                                                                                                                                                                             cctttGTCCGCACCGTGTGTCATATGATGGAAAGCGTGGCGCTCCCGGGGGGCGGCGGCGTGAAGTCGGCGAGTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCTGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGt                                        <  2:778154/247‑1 (MQ=255)
                                                                                                                                                                                             ccATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAATACTACGGCGGCGCCggt                                       >  1:229783/1‑249 (MQ=255)
                                                                                                                                                                                                        accgtggtgCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGGAAATACAACAGCGGAGCCGttttcccgcccgg                            >  1:678070/4‑241 (MQ=255)
                                                                                                                                                                                                             ggtgCATATGAGGGAAAGCGTGGCGCTACCGGGCGGCGGCTGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGGCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCttt                        <  2:576681/250‑1 (MQ=255)
                                                                                                                                                                                                                        atgGAAAGCGTGGCGCTCCCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCcatcat             <  1:167928/251‑1 (MQ=255)
                                                                                                                                                                                                                        agggAAAGCGGTGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTAGCGTATCGCCATTAAAACCGGTACCGAGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCcatcat             <  2:84605/249‑1 (MQ=255)
                                                                                                                                                                                                                                cGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGAGCCGTTTCCCCGCCGGGCCTTTGTGGCCTCATGGggggcg     >  1:599203/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTTTCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGTGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTAt  <  2:616630/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
GCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTAT  >  NC_000913/92608‑93084

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