Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,172 T→A 15.8% intergenic (+44/‑186) 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,1720TA15.8% 256.6 / 41.3 108intergenic (+44/‑186)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (41/50);  new base A (8/9);  total (49/59)
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.16e-01

Only 100 of 113 total aligned reads displayed.

TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAAT  >  NC_000913/4177073‑4177271
                                                                                                   |                                                                                                   
tCCGTGAAGGCGGCCGTACCGATGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCCATAACATTTGACGCAATGCGCACTAAAAGGGGATCAtt                                                                                                     >  2:1891937/1‑100 (MQ=255)
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGATCTAATTGCTGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                    >  2:538049/1‑100 (MQ=255)
   gTGAGGGCGGCCGTACCGTTTGAGCGGTCGTTGTAGCAAATGTTCTGAGCTAAGTGCCGATAACATTTGACGTAATGCGCACTAAATGGGCATCAttt                                                                                                    <  2:189216/98‑1 (MQ=37)
    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                                 >  1:2123902/1‑100 (MQ=255)
     gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                <  2:1159842/100‑1 (MQ=255)
      aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGc                                                                                               <  1:745795/100‑1 (MQ=255)
      aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATACACACTAAAAAGGCATAAtttaa                                                                                                  >  1:113997/1‑95 (MQ=25)
        ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCTGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                             >  2:1745712/1‑100 (MQ=255)
        ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                             <  1:1982361/100‑1 (MQ=255)
        ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                             <  2:402532/100‑1 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTATCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                            <  1:532631/99‑1 (MQ=38)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtt                                                                                           >  1:2011046/1‑100 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                           <  2:1698944/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTTAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                         <  1:141557/100‑1 (MQ=38)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                         <  2:1098974/100‑1 (MQ=255)
            gcCGGACCGTTGACGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                          <  2:366704/99‑1 (MQ=39)
             ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTTATTGCCGATAACATTTAACGCAATGCGCACTAAATGGGCATCAAATGATGCCCtttt                                                                                         <  2:1104517/99‑1 (MQ=38)
             ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                        <  2:1247723/100‑1 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                     <  1:1091939/100‑5 (MQ=38)
                 atcGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAAGTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtgc                                                                                    <  1:2092610/98‑6 (MQ=38)
                 aCCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCAc                                                                                    <  2:1456909/100‑1 (MQ=255)
                  ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                   <  2:1626822/100‑1 (MQ=255)
                  ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                   <  2:67504/100‑1 (MQ=255)
                    gTTGGCGCGGGTGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtgcgca                                                                                 <  1:2729166/100‑9 (MQ=37)
                    gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCt                                                                                 <  1:2374017/100‑1 (MQ=255)
                       ggCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                              <  1:2032622/100‑1 (MQ=255)
                         cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCttt                                                                               >  2:1772705/1‑97 (MQ=255)
                          gcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                            >  2:1063779/1‑99 (MQ=255)
                           cgGGCGTTGTAGCAAAAGTTATGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGAGCATCATTTGATGCCATTTTTGCACGCTTTCGTAc                                                                          <  2:2538894/100‑1 (MQ=255)
                              gCGTTGTAGCAAAAGTTCTGAGATAATTGCCGATAACATTTGACGCAATGCGAACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                            <  2:216032/95‑1 (MQ=255)
                                  tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACCTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGAACGCTTTGATACCAGAAcc                                                                   >  1:21408/1‑100 (MQ=255)
                                  tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCTCTTTTTGCACGCTTTCGTACCAGAAcc                                                                   <  2:113996/100‑1 (MQ=255)
                                  tGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                   <  1:652061/100‑1 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                                >  1:2681470/1‑100 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                                >  2:2961892/1‑100 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                                >  2:2397126/1‑100 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATAATTTGATACCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                                 >  1:1081649/1‑99 (MQ=255)
                                      tcttccgatCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                               <  2:2546236/92‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                               <  1:1881887/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                               >  2:2210291/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                               >  2:1721854/1‑100 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                               <  2:866223/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                               <  1:1453245/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGAGCACTAAAAGTGAATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                   >  2:871956/1‑96 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGTAGCAATGCGCACTAAAAAGGCATCATTTGATGCCCTTTTTGCACGCCTTCGTACCAGAACCTCGCt                                                              >  2:1977468/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCAGATAACATTTGACGCAATGCGCGCTAAAAGGGCATCATTTGATGCTCATGTTGCACGCTGTCGTACGAGAACCTgg                                                                >  2:2639686/1‑98 (MQ=255)
                                            gTTCTGATCTAATTGCCGATAACATTTGACGCAATGCGCACTGAAAGCGCATCATTTGATGCTCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcagc                                                          <  2:775547/99‑3 (MQ=255)
                                             ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                        <  2:1629446/100‑1 (MQ=255)
                                               cTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTg                                                      >  2:2074728/1‑100 (MQ=255)
                                                   gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAgatcgga                                                  >  1:1170615/1‑94 (MQ=255)
                                                   gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAttt                                                  <  2:1024789/100‑1 (MQ=255)
                                                    cTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtttt                                                 >  1:2325112/1‑100 (MQ=255)
                                                      aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                               <  1:690731/100‑1 (MQ=255)
                                                        ttGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                             <  1:2357043/100‑1 (MQ=255)
                                                         tGCTGATATCATTTGACGTAATGCGCACTAAAAGTGTATCATTTGATGTTCTTTTTGTACGCTTTCGTACCAGAACTTGGCTCATCAGTGATTTTCTTTg                                            <  2:1627486/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                            >  1:982556/1‑100 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGTGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTTTCTCATCAGTGATTATCTTTg                                            >  1:1329641/1‑99 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                             <  2:2578626/98‑1 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                           >  2:1245876/1‑100 (MQ=255)
                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                           >  2:1941426/1‑100 (MQ=255)
                                                           ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                          >  2:1997725/1‑100 (MQ=255)
                                                           ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                          <  2:2325086/100‑1 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                         <  2:2647986/100‑1 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGATCATCAGTGATTTTCTTTGTCa                                         >  2:829334/1‑100 (MQ=255)
                                                               tAACATTTGACGTAATGCGCACTAAAAAGACATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAATCTAGCTCATCAGTGATTTTCTTTGTCATa                                       >  2:748549/1‑99 (MQ=255)
                                                                aaCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAACAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAAt                                     >  2:1743461/1‑100 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                    >  1:2292884/1‑100 (MQ=255)
                                                                    ttggACGCAATGCTCACTAACAGGGCATCAGTTGATGCTCTTTTTTCATGTTTTCGGACCAGAACCTGGTTCATCAGTGATTTTCTTTGTCATAATCAtt                                 <  2:1010760/97‑1 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCt                              >  1:3028223/1‑100 (MQ=255)
                                                                         cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCAAGCTTTCGTACCAGAAACTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                            >  1:1530011/1‑100 (MQ=255)
                                                                         agCAATGCGCACTAAAAAAACATCATTTGATGCCATTTTTGCACGCTTTCGTACCAGAACCTGGCTAATCAGTGATTTTCTTTGTCATAATCATTGCTga                            >  2:1296120/2‑100 (MQ=255)
                                                                          gCAATGCGCACTAAAAGGGTATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                            <  1:1745732/99‑1 (MQ=255)
                                                                          gCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                           <  2:2823026/100‑1 (MQ=255)
                                                                            aaagcgtgcaaaaaGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                         >  2:50479/11‑100 (MQ=18)
                                                                              tGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                       >  1:138641/1‑100 (MQ=255)
                                                                              tGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                       >  2:371031/1‑100 (MQ=255)
                                                                                ctcACTAATAGGGTATCATTTGATGCCCTTTTTGTTCGTTTTCGTACTATAACCTGTGTCATAAGTGATTGTCTTTTTCATAATCATTGCTga                            <  2:2867413/91‑1 (MQ=255)
                                                                                cgtgcacacaGGGCATCAAATGATGACCTTTTTGCACGCTTTCGTACCAGAACCAGGCTCAAAAGTGATTCTCTTTGTCATAATCATTGCTga                            >  2:14289/10‑93 (MQ=255)
                                                                                 gcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                      <  2:3044939/98‑1 (MQ=255)
                                                                                 gcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                     >  1:1200032/1‑99 (MQ=255)
                                                                                  cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                   >  1:1309879/1‑100 (MQ=255)
                                                                                   aCTAAAAGGGCATCATTGTATGCCCTTTTTGTACGCTTTCGTACCAGAACCTGGCGCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                   <  2:2442296/99‑1 (MQ=255)
                                                                                    caaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAAGctct                  >  1:2242441/3‑99 (MQ=38)
                                                                                    cTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg                 >  1:2656400/1‑100 (MQ=255)
                                                                                     aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt                >  2:1808330/2‑100 (MQ=255)
                                                                                      aaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt               >  2:1245886/1‑100 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTTTTTTTTTTTGTACCAGAACCTGGCTCATCAGTGTTTTTCTTTGTCATAATCATTGCTTAGACAGGCTCTGtt               >  2:930223/1‑99 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg              <  2:646425/100‑1 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg              <  2:2294251/100‑1 (MQ=255)
                                                                                        aaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa             <  1:1263939/100‑1 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTTTTTTTTTTTTTACCAGAACCTGGCTCATCAGTTATTTTCTTTGTCTTAATCATTGCTTATACAGGCTCTGttta             >  2:731088/1‑98 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa             >  1:1246750/1‑100 (MQ=255)
                                                                                         aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg            <  1:2479866/100‑1 (MQ=255)
                                                                                           ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg          >  1:2531185/1‑100 (MQ=255)
                                                                                             catcgtTTGATGCCCTTTTTGCACGCTTTAGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg        <  2:131381/100‑1 (MQ=255)
                                                                                              atcatTTGATGCCCTTTTTGCACGCTTTCTTATCAGAACCTGGTTCATCAGTAATTTTCTTTGTGATTATCATTGCTTAGATATGCTGTGTTGAGGGCg        <  2:2451468/99‑1 (MQ=255)
                                                                                                catTTGATGCCCTTTTTGCACGCTTTCGTACCAGTTCCTGGCTCATCAGTGATTTTCTTTGTCATAATCATGGCTGAGACAGGCTCTGTTGAGGGCGtat     <  1:192186/100‑1 (MQ=255)
                                                                                                catTTGATGCCCTTTTTGCACGCTTTCGTACCAGATCCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtat     <  1:1772725/100‑1 (MQ=255)
                                                                                                 atTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtata    >  2:189089/1‑100 (MQ=255)
                                                                                                   ttGATGCCCTTTTTGCACGCTTTCGTACCAGATCCTGGCTTATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAAt  <  1:1153634/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAAT  >  NC_000913/4177073‑4177271

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