Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,917,424 A→G G389G (GGA→GGG yebS → intermembrane transport protein YebS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,917,4240AG99.7% 1250.9 / ‑4.3 304G389G (GGA→GGGyebSintermembrane transport protein YebS
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (153/150);  total (154/150)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

Only 100 of 308 total aligned reads displayed.

GTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAA  >  NC_000913/1917326‑1917521
                                                                                                  |                                                                                                 
gTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGcc                                                                                                 <  1:3210938/101‑1 (MQ=255)
gTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGcc                                                                                                 <  2:3211125/101‑1 (MQ=255)
 tGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCg                                                                                                >  2:3211468/1‑101 (MQ=255)
  gACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCgg                                                                                               <  1:3211005/101‑1 (MQ=255)
     cTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTg                                                                                            >  2:3211472/1‑101 (MQ=255)
     cTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTg                                                                                            <  1:3211078/101‑1 (MQ=255)
     cTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTg                                                                                            >  1:3211473/1‑101 (MQ=255)
      tGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGc                                                                                           <  2:3210838/101‑1 (MQ=255)
      tGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGc                                                                                           <  1:3210969/101‑1 (MQ=255)
          ttGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGttt                                                                                       <  2:3211177/101‑1 (MQ=255)
          ttGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGttt                                                                                       <  2:3210984/101‑1 (MQ=255)
           tGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGtttt                                                                                      >  2:3211482/1‑101 (MQ=255)
             gTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTAt                                                                                    >  1:3211484/1‑101 (MQ=255)
              tCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTAtt                                                                                   >  2:3211487/1‑101 (MQ=255)
                 cTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCg                                                                                <  1:3211032/101‑1 (MQ=255)
                 cTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCg                                                                                >  2:3211491/1‑101 (MQ=255)
                   ggTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGgc                                                                              <  2:3210875/101‑1 (MQ=255)
                    gTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGgcg                                                                             <  1:3211132/101‑1 (MQ=255)
                      cGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCa                                                                           >  1:3211493/1‑101 (MQ=255)
                        aTGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGc                                                                         >  1:3211497/1‑101 (MQ=255)
                        aTGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGc                                                                         <  2:3211109/101‑1 (MQ=255)
                          gCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGTGGCCGGCTGCGTTTTATTTCGGCGCAGCgg                                                                       >  1:3211501/1‑101 (MQ=255)
                          gCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCgg                                                                       <  2:3210711/101‑1 (MQ=255)
                           cTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                                      >  2:3211502/1‑101 (MQ=255)
                           cTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                                      >  2:3211506/1‑101 (MQ=255)
                           cTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                                      <  2:3210999/101‑1 (MQ=255)
                            tCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTa                                                                     <  1:3211038/101‑1 (MQ=255)
                             cGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTaa                                                                    >  1:3211510/1‑101 (MQ=255)
                              gACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAAt                                                                   >  2:3211511/1‑101 (MQ=255)
                              gACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAAt                                                                   <  2:3210917/101‑1 (MQ=255)
                               aCCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAAtt                                                                  >  1:3211514/1‑101 (MQ=255)
                                ccTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAAttt                                                                 <  1:3210857/101‑1 (MQ=255)
                                 cTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAAtttt                                                                <  2:3211214/101‑1 (MQ=255)
                                 cTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAAtttt                                                                <  1:3210985/101‑1 (MQ=255)
                                   gTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGa                                                              >  1:3211518/1‑101 (MQ=255)
                                     ttGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACt                                                            <  1:3211053/101‑1 (MQ=255)
                                      tGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTa                                                           <  1:3211016/101‑1 (MQ=255)
                                      tGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTa                                                           >  2:3211520/1‑101 (MQ=255)
                                       gTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTAt                                                          <  2:3211243/101‑1 (MQ=255)
                                       gTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTAt                                                          <  1:3210908/101‑1 (MQ=255)
                                         cATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTc                                                        >  1:3211525/1‑101 (MQ=255)
                                          atatCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCt                                                       >  2:3211531/1‑101 (MQ=255)
                                          atatCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCt                                                       >  1:3211529/1‑101 (MQ=255)
                                           tatCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCtt                                                      >  1:3211534/1‑101 (MQ=255)
                                            atCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGgc                                                                              >  1:3211538/1‑76 (MQ=255)
                                            atCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTg                                                     <  2:3211151/101‑1 (MQ=255)
                                            atCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTg                                                     <  2:3211099/101‑1 (MQ=255)
                                             tCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGc                                                    >  2:3211539/1‑101 (MQ=255)
                                                ttAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtgt                                                 <  2:3211036/101‑1 (MQ=255)
                                                  aaCCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTgg                                               >  1:3211544/1‑101 (MQ=255)
                                                     cATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAAt                                            >  2:3211546/1‑101 (MQ=255)
                                                     cATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAAt                                            <  2:3211358/101‑1 (MQ=255)
                                                      aTGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATg                                           <  2:3211183/101‑1 (MQ=255)
                                                       tGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtgt                                                 >  2:3211552/1‑94 (MQ=255)
                                                       tGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATgg                                          <  2:3211077/101‑1 (MQ=255)
                                                       tGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATgg                                          >  1:3211551/1‑101 (MQ=255)
                                                        gTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGc                                         <  2:3211141/101‑1 (MQ=255)
                                                         tCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCt                                        <  1:3211090/101‑1 (MQ=255)
                                                          cGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTg                                       >  1:3211557/1‑101 (MQ=255)
                                                          cGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTg                                       <  1:3210930/101‑1 (MQ=255)
                                                           gCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTa                                                                     >  1:3211561/1‑70 (MQ=255)
                                                           gCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTgg                                      >  2:3211558/1‑101 (MQ=255)
                                                              gATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGTTGGACa                                   <  1:3211227/101‑1 (MQ=255)
                                                              gATTAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACa                                   >  1:3211565/1‑101 (MQ=255)
                                                                ttAATCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGc                                 >  2:3211567/1‑101 (MQ=255)
                                                                   aTCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGc                              >  1:3211571/1‑101 (MQ=255)
                                                                   aTCGCGATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGc                              >  1:3211574/1‑101 (MQ=255)
                                                                     cgcgATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCtt                            >  2:3211578/1‑101 (MQ=255)
                                                                     cgcgATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCtt                            >  1:3211579/1‑101 (MQ=255)
                                                                     cgcgATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCtt                            <  1:3211116/101‑1 (MQ=255)
                                                                      gcgATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCg                                                                                >  2:3211583/1‑48 (MQ=255)
                                                                       cgATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTAc                          <  1:3210737/101‑1 (MQ=255)
                                                                        gATCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACt                         >  1:3211585/1‑101 (MQ=255)
                                                                         aTCAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACtt                        >  1:3211588/1‑101 (MQ=255)
                                                                           cAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTGTTCTTGCTGTGGAATGGCTGGACAGCCGCTTACtt                        <  1:3211182/99‑1 (MQ=255)
                                                                           cAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACttt                       <  2:3210965/100‑1 (MQ=255)
                                                                           cAGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTg                      <  2:3210957/101‑1 (MQ=255)
                                                                            aGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTgg                     >  2:3211591/1‑101 (MQ=255)
                                                                            aGATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTgg                     <  1:3211013/101‑1 (MQ=255)
                                                                             gATCCTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTggg                    <  1:3211300/101‑1 (MQ=255)
                                                                                ccTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATg                 >  1:3211596/1‑101 (MQ=255)
                                                                                 cTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGc                <  2:3211270/101‑1 (MQ=255)
                                                                                 cTCGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGc                <  2:3211192/101‑1 (MQ=255)
                                                                                   cGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGcac              >  2:3211599/1‑101 (MQ=255)
                                                                                   cGCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGcac              <  1:3211026/101‑1 (MQ=255)
                                                                                   cGCTTTTACTAGGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCt                             <  1:3211061/86‑1 (MQ=255)
                                                                                    gCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGcaca             <  1:3211166/101‑1 (MQ=255)
                                                                                    gCTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGcaca             >  1:3211600/1‑101 (MQ=255)
                                                                                     cTTTTACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACAt            >  1:3211605/1‑101 (MQ=255)
                                                                                      ttttACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATg           >  1:3211610/1‑101 (MQ=255)
                                                                                      ttttACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATg           >  2:3211609/1‑101 (MQ=255)
                                                                                       tttACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGa          <  1:3211142/101‑1 (MQ=255)
                                                                                         tACTATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGt        >  1:3211614/1‑101 (MQ=255)
                                                                                            tATGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAg     >  1:3211617/1‑101 (MQ=255)
                                                                                             aTGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAgg    <  2:3211048/101‑1 (MQ=255)
                                                                                             aTGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAgg    >  1:3211621/1‑101 (MQ=255)
                                                                                             aTGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAgg    >  2:3211622/1‑101 (MQ=255)
                                                                                              tGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGa   <  2:3211307/101‑1 (MQ=255)
                                                                                              tGGGGCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGa   <  1:3211478/101‑1 (MQ=255)
                                                                                               ggggCCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGaa  <  1:3211110/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAA  >  NC_000913/1917326‑1917521

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