Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,357 C→T 15.3% intergenic (‑28/+43) tufA ← / ← fusA translation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,471,3570CT15.3% 552.9 / 72.7 220intergenic (‑28/+43)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (100/85);  new base T (16/18);  total (117/103)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

Only 100 of 232 total aligned reads displayed.

ACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATA  >  NC_000913/3471258‑3471456
                                                                                                   |                                                                                                   
aCCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcc                                                                                                     <  2:3550716/100‑1 (MQ=255)
aCCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGGTGTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCGTTACTATAGTGATCTcc                                                                                                     >  1:449663/1‑100 (MQ=21)
 ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTccc                                                                                                    >  1:598305/1‑100 (MQ=255)
 ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcc                                                                                                     <  2:89114/99‑1 (MQ=255)
  cGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcccc                                                                                                   <  1:1206675/100‑1 (MQ=255)
     ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                >  2:1081520/1‑100 (MQ=255)
      gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                               >  2:2242176/1‑100 (MQ=255)
      gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                               >  1:1680433/1‑100 (MQ=255)
       tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                              >  1:1941779/1‑100 (MQ=255)
        cAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAgg                                                                                             >  1:3759657/1‑100 (MQ=255)
         aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGga                                                                                            <  2:1534556/100‑1 (MQ=255)
          aCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                           >  2:1644417/1‑100 (MQ=255)
           cGTGGCCGATAGTACCAACGTTAACGTGCGGTTATGTACGTTCAAATTTTTCTTTAGACACGGCTATATTACTTACTATAGTGCTCTCCCCTTCAGgaga                                                                                          >  1:2651704/1‑100 (MQ=255)
            gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                         <  2:3764256/100‑1 (MQ=255)
             tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                        <  2:781663/100‑1 (MQ=255)
             tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                        <  2:1015971/100‑1 (MQ=255)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGc                                                                                      <  2:4347071/100‑1 (MQ=25)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGc                                                                                      <  2:2366776/100‑1 (MQ=25)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGc                                                                                      <  2:1166545/100‑1 (MQ=255)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGc                                                                                      <  2:642030/100‑1 (MQ=255)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGc                                                                                      <  1:958744/100‑1 (MQ=255)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACCATAGTGCTCTCTCCTGAAGGGGAGAGc                                                                                      <  1:3865157/100‑1 (MQ=21)
                ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                     >  2:3705059/1‑100 (MQ=255)
                  gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                   <  1:1456242/100‑2 (MQ=25)
                  gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                   <  1:3020998/100‑2 (MQ=25)
                  gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACg                                                                                   >  2:1803862/1‑100 (MQ=255)
                    tAGTACCAATGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                   <  1:4035638/98‑2 (MQ=21)
                    tAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                 >  2:512070/1‑100 (MQ=255)
                    tAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                 <  1:2250209/100‑1 (MQ=255)
                      gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatag                                                                               <  1:2229040/100‑6 (MQ=25)
                      gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGAc                                                                               <  2:487077/100‑1 (MQ=255)
                         ccAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                            >  2:4025941/1‑100 (MQ=255)
                         ccAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                            >  1:945740/1‑100 (MQ=255)
                              gTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAt                                                                       >  2:3412615/1‑100 (MQ=255)
                               ttAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAtt                                                                      <  1:2748102/100‑1 (MQ=255)
                                tAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGGGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTa                                                                     <  2:1991277/100‑1 (MQ=255)
                                 aaCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTaa                                                                    <  2:1320355/100‑1 (MQ=255)
                                  aCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAc                                                                   <  2:2557357/100‑1 (MQ=255)
                                  aCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAc                                                                   <  2:3333493/100‑1 (MQ=255)
                                    gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                 >  2:3838207/1‑100 (MQ=255)
                                     tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                >  2:2803663/1‑100 (MQ=255)
                                     tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                >  1:3496688/1‑100 (MQ=255)
                                       cGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTa                                                              >  2:1575671/1‑100 (MQ=255)
                                       cGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTa                                                              >  1:2450017/1‑100 (MQ=255)
                                        ggTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTa                                                              >  2:3163774/1‑99 (MQ=255)
                                        ggTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAg                                                             <  1:2524149/100‑1 (MQ=255)
                                           tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTTTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCt                                                          >  1:983610/1‑100 (MQ=255)
                                              gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                       >  1:1712413/1‑100 (MQ=255)
                                              gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                       <  2:3282110/100‑1 (MQ=255)
                                               tACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAtt                                                      <  2:632480/100‑1 (MQ=255)
                                                 cGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTa                                                    >  2:3053699/1‑100 (MQ=255)
                                                 cGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTa                                                    >  2:2528817/1‑100 (MQ=255)
                                                     cAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCCCGGGACTTTGGTATTAACCCTTAGGCTTATCTCCCAc                                                >  1:562949/1‑100 (MQ=255)
                                                     cAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCAc                                                <  1:621378/100‑1 (MQ=255)
                                                     cAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCAc                                                <  1:4198495/100‑1 (MQ=255)
                                                      aaaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACg                                               >  1:3391781/1‑100 (MQ=255)
                                                        aTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                             <  2:277304/100‑1 (MQ=255)
                                                        aTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                             >  2:3360343/1‑100 (MQ=255)
                                                        aTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                              >  1:932548/1‑99 (MQ=255)
                                                         tttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                            >  2:1405128/1‑100 (MQ=255)
                                                            ttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTc                                         <  1:862020/100‑1 (MQ=255)
                                                               tttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAAt                                      >  2:2895196/1‑100 (MQ=255)
                                                               tttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAAt                                      >  2:574576/1‑100 (MQ=255)
                                                                 tAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTa                                    <  2:3648521/100‑1 (MQ=255)
                                                                  aGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTAc                                   >  2:831224/1‑100 (MQ=255)
                                                                   gACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACg                                  >  2:2235199/1‑100 (MQ=255)
                                                                     cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                >  1:1503945/1‑100 (MQ=255)
                                                                       cGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCt                              >  2:1822719/1‑100 (MQ=255)
                                                                       cGGCTAGATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCt                              >  2:1197860/1‑100 (MQ=255)
                                                                        ggCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCATTAGGCTTATTTACCACGGGCTTCAATTACg                                  >  2:3093218/1‑95 (MQ=255)
                                                                          cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                           >  1:799248/1‑100 (MQ=255)
                                                                          cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                           >  1:3604074/1‑100 (MQ=255)
                                                                           tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                          >  2:4228966/1‑100 (MQ=255)
                                                                           tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                          <  2:2289448/100‑1 (MQ=255)
                                                                              atTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACAACGGGCTTCAATTACGGCCTGAGCAAc                       >  1:4221183/1‑100 (MQ=255)
                                                                                tCCTTACTATAGTGCTCTCCCCTTTAGTAGAGAGCACGGGTCTTTGTTATTAACCCTTCGGCTTATTTACCACGGGCTTCAATTACGGCCTGTGCAACg                      <  1:1987067/99‑1 (MQ=255)
                                                                                 ccTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCCCCGGGCTTTGGTATTCACCCTTAGGCTTATTTACCACGGTCTTCCATTACGGCCTTAGCAACGtt                    >  1:1908296/1‑100 (MQ=255)
                                                                                  cTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTa                   <  2:2127149/100‑1 (MQ=255)
                                                                                  cTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGCTATTAACACTTAGGCTTATTTACCACGGGCTTCAATTACAGCCTGAGCAACGtt                    >  1:3897481/1‑99 (MQ=255)
                                                                                      cTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCg               >  2:1279248/1‑100 (MQ=255)
                                                                                       tataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTGGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTc                <  1:2860221/98‑1 (MQ=255)
                                                                                        cccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             >  2:3890173/4‑100 (MQ=38)
                                                                                        ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             >  2:4091919/1‑100 (MQ=255)
                                                                                        ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             <  1:3695949/100‑1 (MQ=255)
                                                                                         ccgTGCTCTCTCCTGAAGGGGAGAGCACTGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg            >  1:4089907/3‑100 (MQ=37)
                                                                                         ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg            >  1:1239933/3‑100 (MQ=38)
                                                                                          cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAACAACGTTACACGgcgc           >  1:2324790/2‑100 (MQ=38)
                                                                                          cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAAGCTTATTTACCACGGGCTTCAATTACGGCCTAAGCAACGTTACTCGgc             >  1:315797/2‑98 (MQ=37)
                                                                                           gTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt          >  1:3353836/1‑100 (MQ=39)
                                                                                           gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt          >  2:4169759/1‑100 (MQ=255)
                                                                                           gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt          <  1:1009266/100‑1 (MQ=255)
                                                                                            tGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCtt         <  2:1941785/100‑1 (MQ=255)
                                                                                             gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc        >  1:1317004/1‑100 (MQ=255)
                                                                                             gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTCAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc        >  2:4380928/1‑100 (MQ=255)
                                                                                               tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat      <  1:1192948/100‑1 (MQ=255)
                                                                                               tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat      >  1:2535970/1‑100 (MQ=255)
                                                                                                ctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc     <  2:586932/100‑1 (MQ=255)
                                                                                                  ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat   >  1:4156614/1‑100 (MQ=255)
                                                                                                   cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat      <  1:1654484/96‑1 (MQ=255)
                                                                                                   cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat   >  2:1740456/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
ACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATA  >  NC_000913/3471258‑3471456

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: