Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,361 T→G 16.3% intergenic (‑32/+39) 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,3610TG16.3% 291.0 / 55.9 124intergenic (‑32/+39)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base T (46/57);  new base G (14/6);  total (61/63)
Fisher's exact test for biased strand distribution p-value = 5.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-02

Only 100 of 126 total aligned reads displayed.

TGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACT  >  NC_000913/3471262‑3471458
                                                                                                   |                                                                                                 
tGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCtt                                                                                                   <  1:798837/100‑1 (MQ=255)
  gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                                 >  1:2247428/1‑100 (MQ=255)
  gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACGATAGTGCTCTCCCCTTCa                                                                                                 <  2:578998/100‑1 (MQ=255)
   tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                                >  1:2910123/1‑100 (MQ=255)
        gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                           <  1:454345/100‑1 (MQ=255)
          ggCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgaga                                                                                            >  1:125282/1‑97 (MQ=255)
           gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGc                                                                                        <  2:2449572/100‑1 (MQ=255)
            ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                       >  2:853928/1‑100 (MQ=255)
             cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                      <  2:2373576/100‑1 (MQ=25)
             cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                      <  2:2795188/100‑1 (MQ=25)
               aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                    <  2:2265659/100‑1 (MQ=255)
               aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                    <  2:1312596/100‑1 (MQ=255)
                tAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATGTTCCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                     <  1:948795/98‑2 (MQ=25)
                  gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGAc                                                                                 <  1:2233252/100‑1 (MQ=255)
                   tgccAACGTTAACGTGCGGTTTTGTACGTTCACTTTTTTCTTTAACCACGGCTATCATCCTTACTATAGTGCTCTCCCCTTCAGAAGTCAGCACgttac                                                                                 >  1:1859283/3‑95 (MQ=21)
                   tACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatagt                                                                                <  2:742729/100‑7 (MQ=25)
                   tACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACt                                                                                >  2:435086/1‑100 (MQ=255)
                    aCCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACtt                                                                               >  1:2287644/1‑100 (MQ=255)
                    aCCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACtt                                                                               >  2:1643716/1‑100 (MQ=255)
                     ccAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACtt                                                                               >  1:1714480/1‑99 (MQ=255)
                      cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCATGAGAGAGCACGGGACTTTg                                                                             <  1:563987/100‑1 (MQ=255)
                      cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                             <  2:351645/100‑1 (MQ=255)
                       aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  2:1417516/100‑1 (MQ=255)
                       aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  2:1879856/100‑1 (MQ=255)
                       aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  1:334624/100‑1 (MQ=255)
                       aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  1:1113375/100‑1 (MQ=255)
                        aCGTTAACGTGTTGTTTTGTTCGTTCAAATTATAAATTAAAAAAAGATATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                           <  1:1759235/100‑1 (MQ=255)
                        aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                           <  1:868680/100‑1 (MQ=255)
                        aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCGATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                           <  1:1176569/100‑1 (MQ=255)
                         cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  1:334714/98‑1 (MQ=255)
                         cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                          >  2:406370/1‑100 (MQ=255)
                          gTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTCGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAt                                                                         <  1:2719822/100‑1 (MQ=255)
                          gTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAt                                                                         <  2:2829742/100‑1 (MQ=255)
                          gTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAt                                                                         <  2:1487977/100‑1 (MQ=255)
                                gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                   >  2:2022137/1‑100 (MQ=255)
                                gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                   >  2:1151741/1‑100 (MQ=255)
                                gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                   >  2:1465090/1‑100 (MQ=255)
                                   cGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTa                                                                <  2:2657413/100‑1 (MQ=255)
                                    ggTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAg                                                               >  2:2106978/1‑100 (MQ=255)
                                     gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGATAGAGCACGGGACTTTGGTATTAACCCTTAgg                                                              >  2:2175614/1‑100 (MQ=255)
                                     gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAgg                                                              <  2:772803/100‑1 (MQ=255)
                                         tGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTa                                                          >  1:164342/1‑100 (MQ=255)
                                         tGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTa                                                          >  2:2529158/1‑100 (MQ=255)
                                              gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                     <  1:2855966/100‑1 (MQ=255)
                                              gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAAAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                     >  1:1582939/1‑100 (MQ=255)
                                               ttCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAcc                                                    <  1:1072030/100‑1 (MQ=255)
                                                 cAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCa                                                   >  2:2911086/1‑99 (MQ=255)
                                                   aaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTTTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACg                                                 >  1:192655/1‑99 (MQ=255)
                                                   aaTTTTTCTTTAGACACGGCTATATTCCTTACGATAGTGATCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                                <  1:2050635/100‑1 (MQ=255)
                                                    aTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                               <  2:1394420/100‑1 (MQ=255)
                                                    aTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                                <  1:1410142/99‑1 (MQ=255)
                                                      ttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCt                                             <  2:137455/100‑1 (MQ=255)
                                                      ttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGACAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCt                                             <  1:1612479/100‑1 (MQ=255)
                                                        ttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAATCCTTAGGCTTATTTACCACGGGCTTc                                           <  2:1776576/100‑1 (MQ=255)
                                                          cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCaa                                         >  1:2283561/1‑100 (MQ=255)
                                                           tttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAAt                                        >  2:729775/1‑100 (MQ=255)
                                                                acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                   <  1:866035/100‑1 (MQ=255)
                                                                acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                   <  2:570713/100‑1 (MQ=255)
                                                                 cgcGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                  >  2:1092744/3‑100 (MQ=255)
                                                                  acGGCTATATTCCTTACTATAGTGCTCTCCCCTTCTGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACg                                    <  1:1847797/97‑1 (MQ=255)
                                                                      cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGATAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                             <  1:1303746/100‑1 (MQ=255)
                                                                      cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                             >  2:467487/1‑100 (MQ=255)
                                                                      cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                             >  2:656320/1‑100 (MQ=255)
                                                                           ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                        >  2:639892/1‑100 (MQ=255)
                                                                            tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                       >  1:59068/1‑100 (MQ=255)
                                                                            tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                       <  1:424070/100‑1 (MQ=255)
                                                                            tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                       >  1:162229/1‑100 (MQ=255)
                                                                              cTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGTACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTa                     >  2:1235134/1‑100 (MQ=255)
                                                                               ttACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTa                     <  2:1163484/99‑1 (MQ=255)
                                                                                 aCTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTc                  <  1:1495846/100‑1 (MQ=255)
                                                                                 aCTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTc                  <  1:36896/100‑1 (MQ=255)
                                                                                  cTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCg                 <  2:2374812/100‑1 (MQ=255)
                                                                                  cTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCg                 >  1:677832/1‑100 (MQ=255)
                                                                                   tcccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCgg                <  1:247598/96‑1 (MQ=37)
                                                                                    cccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               >  2:2892252/4‑100 (MQ=38)
                                                                                    ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGTCTTCAATTACGGCCTGAGCAACGTTACTCGgc               >  1:2488927/1‑100 (MQ=255)
                                                                                    ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               <  2:1035282/100‑1 (MQ=255)
                                                                                    ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               <  1:2257877/100‑1 (MQ=255)
                                                                                     ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  1:488687/3‑100 (MQ=38)
                                                                                     ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  2:2033924/3‑100 (MQ=38)
                                                                                     ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  1:2227001/3‑100 (MQ=38)
                                                                                     ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               >  1:60993/3‑99 (MQ=38)
                                                                                      cgTGCTCTCTCCTGAGGGGGAGAGCACGGGACTTTGGTATTACCCCTTAGGCTTGTGTACCACGGGCTTCAATGACGGCCTGAACAACGTTACTCGgcgc             >  1:1698401/2‑100 (MQ=255)
                                                                                      cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  2:2285654/2‑100 (MQ=38)
                                                                                      cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  2:497726/2‑100 (MQ=38)
                                                                                      cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  1:1476327/2‑100 (MQ=38)
                                                                                      cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  2:473132/2‑100 (MQ=38)
                                                                                       gTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  1:1659630/1‑98 (MQ=39)
                                                                                       gTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt            >  1:2834228/1‑100 (MQ=39)
                                                                                         gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCtt           >  1:879077/1‑99 (MQ=255)
                                                                                         gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc          <  2:2284327/100‑1 (MQ=255)
                                                                                            ctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc       <  1:2355849/100‑1 (MQ=255)
                                                                                              ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat     >  1:498194/1‑100 (MQ=255)
                                                                                              ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat     <  1:2897866/100‑1 (MQ=255)
                                                                                              ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat     <  1:650437/100‑1 (MQ=255)
                                                                                              ccccTTCAGGAGAGAGCACGGGACGTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc       <  2:2668909/98‑1 (MQ=255)
                                                                                                ccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATAc   <  2:1979926/100‑1 (MQ=255)
                                                                                                 cTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACt  <  2:492414/100‑1 (MQ=255)
                                                                                                 cTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACt  <  2:795294/100‑1 (MQ=255)
                                                                                                  ttCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATAc   <  1:2505000/98‑1 (MQ=255)
                                                                                                   |                                                                                                 
TGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACT  >  NC_000913/3471262‑3471458

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: