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

TTTCCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCT  >  NC_000913/3770294‑3770784
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
tttCCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACtt                                                                                                                                                                                                                                                  >  1:112385/1‑251 (MQ=255)
         cgccGATACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGTCGTAGATGCCGTCGGGCAAGGCATCGATATCATCGTTAGGGCCCAGTTCAATGGTTCCCCGCACACCGTCCGAGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAGCAGCCTGGCCCGACAGGTAAAATACTAGTCAGATTGACGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGcatcat                                                                                                                                                                                                                                         >  2:54061/3‑251 (MQ=255)
         cACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGGAGCGCGTTAAACATCACTTCCGaatcat                                                                                                                                                                                                                                         >  2:32177/1‑251 (MQ=255)
          aCCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCatcatc                                                                                                                                                                                                                                        >  1:201888/1‑251 (MQ=255)
             gAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGCTGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACCCCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCCCGACCGGTAAAATACTAGTCAGTTTGCCGTGAAACCCCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGcc                                                                                                                                                                                                                                     >  1:544789/1‑251 (MQ=255)
                  aTGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGtt                                                                                                                                                                                                                                >  1:210899/1‑251 (MQ=255)
                       tGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTc                                                                                                                                                                                                                              >  1:18902/1‑248 (MQ=255)
                         aaaTGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACg                                                                                                                                                                                                                         >  2:115681/1‑251 (MQ=255)
                                cGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAg                                                                                                                                                                                                                   <  1:275440/250‑1 (MQ=255)
                                 ggAGTAAACCGCCACCTGGGCGTCGATGCCGACGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTAGCCGTGAAACACCTGCCCAGGTAGCGCGTTACACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAg                                                                                                                                                                                                                   >  1:144451/1‑249 (MQ=255)
                                    gTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACCCCTGCCCAGGGAGCGCGTTAAACATCACTTCAGCATCATCGCCAGGTTTCAGACGTAACAGCGAGt                                                                                                                                                                                                              >  1:449923/1‑251 (MQ=255)
                                      aaaCCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGttt                                                                                                                                                                                                            >  2:38513/1‑251 (MQ=255)
                                      aaaCCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATAGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCTGACGTAACAGCGAGttt                                                                                                                                                                                                            >  1:372898/1‑251 (MQ=255)
                                                    gCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGc                                                                                                                                                                                              <  1:189754/251‑1 (MQ=255)
                                                       gagacgCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGAc                                                                                                                                                                                           <  1:304951/246‑1 (MQ=255)
                                                         gATGCCGTCGGGTAAGGCATCGATCTCATCGTTCGGGTCCAGTTCAATGGTTCCCAGCCCACCGTCCGTGCCGGGC‑ACGCCCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGACCCGCCTGGCACGACAGGTAAAATCCTCGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTACACATCACTTCCGCATCATCGCCAGGTTTAAGACGTAACAGCGAGTTTTGCCGAAATTGGTCGACa                                                                                                                                                                                          >  2:552935/1‑250 (MQ=255)
                                                              cGTCGGGGAAGGCATCGATATCATCTTTAGGGTCCAGTTCAATTGTTCCCAGCACACCGTCCGGGCCGGGC‑ACGACCGTTAATGATTGCAAGACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGGTGGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCGCTGCCGCATCATCGCCAGGTTTCAGGCGTAACAGCGAGTTTTGCCGAAATTGGTCGACAATTTGc                                                                                                                                                                                    <  2:219926/251‑1 (MQ=255)
                                                                  gggTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGtt                                                                                                                                                                                <  1:226942/251‑1 (MQ=255)
                                                                    gTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGtttt                                                                                                                                                                              >  1:259997/1‑251 (MQ=255)
                                                                     tAAGGCATCGATAACATCGTTAGGCTCCAGTTCAATGGTTCCCAGCACAACGTCCGTGCCGGGC‑ACGCCCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAACCATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGttttt                                                                                                                                                                             >  2:285740/1‑251 (MQ=255)
                                                                        ggCATCGATATCATCGTTAGGTTACAGTTCAATGGTTCCCCGCAACCCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACACCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGACGTGAAACACCTGCCCCGGTACCGCGTTAAACCTCACTTCCGCATCATCGCCAGGCTTCAGACGTAACAGAGCGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCt                                                                                                                                                                          >  2:549448/1‑251 (MQ=255)
                                                                             cGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCgggg                                                                                                                                                                     >  1:417306/1‑251 (MQ=255)
                                                                                tatCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAACCACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGTGCGACAATTTGCCGTTTTTGCTCGGGGATg                                                                                                                                                                  >  1:530541/1‑251 (MQ=255)
                                                                                       gTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTACGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGCTGAAGAcc                                                                                                                                                            >  1:22438/1‑250 (MQ=255)
                                                                                            ggTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCGGGGC‑ACGACCGTTAATGATTGCAATTCCCCCTGCGCCTGATAATAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACTCCTTCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCAGTTTTTGCTCGGGGATGAAGACCATCAcc                                                                                                                                                      <  1:482496/251‑1 (MQ=255)
                                                                                              tCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCgg                                                                                                                                                    >  1:334007/1‑251 (MQ=255)
                                                                                                 aGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTATACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACg                                                                                                                                                 <  1:192536/251‑1 (MQ=255)
                                                                                                    tCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGc                                                                                                                                                >  1:307793/1‑249 (MQ=255)
                                                                                                    tCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCa                                                                                                                                               <  2:417306/250‑1 (MQ=255)
                                                                                                      aaTGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAAt                                                                                                                                                                                                   <  2:116783/196‑1 (MQ=255)
                                                                                                      aaTGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAAt                                                                                                                                                                                                   >  1:116783/1‑196 (MQ=255)
                                                                                                      aaTGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCg                                                                                                                                            >  1:156920/1‑251 (MQ=255)
                                                                                                             cccAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                     >  1:262932/1‑251 (MQ=255)
                                                                                                               cAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGc                                                                                                                                     >  2:205855/1‑249 (MQ=255)
                                                                                                                agcccgaagtccgtGCCGGGCGAC‑ACCGTTAGTGATTGCAATACCCCCTACGCCTGATACGTACCGCCTGGCACGACCGGTAACATCCTACTCAGTTTGCCGTGACACACCTGCCCACGTAGAGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTAGGCGACAATTTGCCGTTTTTGCTCGGGGATGATGAGCATCACCGGACGCAGCGGCACGGCAGc                                                                                                                                  >  2:226953/9‑251 (MQ=255)
                                                                                                                       cgtccgtGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCCCGACAGGAAAAATACTAGTCAGTTTGCCGTGAACCACCTGCCCAGGTAGCGCGTTAAACCTCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGCGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGAACATCACCGGACGCAGCGGCACGGAAGCTGCGTAt                                                                                                                           >  1:372627/1‑251 (MQ=255)
                                                                                                                           cgtGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGCTAAGAACCGACTGGCACGACAGGAAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTCGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGCAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTAc                                                                                                                       >  2:325899/1‑251 (MQ=255)
                                                                                                                            gtGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTAcc                                                                                                                      >  1:401012/1‑251 (MQ=255)
                                                                                                                                cGGGC‑ACGACCTTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCTCATCATCGCCAGGGTTCAGACTTAACATCGAGTTTTGGCGAAATTGGGCTACAATTTGCCGTTTTTGCTCGTGGAGGAAGACCATCACCGGACGCAGCGGCAAGGCATCGGCGTATGTACCTggg                                                                                                                  <  2:503158/251‑1 (MQ=255)
                                                                                                                                  ggC‑ACGACCGTTAATGATTGCAATACCCCCGGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCg                                                                                                                <  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)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
TTTCCGCATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGC‑ACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCT  >  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: ‑