Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,362 C→A 20.1% intergenic (‑33/+38) 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,3620CA20.1% 393.0 / 79.7 157intergenic (‑33/+38)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (56/69);  new base A (17/15);  total (73/84)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 163 total aligned reads displayed.

GGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCA  >  NC_000913/3471263‑3471461
                                                                                                   |                                                                                                   
ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                     >  1:3486550/1‑100 (MQ=255)
ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAATGCTCTCCCCTTc                                                                                                     >  2:2348504/1‑100 (MQ=255)
 gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                                    <  1:824660/100‑1 (MQ=255)
  tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                                   <  1:1910880/100‑1 (MQ=255)
   cAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAgg                                                                                                  <  1:982082/100‑1 (MQ=255)
     aCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                                <  1:1456832/100‑1 (MQ=255)
     aCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                                >  1:3797427/1‑100 (MQ=255)
       gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                              <  2:3198728/100‑1 (MQ=255)
       gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                              <  2:1324024/100‑1 (MQ=255)
        tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                             <  2:1764466/100‑1 (MQ=255)
         ggCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                            <  1:1764121/100‑1 (MQ=255)
         ggCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                            >  1:3798810/1‑100 (MQ=255)
          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)
            cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                         <  2:2241940/100‑1 (MQ=25)
            cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCAc                                                                                         <  2:2713132/100‑1 (MQ=255)
            cGATAGTACCAACGTTAACGTGCGGTTTCGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTCCAGGAGAGAGCAc                                                                                         >  2:2050323/1‑100 (MQ=255)
             gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                        <  2:1107264/100‑2 (MQ=25)
             gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                        <  1:291629/100‑2 (MQ=25)
             gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                        <  1:292209/100‑2 (MQ=25)
              aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActa                                                                                       <  2:1495924/100‑3 (MQ=25)
              aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                       <  2:2869139/100‑1 (MQ=255)
                 gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatag                                                                                    <  2:700061/100‑6 (MQ=25)
                 gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatag                                                                                    <  2:3350211/100‑6 (MQ=25)
                     cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                                >  2:254264/1‑100 (MQ=255)
                      aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                               <  2:1026077/100‑1 (MQ=255)
                      aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                                >  2:3137563/1‑99 (MQ=255)
                       aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                               <  2:347936/99‑1 (MQ=255)
                       aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                              <  1:2702819/100‑1 (MQ=255)
                       aCGTTAACGTGCGGTTTTGTACGGTCAAATTTTTCTTGAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActatagtaa                                                                                 <  2:3197945/97‑9 (MQ=25)
                        cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                             <  2:2355282/100‑1 (MQ=255)
                         gTTAACGTGCGGTTTTTTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAt                                                                            <  2:3592549/100‑1 (MQ=255)
                          ttAACGTGCGGTTTTGTACGTTCAAATCTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                             <  2:1733722/98‑1 (MQ=255)
                           tAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTa                                                                          <  2:1370918/100‑1 (MQ=255)
                              cGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                       <  1:2819795/100‑1 (MQ=255)
                              cGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                       >  1:3226241/1‑100 (MQ=255)
                              cGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                       <  2:1989960/100‑1 (MQ=255)
                              cGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                       <  1:2330691/100‑1 (MQ=255)
                               gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                      <  2:756245/100‑1 (MQ=255)
                               gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                      >  1:2345606/1‑100 (MQ=255)
                                tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                     <  2:1593935/100‑1 (MQ=255)
                                tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                     >  1:2470199/1‑100 (MQ=255)
                                tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                     >  2:291073/1‑100 (MQ=255)
                                 gCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCtt                                                                    >  2:2782183/1‑100 (MQ=255)
                                    gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAgg                                                                 >  1:1375083/1‑100 (MQ=255)
                                      tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCt                                                               <  1:2206273/100‑1 (MQ=255)
                                      tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTACGCTTCAGGAGAGAGCACGGTACTTTGGTATTAACCCTTAg                                                                  >  2:558325/1‑97 (MQ=255)
                                         gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                            <  2:708408/100‑1 (MQ=255)
                                         gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAt                                                            <  2:3659934/100‑1 (MQ=255)
                                           aCGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAttt                                                          >  1:3079216/1‑100 (MQ=255)
                                             gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                        >  2:3774436/1‑100 (MQ=255)
                                             gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                        >  1:111046/1‑100 (MQ=255)
                                                cAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCAc                                                     >  1:3251860/1‑100 (MQ=255)
                                                    tttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                 <  1:2497568/100‑1 (MQ=255)
                                                      tttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCtt                                               <  2:1657537/100‑1 (MQ=255)
                                                      tttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCtt                                               >  1:3863442/1‑100 (MQ=255)
                                                               acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACgg                                      <  1:3335414/100‑1 (MQ=255)
                                                                 acGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGcc                                    <  1:3431653/100‑1 (MQ=255)
                                                                    gCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGa                                 >  2:1062305/1‑100 (MQ=255)
                                                                     cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                                >  1:2997831/1‑100 (MQ=255)
                                                                      tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                               >  1:1023873/1‑100 (MQ=255)
                                                                      tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                               <  1:2303620/100‑1 (MQ=255)
                                                                      tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                               <  2:878202/100‑1 (MQ=255)
                                                                           tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                          <  2:1504263/100‑1 (MQ=255)
                                                                               tACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACt                      <  1:2019649/100‑1 (MQ=255)
                                                                                 gtcccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCg                    >  1:383564/6‑100 (MQ=37)
                                                                                  tcccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCgg                   >  1:3318898/5‑100 (MQ=37)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                  >  2:2144388/1‑100 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                  <  2:880553/100‑1 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                  <  1:551949/100‑1 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                  >  2:1004469/1‑100 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                  >  2:1326057/1‑100 (MQ=255)
                                                                                   ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                  <  1:1062312/100‑1 (MQ=255)
                                                                                    ccgTGCTCTCTCCTGAAGGTGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg                 >  1:183453/3‑100 (MQ=255)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc                >  1:1897337/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)
                                                                                     cgTGCTCTCTCCTGAAGGGGAGAGCACGGGAATTTTGTATTAGCCCTTAGGGGTATTTGACACGGGCGTCAATTACGGCCTGAGCAACGTTACTCGgcgc                >  2:2452127/2‑100 (MQ=21)
                                                                                     cgTGCTCTCTCCTGAAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCATATTTACCACGGGCGTCAATTACGGCCTGAGCAACGTTACTTGgcgc                >  2:717008/2‑100 (MQ=255)
                                                                                      gTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt               >  2:1829416/1‑100 (MQ=38)
                                                                                      gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt               <  1:3613568/100‑1 (MQ=255)
                                                                                      gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt               >  1:2752679/1‑100 (MQ=255)
                                                                                       tGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCtt              >  2:3397660/1‑100 (MQ=255)
                                                                                       tGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCtt              >  2:1042747/1‑100 (MQ=255)
                                                                                       tGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCtt              >  1:2337144/1‑100 (MQ=255)
                                                                                        gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc             >  1:698459/1‑100 (MQ=255)
                                                                                        gCTCTCCCCTTCAGGAGAGAGAACTGGACTTTGCTATTAACCCGGAGGCTCATATACCAGCGGCGTCAACTACGGACCGCCCAACGTTACTc                     >  2:3032731/1‑92 (MQ=255)
                                                                                          tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat           <  1:3567557/100‑1 (MQ=255)
                                                                                            tcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatca         <  1:3774464/100‑1 (MQ=255)
                                                                                             ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat        >  1:1160525/1‑100 (MQ=255)
                                                                                             ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat        <  1:2887215/100‑1 (MQ=255)
                                                                                              ctcTTCAGGAGAGCGCGCGGGACAGTGGTATTAACCCTTAGGCTTATTGCCCACGAACTTGGATTACGTCTTGAGCAACGTTACTCTGCGCttactc          >  2:3827808/3‑93 (MQ=255)
                                                                                              cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa       >  1:2634033/1‑100 (MQ=255)
                                                                                                cTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACt     <  1:1722193/100‑1 (MQ=255)
                                                                                                   cAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCa  <  1:1201522/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCA  >  NC_000913/3471263‑3471461

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: