Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,361 T→G 20.0% 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,3610TG20.0% 351.0 / 90.1 161intergenic (‑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 (57/71);  new base G (17/15);  total (75/86)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-02

Only 100 of 161 total aligned reads displayed.

GGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACT  >  NC_000913/3471263‑3471458
                                                                                                  |                                                                                                 
ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTGCGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                  <  2:2524470/100‑1 (MQ=255)
ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                  <  1:2758151/100‑1 (MQ=255)
ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                  >  1:3486550/1‑100 (MQ=255)
ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAATGCTCTCCCCTTc                                                                                                  >  2:2348504/1‑100 (MQ=255)
 gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                                 <  1:824660/100‑1 (MQ=255)
 gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                                 <  1:1846890/100‑1 (MQ=255)
  tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                                <  1:2316982/100‑1 (MQ=255)
  tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                                <  1:1910880/100‑1 (MQ=255)
     aCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                             <  1:1456832/100‑1 (MQ=255)
       gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                           <  2:3222785/100‑1 (MQ=255)
       gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                           <  1:2824011/100‑1 (MQ=255)
         ggCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                         <  1:1764121/100‑1 (MQ=255)
         ggCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                         >  1:3798810/1‑100 (MQ=255)
          gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATTGTGCTCTCTCCTGAAGGGGAGAGc                                                                                        <  2:2039198/100‑1 (MQ=21)
          gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGc                                                                                        <  1:2458295/100‑1 (MQ=255)
           ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                       >  1:2962559/1‑100 (MQ=255)
           ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                       >  1:2808434/1‑100 (MQ=255)
            cGATAGTACCAATGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                      <  2:921/100‑1 (MQ=21)
            cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                      <  1:34599/100‑1 (MQ=25)
            cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCAc                                                                                      <  2:2713132/100‑1 (MQ=255)
            cGATAGTACCAACGTTAACGTGCGGTTTCGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTCCAGGAGAGAGCAc                                                                                      >  2:2050323/1‑100 (MQ=255)
             gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                     <  1:291629/100‑2 (MQ=25)
             gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                     <  1:3503774/100‑2 (MQ=25)
              aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActa                                                                                    <  2:1495924/100‑3 (MQ=25)
              aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                    <  2:2869139/100‑1 (MQ=255)
                 gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatag                                                                                 <  2:3350211/100‑6 (MQ=25)
                 gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatag                                                                                 <  2:700061/100‑6 (MQ=25)
                      caCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  2:3459533/99‑1 (MQ=255)
                      aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  2:1026077/100‑1 (MQ=255)
                      aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            >  2:1195086/1‑100 (MQ=255)
                      aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  1:585029/100‑1 (MQ=255)
                      aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                             >  2:3137563/1‑99 (MQ=255)
                       aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                            <  2:347936/99‑1 (MQ=255)
                       aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                           <  1:2707474/100‑1 (MQ=255)
                       aCGTTAACGTGCGGTTTTGTACGGTCAAATTTTTCTTGAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatagtaa                                                                              <  2:3197945/97‑9 (MQ=25)
                              cGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                    >  1:3226241/1‑100 (MQ=255)
                               gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                   <  2:756245/100‑1 (MQ=255)
                                tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                  >  2:291073/1‑100 (MQ=255)
                                    gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAgg                                                              <  1:3793869/100‑1 (MQ=255)
                                     ttttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGc                                                             <  2:139032/100‑1 (MQ=255)
                                      tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCt                                                            >  2:1655112/1‑100 (MQ=255)
                                      tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCt                                                            <  1:2206273/100‑1 (MQ=255)
                                      tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCCTAGGCt                                                            >  2:617779/1‑100 (MQ=255)
                                         gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                         <  2:3659934/100‑1 (MQ=255)
                                         gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                         <  1:1655121/100‑1 (MQ=255)
                                         gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                         <  2:708408/100‑1 (MQ=255)
                                          tACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                         <  2:2345594/99‑1 (MQ=255)
                                           aCGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAttt                                                       <  2:2566725/100‑1 (MQ=255)
                                           aCGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAttt                                                       >  1:3079216/1‑100 (MQ=255)
                                            cGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTa                                                      <  1:1755718/100‑1 (MQ=255)
                                             gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                     >  1:111046/1‑100 (MQ=255)
                                             gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                     >  2:3774436/1‑100 (MQ=255)
                                              ttCAGATTTTTCTTTAGACACGGCTATAGTCCTTACTATAGTGCTCTCTCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCGTTAGGCTTAttt                                                       <  2:1346902/97‑1 (MQ=255)
                                              ttCAAATTTTTCTTTACACACGGCTATATTCCTTACTATATTTCTTTGCCCTTCACGAGAGTGCACGTGACTTTGGTTTTAACCCTTAGTCTTATTTAcc                                                    >  2:2841753/1‑100 (MQ=255)
                                                cAAATTTTTCTTTAGACACGGCTATATTCCTTACTATATTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCa                                                   >  2:1956955/1‑99 (MQ=255)
                                                    tttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                              <  1:192330/100‑1 (MQ=255)
                                                    tttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                              <  1:2332668/100‑1 (MQ=255)
                                                      tttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCtt                                            >  1:3863442/1‑100 (MQ=255)
                                                      tttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCtt                                            <  2:1657537/100‑1 (MQ=255)
                                                           ntAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAAtt                                       >  1:216678/2‑100 (MQ=255)
                                                               acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                   <  1:3335414/100‑1 (MQ=255)
                                                                cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                  >  1:903640/1‑100 (MQ=255)
                                                                cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                  >  1:137254/1‑100 (MQ=255)
                                                                cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                  >  1:2611594/1‑100 (MQ=255)
                                                                cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                  >  2:1059489/1‑100 (MQ=255)
                                                                    gCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGa                              >  2:1062305/1‑100 (MQ=255)
                                                                     cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                             <  1:3360191/100‑1 (MQ=255)
                                                                      tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                            <  2:878202/100‑1 (MQ=255)
                                                                      tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                            >  1:1023873/1‑100 (MQ=255)
                                                                          ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                        >  1:30074/1‑100 (MQ=255)
                                                                              ttACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTAc                    >  2:2193/1‑100 (MQ=255)
                                                                               tACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTCTTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTAc                    >  2:151634/1‑99 (MQ=255)
                                                                                aCTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTc                  <  2:3047489/100‑1 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               <  2:880553/100‑1 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               <  1:1062312/100‑1 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               >  2:1326057/1‑100 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc               <  1:551949/100‑1 (MQ=255)
                                                                                    taGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  1:991823/1‑100 (MQ=255)
                                                                                    ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  2:176210/3‑100 (MQ=38)
                                                                                    ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  1:2357771/3‑100 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  1:3824932/2‑100 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  1:2932642/2‑100 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  2:354777/2‑100 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc             >  1:3285470/2‑100 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg              >  2:534821/2‑99 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCCgcgc             >  2:2451929/2‑100 (MQ=38)
                                                                                     cgTGCTCTCTCCTGAAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCATATTTACCACGGGCGTCAATTACGGCCTGAGCAACGTTACTTGgcgc             >  2:717008/2‑100 (MQ=255)
                                                                                      gTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt            >  2:1829416/1‑100 (MQ=38)
                                                                                      gTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt            >  2:1029166/1‑100 (MQ=38)
                                                                                      gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt            <  1:3613568/100‑1 (MQ=255)
                                                                                      gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt            >  1:3259278/1‑100 (MQ=255)
                                                                                       tGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCtt           >  1:2337144/1‑100 (MQ=255)
                                                                                         ctctCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca         <  1:2491433/100‑1 (MQ=255)
                                                                                          tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat        >  2:2368069/1‑100 (MQ=255)
                                                                                            tcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatca      <  1:3774464/100‑1 (MQ=255)
                                                                                             ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat     >  1:1160525/1‑100 (MQ=255)
                                                                                             ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat     <  1:2887215/100‑1 (MQ=255)
                                                                                             ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat     >  2:3603098/1‑100 (MQ=255)
                                                                                              ctcTTCAGGAGAGCGCGCGGGACAGTGGTATTAACCCTTAGGCTTATTGCCCACGAACTTGGATTACGTCTTGAGCAACGTTACTCTGCGCttactc       >  2:3827808/3‑93 (MQ=255)
                                                                                              cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa    >  1:2634033/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
GGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACT  >  NC_000913/3471263‑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: