Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,171 T→A 18.9% 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,1710TA18.9% 421.0 / 94.6 187intergenic (+43/‑187)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (72/79);  new base A (18/17);  total (90/97)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-02

Only 100 of 199 total aligned reads displayed.

TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAA  >  NC_000913/4177073‑4177270
                                                                                                  |                                                                                                   
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                   <  2:2153450/100‑1 (MQ=255)
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                   >  2:1160142/1‑100 (MQ=255)
 ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                   >  1:3634435/1‑100 (MQ=255)
   gTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                 >  1:4163367/1‑100 (MQ=255)
     gAAGGCGGCCTTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                               <  1:1505936/100‑1 (MQ=255)
     gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                               <  1:278725/100‑1 (MQ=255)
     gAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                               >  1:3904944/1‑100 (MQ=255)
      aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGc                                                                                              >  1:403859/1‑100 (MQ=255)
      aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGc                                                                                              <  2:694973/100‑1 (MQ=255)
      aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGc                                                                                              <  1:2009075/100‑1 (MQ=255)
       aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                             <  2:1254978/100‑1 (MQ=255)
        ggcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                            <  1:1482366/100‑1 (MQ=255)
         gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCt                                                                                           <  2:3693835/100‑1 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtt                                                                                          >  2:4393227/1‑100 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtt                                                                                          <  1:470681/100‑1 (MQ=255)
          cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtt                                                                                          <  1:3382938/100‑1 (MQ=255)
          cggcCGTACCGTTCGCGCGGGCGTGGTCGCAACAGCTCGGAGCTAATCACCGCTAGCATTTGACGCAATGCGCACTAAACGAGCATCATTTGATGCGCtt                                                                                          >  1:1094454/1‑100 (MQ=255)
           ttccgatctgTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                         <  2:4200094/91‑1 (MQ=38)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttt                                                                                         <  1:3209090/100‑1 (MQ=255)
           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                         <  1:1971427/100‑1 (MQ=255)
            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttt                                                                                        <  2:2000675/100‑1 (MQ=255)
             ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                       >  2:1907174/1‑100 (MQ=255)
             ccGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAAATGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttta                                                                                       <  1:822686/100‑2 (MQ=38)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                      <  1:1495774/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                      <  1:753974/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                      <  1:1943205/100‑1 (MQ=255)
              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttag                                                                                      <  2:1772462/100‑3 (MQ=255)
               gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGc                                                                                     <  1:3584734/100‑1 (MQ=255)
               gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagt                                                                                     <  2:1631934/100‑4 (MQ=38)
               gTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagt                                                                                     <  2:2382251/100‑4 (MQ=38)
                tgccGTTGGCGCGAGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                    >  2:2988757/3‑100 (MQ=255)
                tctcGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                    <  2:3345116/97‑1 (MQ=255)
                tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagtg                                                                                    <  1:3367572/100‑5 (MQ=38)
                naCCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                      >  1:1020319/2‑98 (MQ=255)
                  ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                  <  2:3105028/100‑1 (MQ=255)
                   cGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCaaga                                                                                 >  1:3345099/1‑97 (MQ=255)
                       ggCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                             >  2:3877345/1‑100 (MQ=255)
                        gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                            <  2:2865436/100‑1 (MQ=255)
                          gcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTa                                                                          <  2:3159947/100‑1 (MQ=255)
                             ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCa                                                                       >  2:3752032/1‑100 (MQ=255)
                               cGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGa                                                                     >  2:3430434/1‑100 (MQ=255)
                                gTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGaa                                                                    >  2:4226420/1‑100 (MQ=255)
                                     aGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                               <  1:4311977/100‑1 (MQ=255)
                                      gCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                              >  1:2469396/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                             >  1:176615/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                             >  1:1233605/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                             >  1:2261769/1‑100 (MQ=255)
                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                             <  1:863260/100‑1 (MQ=255)
                                        aaaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtc                                                            >  1:1656750/1‑100 (MQ=255)
                                         aaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtca                                                           >  1:71930/1‑100 (MQ=255)
                                           aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                         <  1:1917191/100‑1 (MQ=255)
                                            gTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                        <  2:71930/100‑1 (MQ=255)
                                            gTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                        <  1:2373094/100‑1 (MQ=255)
                                             ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                       <  2:3478115/100‑1 (MQ=255)
                                                tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                    >  2:2329147/1‑100 (MQ=255)
                                                 gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                   >  1:2976795/1‑100 (MQ=255)
                                                       aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGCTGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                             >  1:1164294/1‑100 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                           >  1:3076171/1‑100 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                           <  2:2785215/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                           <  2:2586557/100‑1 (MQ=255)
                                                         tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                           <  2:1811181/100‑1 (MQ=255)
                                                           ccGATAACATTTGACGCAATGCGCACTAAAAGGGGATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGTTTTCCTTTGTc                                         <  1:2310489/100‑1 (MQ=255)
                                                            cGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                        <  2:1047269/100‑1 (MQ=255)
                                                             gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAt                                       >  2:2047752/1‑100 (MQ=255)
                                                             gATAACATTGGACGCAATGCGCACTAAGAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                        <  2:2624502/99‑1 (MQ=255)
                                                              aTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATa                                      >  2:4154897/1‑100 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                   <  1:371573/100‑1 (MQ=255)
                                                                 aCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATc                                   <  2:1074834/100‑1 (MQ=255)
                                                                  cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCa                                  >  2:600935/1‑100 (MQ=255)
                                                                  cATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCa                                  <  1:3545354/100‑1 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                                 >  1:3775704/1‑100 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                                 >  2:3529698/1‑100 (MQ=255)
                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAt                                 >  2:1360077/1‑100 (MQ=255)
                                                                    tttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAtt                                <  2:3121056/100‑1 (MQ=255)
                                                                     ttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTg                               <  2:289557/100‑1 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                              <  1:4103069/99‑1 (MQ=255)
                                                                       gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCACCAGTGATTTTCTTTGTCATAATCATTGCt                             >  1:2653294/1‑100 (MQ=255)
                                                                          gCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgag                          >  2:1550666/1‑100 (MQ=255)
                                                                            aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                        <  1:2288937/100‑1 (MQ=255)
                                                                              tGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                      >  1:1381392/1‑100 (MQ=255)
                                                                               gcgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCGTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                     <  1:5226/100‑1 (MQ=255)
                                                                                cgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCCGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                    >  2:2110152/1‑100 (MQ=255)
                                                                                 gcACTAAAAGGGCATAATTTGATGCCCTTTTTGCACGCATTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                   >  1:3488227/1‑100 (MQ=255)
                                                                                  tgcaaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                  >  2:460183/5‑100 (MQ=38)
                                                                                    caaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTg                >  2:2773505/3‑100 (MQ=255)
                                                                                    cTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgagagatcggaag                >  2:4039622/1‑91 (MQ=255)
                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg             <  1:3578693/100‑1 (MQ=255)
                                                                                       aaaGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt              >  1:1581290/1‑99 (MQ=37)
                                                                                        aaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGGGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg             >  2:2715213/1‑99 (MQ=255)
                                                                                        aaGGGCATCATTTGATGCCCTTTTTGCACGCTATAGTAAGAGGACCAGGATGAACAGTCCTTTTCTTTGTCAGACACATTGCTGAACCAGGCTCTGttaa            >  1:421098/1‑98 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa            >  1:2512157/1‑100 (MQ=255)
                                                                                        aaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa            >  2:3790514/1‑100 (MQ=255)
                                                                                         aGGGCATCATTTGATGCCCTTTTTGTACGCGTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAACCATTGCTGAGACAGGCTCTGTTGAg           <  1:26780/100‑1 (MQ=255)
                                                                                          gggCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg          >  1:3963976/1‑100 (MQ=255)
                                                                                           ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg         <  1:3993014/100‑1 (MQ=255)
                                                                                            gCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc        >  2:3650109/1‑100 (MQ=255)
                                                                                             catcatTTGATGCCCTTTTTGTACGCTTTCGTACCAGAACCTGGCTCATTAGTGATTTTCTTTTTCATAATCATTGCTGAGACAGGACCTGTTGAg           >  1:1865154/1‑96 (MQ=255)
                                                                                              atcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt      >  1:806828/1‑100 (MQ=255)
                                                                                               tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGta     >  2:1854219/1‑100 (MQ=255)
                                                                                                  tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATaa  <  2:493115/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAA  >  NC_000913/4177073‑4177270

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