Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,770,539 C→T V281M (GTG→ATG)  yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,5390CT97.3% 250.7 / NA 74V281M (GTG→ATG) yibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (45/27);  minor base G (0/2);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TTTCCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCT  >  NC_000913/3770294‑3770784
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
tttCCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACtt                                                                                                                                                                                                                                                  >  1:112385/1‑251 (MQ=255)
         cgccGATACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGTCGTAGATGCCGTCGGGCAAGGCATCGATATCATCGTTAGGGCCCAGTTCAATGGTTCCCCGCACACCGTCCGAGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAGCAGCCTGGCCCGACAGGTAAAATACTAGTCAGATTGACGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGcatcat                                                                                                                                                                                                                                         >  2:54061/3‑251 (MQ=255)
         cACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGGAGCGCGTTAAACATCACTTCCGaatcat                                                                                                                                                                                                                                         >  2:32177/1‑251 (MQ=255)
          aCCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCatcatc                                                                                                                                                                                                                                        >  1:201888/1‑251 (MQ=255)
             gAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGCTGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACCCCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCCCGACCGGTAAAATACTAGTCAGTTTGCCGTGAAACCCCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGcc                                                                                                                                                                                                                                     >  1:544789/1‑251 (MQ=255)
                  aTGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGtt                                                                                                                                                                                                                                >  1:210899/1‑251 (MQ=255)
                       tGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTc                                                                                                                                                                                                                              >  1:18902/1‑248 (MQ=255)
                         aaaTGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACg                                                                                                                                                                                                                         >  2:115681/1‑251 (MQ=255)
                                cGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAg                                                                                                                                                                                                                   <  1:275440/250‑1 (MQ=255)
                                 ggAGTAAACCGCCACCTGGGCGTCGATGCCGACGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTAGCCGTGAAACACCTGCCCAGGTAGCGCGTTACACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAg                                                                                                                                                                                                                   >  1:144451/1‑249 (MQ=255)
                                    gTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACCCCTGCCCAGGGAGCGCGTTAAACATCACTTCAGCATCATCGCCAGGTTTCAGACGTAACAGCGAGt                                                                                                                                                                                                              >  1:449923/1‑251 (MQ=255)
                                      aaaCCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGttt                                                                                                                                                                                                            >  2:38513/1‑251 (MQ=255)
                                      aaaCCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATAGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCTGACGTAACAGCGAGttt                                                                                                                                                                                                            >  1:372898/1‑251 (MQ=255)
                                                    gCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGc                                                                                                                                                                                              <  1:189754/251‑1 (MQ=255)
                                                       gagacgCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGAc                                                                                                                                                                                           <  1:304951/246‑1 (MQ=255)
                                                         gATGCCGTCGGGTAAGGCATCGATCTCATCGTTCGGGTCCAGTTCAATGGTTCCCAGCCCACCGTCCGTGCCGGGCACGCCCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGACCCGCCTGGCACGACAGGTAAAATCCTCGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTACACATCACTTCCGCATCATCGCCAGGTTTAAGACGTAACAGCGAGTTTTGCCGAAATTGGTCGACa                                                                                                                                                                                          >  2:552935/1‑250 (MQ=255)
                                                              cGTCGGGGAAGGCATCGATATCATCTTTAGGGTCCAGTTCAATTGTTCCCAGCACACCGTCCGGGCCGGGCACGACCGTTAATGATTGCAAGACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGGTGGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCGCTGCCGCATCATCGCCAGGTTTCAGGCGTAACAGCGAGTTTTGCCGAAATTGGTCGACAATTTGc                                                                                                                                                                                    <  2:219926/251‑1 (MQ=255)
                                                                  gggTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGtt                                                                                                                                                                                <  1:226942/251‑1 (MQ=255)
                                                                    gTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGtttt                                                                                                                                                                              >  1:259997/1‑251 (MQ=255)
                                                                     tAAGGCATCGATAACATCGTTAGGCTCCAGTTCAATGGTTCCCAGCACAACGTCCGTGCCGGGCACGCCCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAACCATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGttttt                                                                                                                                                                             >  2:285740/1‑251 (MQ=255)
                                                                        ggCATCGATATCATCGTTAGGTTACAGTTCAATGGTTCCCCGCAACCCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACACCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGACGTGAAACACCTGCCCCGGTACCGCGTTAAACCTCACTTCCGCATCATCGCCAGGCTTCAGACGTAACAGAGCGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCt                                                                                                                                                                          >  2:549448/1‑251 (MQ=255)
                                                                             cGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCgggg                                                                                                                                                                     >  1:417306/1‑251 (MQ=255)
                                                                                tatCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAACCACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGTGCGACAATTTGCCGTTTTTGCTCGGGGATg                                                                                                                                                                  >  1:530541/1‑251 (MQ=255)
                                                                                       gTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTACGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGCTGAAGAcc                                                                                                                                                            >  1:22438/1‑250 (MQ=255)
                                                                                            ggTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCGGGGCACGACCGTTAATGATTGCAATTCCCCCTGCGCCTGATAATAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACTCCTTCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCAGTTTTTGCTCGGGGATGAAGACCATCAcc                                                                                                                                                      <  1:482496/251‑1 (MQ=255)
                                                                                              tCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCgg                                                                                                                                                    >  1:334007/1‑251 (MQ=255)
                                                                                                 aGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTATACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACg                                                                                                                                                 <  1:192536/251‑1 (MQ=255)
                                                                                                    tCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGc                                                                                                                                                >  1:307793/1‑249 (MQ=255)
                                                                                                    tCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCa                                                                                                                                               <  2:417306/250‑1 (MQ=255)
                                                                                                      aaTGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAAt                                                                                                                                                                                                   <  2:116783/196‑1 (MQ=255)
                                                                                                      aaTGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAAt                                                                                                                                                                                                   >  1:116783/1‑196 (MQ=255)
                                                                                                      aaTGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCg                                                                                                                                            >  1:156920/1‑251 (MQ=255)
                                                                                                             cccAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                     >  1:262932/1‑251 (MQ=255)
                                                                                                               cAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                     >  2:205855/1‑249 (MQ=255)
                                                                                                                agcccgaagtccgtGCCGGGCGACACCGTTAGTGATTGCAATACCCCCTACGCCTGATACGTACCGCCTGGCACGACCGGTAACATCCTACTCAGTTTGCCGTGACACACCTGCCCACGTAGAGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTAGGCGACAATTTGCCGTTTTTGCTCGGGGATGATGAGCATCACCGGACGCAGCGGCACGGCAGc                                                                                                                                  >  2:226953/9‑251 (MQ=255)
                                                                                                                       cgtccgtGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCCCGACAGGAAAAATACTAGTCAGTTTGCCGTGAACCACCTGCCCAGGTAGCGCGTTAAACCTCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGCGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGAACATCACCGGACGCAGCGGCACGGAAGCTGCGTAt                                                                                                                           >  1:372627/1‑251 (MQ=255)
                                                                                                                           cgtGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGCTAAGAACCGACTGGCACGACAGGAAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTCGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGCAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTAc                                                                                                                       >  2:325899/1‑251 (MQ=255)
                                                                                                                            gtGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTAcc                                                                                                                      >  1:401012/1‑251 (MQ=255)
                                                                                                                                cGGGCACGACCTTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCTCATCATCGCCAGGGTTCAGACTTAACATCGAGTTTTGGCGAAATTGGGCTACAATTTGCCGTTTTTGCTCGTGGAGGAAGACCATCACCGGACGCAGCGGCAAGGCATCGGCGTATGTACCTggg                                                                                                                  <  2:503158/251‑1 (MQ=255)
                                                                                                                                  ggCACGACCGTTAATGATTGCAATACCCCCGGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCg                                                                                                                <  2:340750/251‑1 (MQ=255)
                                                                                                                                       cGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATc                                                                                                            >  2:32553/1‑251 (MQ=255)
                                                                                                                                        gACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCa                                                                                                           <  1:5257/251‑1 (MQ=255)
                                                                                                                                         aCCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACAGCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTGGCCGAGATTGGGCGACAATTTGCCGTTTTGGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGGATGTACCTGGGCGGATCAg                                                                                                          <  2:544789/251‑1 (MQ=255)
                                                                                                                                                  gATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACTCCTGCCCAGGTAGCGCGTTAAACAGCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATGTCCGGACGCAGCGGCAAGGCAGCTGCGTATGTTCCTGGGCGGATCAGTACCTGAGt                                                                                                 <  2:93811/251‑1 (MQ=255)
                                                                                                                                                        aaTACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGt                                                                                            >  2:242553/1‑250 (MQ=255)
                                                                                                                                                           ccccccTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGcc                                                                                        <  1:38513/250‑1 (MQ=255)
                                                                                                                                                           aCCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGcc                                                                                        <  1:200062/251‑1 (MQ=255)
                                                                                                                                                                  gcgcCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTc                                                                                 >  1:130106/1‑251 (MQ=255)
                                                                                                                                                                   cgcCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCg                                                                                <  1:463644/251‑1 (MQ=255)
                                                                                                                                                                            aGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGa                                                                        >  1:4599/1‑250 (MQ=255)
                                                                                                                                                                             gAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGCTCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                      >  1:421914/1‑251 (MQ=255)
                                                                                                                                                                             gAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                      <  2:196563/251‑1 (MQ=255)
                                                                                                                                                                             gAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                      <  2:439378/251‑1 (MQ=255)
                                                                                                                                                                             aaaCCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAAt                                                                      <  1:327986/250‑1 (MQ=255)
                                                                                                                                                                              aaacGCCTGGCACGACAGGTAAAATACTAGGCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGATTTTTGCCGAAATTGGGCGACAATTGGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATg                                                                     <  2:67379/248‑1 (MQ=255)
                                                                                                                                                                               aCCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACCTAGCCATTGCTCGGCGCGCGAATGa                                                                    <  2:431486/251‑1 (MQ=255)
                                                                                                                                                                                         aCGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTGCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGAGGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGGAACGTAGCCATTGCTTGGCGCGCGAATGACAGTCTGCTc                                                          <  2:394967/251‑1 (MQ=255)
                                                                                                                                                                                            aCAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCAGCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCAGTGCTCGGCGCGCGAATGACAGTCTGCTCAAg                                                       <  2:359909/251‑1 (MQ=255)
                                                                                                                                                                                              aGGTAAAATACTAGTCAGGTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGt                                                     <  2:210899/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaaTACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTtat                                                  >  2:133950/1‑250 (MQ=255)
                                                                                                                                                                                                       aCTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTg                                            >  1:400354/1‑251 (MQ=255)
                                                                                                                                                                                                          aGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCACGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTgatt                                         >  2:404803/1‑248 (MQ=255)
                                                                                                                                                                                                           gTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACCGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTAAGGGTTATATTTTGCTTc                                        >  2:518972/1‑251 (MQ=255)
                                                                                                                                                                                                             cAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTc                                        >  1:419225/1‑249 (MQ=255)
                                                                                                                                                                                                              aGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGt                                     >  1:548910/1‑251 (MQ=255)
                                                                                                                                                                                                              aGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGt                                     >  1:521607/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCgg                                                                               >  2:93728/1‑206 (MQ=255)
                                                                                                                                                                                                                 ttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCgg                                                                               <  1:93728/206‑1 (MQ=255)
                                                                                                                                                                                                                 ttGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTa                                    >  1:458191/1‑249 (MQ=255)
                                                                                                                                                                                                                             acCTGCCCAGGTAGAGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGCGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTa                       >  1:79073/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGGCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCt                   <  1:552935/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATGTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCt           <  2:201888/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cgcgGTAAACATCACTTCCGCAGCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCg        <  1:406195/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaaCATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGCGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCCGTACCTGAGTAACGTCGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCt  >  1:552788/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
TTTCCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCT  >  NC_000913/3770294‑3770784

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