Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,171 T→A 19.4% intergenic (+43/‑187) tufB → / → secE translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1710TA19.4% 340.4 / 76.2 154intergenic (+43/‑187)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (60/64);  new base A (14/16);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

Only 100 of 167 total aligned reads displayed.

TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTAT  >  NC_000913/4177073‑4177268
                                                                                                  |                                                                                                 
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                 <  2:1760580/100‑1 (MQ=255)
   gTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                               <  2:2064341/100‑1 (MQ=255)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                              <  1:3438766/100‑1 (MQ=255)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                              <  1:3938030/100‑1 (MQ=255)
     gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                             >  1:1321411/1‑100 (MQ=255)
         gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                         <  2:4131091/100‑1 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                        <  2:2857075/100‑1 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttt                                                                                       <  1:1021474/100‑1 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttt                                                                                       >  1:1392336/1‑100 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                       <  1:3157861/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCTGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                       <  2:471822/99‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtttt                                                                                      <  2:1878428/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtttt                                                                                      >  2:1280021/1‑100 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                      <  1:145082/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                    <  1:608747/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                    <  1:3731095/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                    <  2:3043444/100‑3 (MQ=255)
               gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGc                                                                                   <  1:4145918/100‑1 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                  <  1:3476300/100‑5 (MQ=38)
                 aCCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCAc                                                                                 <  1:2438055/100‑1 (MQ=255)
                 aCCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtgc                                                                                 <  1:995596/100‑6 (MQ=38)
                  ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                >  2:2897653/1‑100 (MQ=255)
                      tGGCGCGGGCGTTGTAGCAAAAGATCTGAGCTAATTGCCGATAACATATGACGAAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCttt                                                                            >  2:3395673/1‑100 (MQ=255)
                        gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                          <  2:2388130/100‑1 (MQ=255)
                         cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCAATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                         <  1:2984209/100‑1 (MQ=255)
                         cgcgGGCGGTGTAGCAAAAGTTCAGAACTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCACTTTTTGCACGCTTTCGt                                                                         >  2:1421881/1‑100 (MQ=255)
                          gcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTa                                                                        <  1:1245370/100‑1 (MQ=255)
                           cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                       <  1:3099029/100‑1 (MQ=255)
                            gggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                      <  1:3769233/100‑1 (MQ=255)
                            gggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                      <  1:3110035/100‑1 (MQ=255)
                             ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                     <  1:1637875/100‑1 (MQ=255)
                             ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                     <  1:3996367/100‑1 (MQ=255)
                              gCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAg                                                                    <  1:3337102/100‑1 (MQ=255)
                                   gTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCt                                                               >  2:4162399/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            >  2:2186962/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            <  1:137841/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            <  2:1226781/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            <  1:1686120/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            <  1:1986386/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                            >  1:2672344/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                           >  2:963991/1‑100 (MQ=255)
                                             ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                     <  1:3200847/100‑1 (MQ=255)
                                                tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                  >  1:2574718/1‑100 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                 >  1:47281/1‑100 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                 >  1:1550308/1‑100 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGAACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCGTTCGTACCAGAACCTGGCTCATCAGTGAt                                                 >  2:2523401/1‑100 (MQ=255)
                                                     gatctGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                             <  2:1851083/96‑1 (MQ=255)
                                                      aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                            >  1:1834480/1‑100 (MQ=255)
                                                       atcgCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                           >  2:3931077/4‑100 (MQ=255)
                                                       aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                           <  2:1554147/100‑1 (MQ=255)
                                                        ttGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                          <  2:3827775/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                         <  2:1594067/100‑1 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                        >  1:4118981/1‑100 (MQ=255)
                                                           ccGAGAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                       <  2:1257771/100‑1 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                 >  1:3595923/1‑100 (MQ=255)
                                                                  cATTTGACGCAATGCGCACTAAACGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCa                                >  2:359797/1‑100 (MQ=255)
                                                                  cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCCTCAGTGATTTTCTTTGTCATAATCa                                >  2:4172310/1‑100 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                               >  2:2352356/1‑100 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                               <  1:647547/100‑1 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                               >  1:629235/1‑100 (MQ=255)
                                                                     ttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATACTCAtt                              >  2:1992865/1‑99 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                            <  2:1130849/100‑1 (MQ=255)
                                                                      tGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                            <  2:1130824/100‑1 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                           <  1:230755/100‑1 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                         <  2:160202/100‑1 (MQ=255)
                                                                          gCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                        >  2:140038/1‑100 (MQ=255)
                                                                          gCAATGCGCACGAAAAGGGCATCATCTTGTGTCTTTTTTGTACGCTTGCGTACCAGAACCTGGCGCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                        <  2:3325813/100‑1 (MQ=255)
                                                                           cAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgaga                       >  1:1591035/1‑100 (MQ=255)
                                                                              tGCTCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                    <  2:1382464/100‑1 (MQ=255)
                                                                               gcgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                   <  1:3384284/100‑1 (MQ=255)
                                                                                tgtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGACCCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                  <  1:4166572/94‑1 (MQ=38)
                                                                                 gtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                 >  1:3362624/6‑100 (MQ=38)
                                                                                 gtgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                 >  1:2020511/6‑100 (MQ=38)
                                                                                  cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                >  1:2607891/1‑100 (MQ=255)
                                                                                    caaaagGGGCATCAAATGATGCCCGTTTTGCACGATTTTGTACCAGAACCTGGCTCATCAGTGATTTGCGTTGTCATAATCATTGCTGAGACAGGctct               >  2:3561168/3‑99 (MQ=255)
                                                                                    cTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg              >  2:4189848/1‑100 (MQ=255)
                                                                                     aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt             >  1:2649242/2‑100 (MQ=255)
                                                                                      aaaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt            <  1:951468/100‑1 (MQ=255)
                                                                                      aaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt            >  2:3662494/1‑100 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTTTTTTTTTTTGTACCAGAACATGGCTCATAAGTGATTTTTTTTGTTATAATTCTTGCTGAGTCTGTCTttgttg           >  2:985504/1‑100 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg           >  1:1758069/1‑100 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg           >  1:3571828/1‑100 (MQ=255)
                                                                                       aaaGGGCATCAAATGATGCCCTTTTTTTTTTTTTTTTTACCAGAACATGGCTCATCAGTGATTTTTTTTGTTATAATCATTGCTGAGACAGGCTCTGTTg           >  2:2010247/1‑100 (MQ=255)
                                                                                        aaGGGCATCATTTGATGCCCTTTTTTTTTGTTTTCGTACCAGAACCTGGCTCATCAGTTATTTTCTTTGTCATAATCATTGCTGAGTCAGGCTCTGtt            >  2:642125/1‑98 (MQ=255)
                                                                                        aaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa          <  2:3892327/100‑1 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTTTACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACCGGCTCTGTTGa          >  2:3325079/1‑100 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa          >  1:342654/1‑100 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTAACAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa          >  2:2582514/1‑100 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACCGGCTCTGTTGa          >  1:3758256/1‑100 (MQ=37)
                                                                                          gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg        >  1:2402061/1‑100 (MQ=255)
                                                                                          gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg        <  1:1055727/100‑1 (MQ=255)
                                                                                           ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg       <  2:1074176/100‑1 (MQ=255)
                                                                                           ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg       >  1:1128890/1‑100 (MQ=255)
                                                                                            gCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc      >  1:385075/1‑100 (MQ=255)
                                                                                            gCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc      <  2:530156/100‑1 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg     <  2:2782621/100‑1 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg     <  2:474763/100‑1 (MQ=255)
                                                                                               tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTTTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCgga   <  1:2834003/100‑3 (MQ=255)
                                                                                                catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtat  <  2:2456976/100‑1 (MQ=255)
                                                                                                catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtat  <  1:1213490/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTAT  >  NC_000913/4177073‑4177268

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