Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,831,329 A→G E291E (GAA→GAG yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,831,3290AG99.0% 380.8 / NA 102E291E (GAA→GAGyicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (48/53);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 104 total aligned reads displayed.

TTGCGCCATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAGGG  >  NC_000913/3831181‑3831476
                                                                                                                                                    |                                                                                                                                                   
ttGCGCCATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAggg                                                                                                                                                   >  2:1194370/1‑151 (MQ=255)
  gcgcCATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGc                                                                                                                                                  <  1:1125330/150‑1 (MQ=255)
   cgcCATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCa                                                                                                                                                 >  2:336358/1‑150 (MQ=255)
   cgcCATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCa                                                                                                                                                 >  1:128027/1‑150 (MQ=255)
      cATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACt                                                                                                                                             >  2:859507/1‑151 (MQ=255)
         tACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTc                                                                                                                                          >  1:404353/1‑151 (MQ=255)
          aCGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTc                                                                                                                                          <  4:433158/150‑1 (MQ=255)
          aCGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTc                                                                                                                                          <  3:437464/150‑1 (MQ=255)
          aCGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCg                                                                                                                                         >  1:406003/1‑151 (MQ=255)
           cGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCg                                                                                                                                         >  1:159531/1‑150 (MQ=255)
           cGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGa                                                                                                                                        <  2:485617/151‑1 (MQ=255)
           cGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGa                                                                                                                                        <  1:1739667/151‑1 (MQ=255)
                 gTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGaa                                                                                                                                    <  2:176255/149‑1 (MQ=255)
                 gTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAAc                                                                                                                                   >  1:398233/1‑150 (MQ=255)
                 gTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACg                                                                                                                                  >  2:287810/1‑151 (MQ=255)
                 gTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACg                                                                                                                                  >  3:568911/1‑151 (MQ=255)
                  tCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACg                                                                                                                                  >  1:440913/1‑150 (MQ=255)
                   cccGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCt                                                                                                                                >  4:11049/1‑151 (MQ=255)
                    ccGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTa                                                                                                                               <  2:1655421/151‑1 (MQ=255)
                    ccGTCGTTGCAGATTGGTCACCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTa                                                                                                                               <  1:31294/151‑1 (MQ=255)
                     cgtcgtTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCt                                                                                                                                >  1:809635/1‑149 (MQ=255)
                     cgtcgtTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAg                                                                                                                              <  2:666931/151‑1 (MQ=255)
                      gtcgtTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGc                                                                                                                             >  3:175270/1‑151 (MQ=255)
                      gtcgtTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGc                                                                                                                             <  1:1565319/151‑1 (MQ=255)
                        cgtTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCg                                                                                                                            >  1:148186/1‑150 (MQ=255)
                         gtTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAg                                                                                                                          <  3:349007/151‑1 (MQ=255)
                          ttGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGt                                                                                                                         >  4:409104/1‑151 (MQ=255)
                           tGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGtt                                                                                                                        >  1:979492/1‑151 (MQ=255)
                            gCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTc                                                                                                                       <  3:255759/151‑1 (MQ=255)
                             cAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTc                                                                                                                       >  2:1665624/1‑150 (MQ=255)
                                aTTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCAttt                                                                                                                   >  3:320054/1‑151 (MQ=255)
                                  tGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTAt                                                                                                                 >  2:901915/1‑151 (MQ=255)
                                   ggTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATg                                                                                                                <  2:327696/151‑1 (MQ=255)
                                      cGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGtt                                                                                                             <  2:1371659/151‑1 (MQ=255)
                                       gCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTc                                                                                                            >  2:1207719/1‑151 (MQ=255)
                                        ccTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCg                                                                                                           >  3:470455/1‑151 (MQ=255)
                                         cTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGAGAGTTCATTTATGGTTCGc                                                                                                          >  1:1084337/1‑151 (MQ=255)
                                          tGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCt                                                                                                         >  1:1549035/1‑151 (MQ=255)
                                          tGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCt                                                                                                         >  2:795509/1‑151 (MQ=255)
                                           ggAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGAAGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTg                                                                                                        >  1:404133/1‑151 (MQ=255)
                                             aaGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGc                                                                                                       >  1:1301925/1‑150 (MQ=255)
                                             aaGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGc                                                                                                       >  2:581822/1‑150 (MQ=255)
                                              aGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCat                                                                                                     <  1:685948/151‑1 (MQ=255)
                                              aGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCat                                                                                                     <  1:1625213/151‑1 (MQ=255)
                                                 gATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCattt                                                                                                   >  3:486228/1‑150 (MQ=255)
                                                   tCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCatttatt                                                                                                >  2:700392/1‑151 (MQ=255)
                                                    cGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCatttattt                                                                                               >  1:1699237/1‑151 (MQ=255)
                                                           cGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCg                                                                                         >  2:1623901/1‑150 (MQ=255)
                                                           cGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGa                                                                                        <  2:1169490/151‑1 (MQ=255)
                                                             tttGCAAGGTCCGGACTGGGCGGTGACCGGCCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGAccc                                                                                          <  4:433381/147‑1 (MQ=255)
                                                              ttGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGGTGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGAtt                                                                                      >  2:181173/1‑150 (MQ=255)
                                                               tGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATtat                                                                                    >  1:1075533/1‑151 (MQ=255)
                                                                gCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATtat                                                                                    <  4:372939/150‑1 (MQ=255)
                                                                gCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATtat                                                                                    <  2:1718967/150‑1 (MQ=255)
                                                                 cAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATaa                                                                                  <  4:140951/151‑1 (MQ=255)
                                                                 cAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATaa                                                                                  <  1:1454770/151‑1 (MQ=255)
                                                                  aaGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATaa                                                                                  >  1:416251/1‑150 (MQ=255)
                                                                    ggTCCGGACTGGGCGGTGACCGACCCCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAaca                                                                               <  2:1293238/151‑1 (MQ=255)
                                                                     gTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAacac                                                                              >  2:1238926/1‑151 (MQ=255)
                                                                       ccGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACgg                                                                            >  1:282710/1‑151 (MQ=255)
                                                                         ggACTGGGCGGTGACCGTCCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGa                                                                           <  1:431391/150‑1 (MQ=255)
                                                                         ggACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGaa                                                                          <  4:320551/151‑1 (MQ=255)
                                                                          gACTGGGCGGTGAACGGGCTCAATATCAGCTTGCGCAACATGACCCTCAGGAAAGATGGCTGGCAGACGAAGGAGGGCAAATTGTCGATGAACGCTAGCGAGTTCAGTTATGGTTCGCTGCCTTTAATTGGCCCGATTATaaccttggaa                                                                          <  1:1049983/150‑9 (MQ=255)
                                                                              gggCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACGGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTc                                                                     >  1:1514065/1‑151 (MQ=255)
                                                                                gCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTccc                                                                   <  2:257359/151‑1 (MQ=255)
                                                                                gCGGAGACCGACCTCGATCTCATCTTGCGCAACATGACTTTCAGTAAAGATGACTGGCAGACACAGGAAGGCAAACAGTCGTTGAACGCTAGCGAGCTCATTTATGGTTCTCTGCATTTATTTGACCCGATTATAAACACGGAATTTTccc                                                                   >  1:239024/1‑151 (MQ=255)
                                                                                  ggTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCg                                                                 >  4:450853/1‑151 (MQ=255)
                                                                                   gTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGc                                                                >  2:287663/1‑151 (MQ=255)
                                                                                      accgaccTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAgg                                                             >  1:757813/1‑151 (MQ=255)
                                                                                       ccgaccTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAggg                                                            <  1:1348006/151‑1 (MQ=255)
                                                                                       ccgaccTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAgg                                                             <  1:837637/150‑1 (MQ=255)
                                                                                       ccgaccTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTAGTTGACCCGATTATAAACACGGAATTTTCCCCGCAgg                                                             <  3:136444/150‑1 (MQ=255)
                                                                                         gaccTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCg                                                          <  2:1621163/151‑1 (MQ=255)
                                                                                         gaccTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCg                                                          <  2:610720/151‑1 (MQ=255)
                                                                                                    cAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAgcg                                                     <  2:516951/145‑1 (MQ=255)
                                                                                                         tGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGtt                                          >  1:113618/1‑151 (MQ=255)
                                                                                                         tGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGt                                           <  2:758243/150‑1 (MQ=255)
                                                                                                          gcgcAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGtt                                          <  3:409124/150‑1 (MQ=255)
                                                                                                          gcgcAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGtt                                          >  1:1335161/1‑150 (MQ=255)
                                                                                                             cAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCAcc                                      >  1:1113931/1‑151 (MQ=255)
                                                                                                                     ccTTCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCt                               <  2:1058556/150‑1 (MQ=255)
                                                                                                                       ttCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTggg                            >  2:1697446/1‑151 (MQ=255)
                                                                                                                        tCAGTAAATATGCCTGGCACACACAGAAGGCCACAGTATCGATGAACGATAGCGGGTGCATTTATGGTTCGCTGCATTTATTTGACCCGTTTATAAACACGGAATTTTCCGCGCAGAGCATAACGCTACGCCAGTTAACCAGACGCgggg                            <  1:655637/150‑1 (MQ=255)
                                                                                                                        tCAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTggg                            <  1:287769/150‑1 (MQ=255)
                                                                                                                         cAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGaa                          <  3:200672/151‑1 (MQ=255)
                                                                                                                         cAGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGa                           <  1:1030104/150‑1 (MQ=255)
                                                                                                                          aGTAAAGATGGCTGGCGGGCACAGGAGGGCAAACTGTCGATGAGCGCTAGCGCGTACTTTTATGCTTCGCTGGTTTTATTTGACCCGATTATAAACACGGCATTTTCCCCGCAGGGCGTCGCGCTGCGCCAGTTCACCAGAGGCTGGGaa                          <  1:126744/150‑1 (MQ=255)
                                                                                                                          aGTAAAGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCCTTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACCCGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGaa                          <  1:183743/150‑1 (MQ=255)
                                                                                                                             aaaGATGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAgggg                      >  3:266764/1‑151 (MQ=255)
                                                                                                                                 aTGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTAt                  <  1:859395/151‑1 (MQ=255)
                                                                                                                                 aTGACTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTAt                  <  3:349532/151‑1 (MQ=255)
                                                                                                                                     cTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGt               <  1:1088428/150‑1 (MQ=255)
                                                                                                                                     cTGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATGTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGt               <  4:140232/150‑1 (MQ=255)
                                                                                                                                      tGGCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTc              >  1:1595252/1‑150 (MQ=255)
                                                                                                                                       ggCAGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCa             <  3:379505/150‑1 (MQ=255)
                                                                                                                                          aGACACAGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAAc         <  1:520572/151‑1 (MQ=255)
                                                                                                                                             cacaGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTc      >  2:1071275/1‑151 (MQ=255)
                                                                                                                                               caGGAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGATGGGAAGGGGGTATGGTCAGAACGTCa     <  2:406064/150‑1 (MQ=255)
                                                                                                                                                aggaggGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAgg   >  3:315828/1‑151 (MQ=255)
                                                                                                                                                 ggAGGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAggg  <  2:1084469/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGCGCCATTACGCTCTGTCCCGTCGTTGCAGATTGGTCGCCTGGAAGTGATCGATGCTCGTTTGCAAGGTCCGGACTGGGCGGTGACCGACCTCGATCTCAGCTTGCGCAACATGACCTTCAGTAAAGATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAGGG  >  NC_000913/3831181‑3831476

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