Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 20.4% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG20.4% 495.2 / 124.0 186intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (80/68);  new base G (36/2);  total (116/70)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

Only 100 of 191 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                         
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                           <  1:3060072/100‑1 (MQ=255)
 tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                          >  2:1637750/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  2:340215/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  2:828539/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                       >  2:220179/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                       >  1:2166964/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      >  2:2122519/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      <  1:2812679/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgn                                                                                                    <  1:1959230/100‑2 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                    <  1:608578/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCATCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGca                                                                                                  <  1:678177/100‑2 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  >  2:2892332/1‑100 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                  <  1:2715559/100‑2 (MQ=37)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                            <  1:848951/100‑1 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                            >  1:1615012/1‑100 (MQ=255)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                           >  2:2623264/1‑100 (MQ=255)
                cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                          <  2:1385507/100‑1 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                         <  2:1790019/100‑1 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                        >  1:1799631/1‑100 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCga                                                                                      >  1:2409739/1‑99 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                      <  1:2244110/100‑1 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                      <  2:90682/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                     >  2:481912/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                   >  2:289802/1‑100 (MQ=255)
                        tctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                                  <  2:522781/97‑1 (MQ=255)
                          agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                <  2:1028483/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  2:1515497/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  2:172828/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  2:925055/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  2:2780876/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  2:1570778/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAACTTCa                                                                             >  1:2620095/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTCTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                           <  1:2149943/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                          >  2:2686805/1‑100 (MQ=255)
                                       ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                   <  1:1937546/100‑1 (MQ=255)
                                       ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                   <  2:2508735/100‑1 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                  <  2:343257/100‑1 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑t                                                        <  1:1757849/100‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑tt                                                       <  1:95882/100‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑tt                                                       <  1:172826/100‑1 (MQ=255)
                                               aCAGATTGTTGACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTc                                                      >  1:2211719/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCg                                                     <  2:2982137/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCg                                                     <  2:13215/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGt                                                    <  2:1840849/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGt                                                    <  1:2031732/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGt                                                    <  1:11977/100‑1 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGa                                            >  1:632603/1‑100 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGa                                            >  1:1319955/1‑100 (MQ=255)
                                                          aCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGAt                                           >  1:2478335/1‑100 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCg                                      >  2:1929379/1‑100 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCg                                      >  1:2881728/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCgg                                     >  1:2377124/1‑100 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAg                                   <  1:565655/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAAAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                                 <  2:2453636/100‑1 (MQ=255)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCAt                              >  1:2101650/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                            >  2:538318/1‑100 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAg                           >  2:1198353/5‑100 (MQ=18)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:608637/100‑1 (MQ=255)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1214858/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAA‑‑‑CTTGTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:1037181/3‑100 (MQ=21)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATA‑‑T‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                        >  2:2407039/3‑99 (MQ=21)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  1:113654/1‑100 (MQ=255)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  2:2311869/1‑100 (MQ=255)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AAT‑‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  2:1753466/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2257811/2‑100 (MQ=17)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  1:459257/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGG‑‑C‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1748693/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCag                       >  2:1647062/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑C‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGCag                       >  2:1478872/1‑99 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CACGTTGGCCTGATAAGCGTAGTGCATCAGACa                        >  2:1937509/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:2188726/9‑98 (MQ=17)
                                                                               ccggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:126444/9‑100 (MQ=18)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     >  2:404064/1‑100 (MQ=255)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑C‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1051034/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCT‑‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGACagctt                    >  1:944523/7‑96 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCT‑‑C‑‑‑AACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:2146338/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCAT‑‑C‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                    >  2:1588636/7‑99 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑C‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGCagctt                    >  2:1958129/7‑96 (MQ=21)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                   >  1:817773/6‑100 (MQ=17)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 >  1:2781560/1‑100 (MQ=255)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 >  2:1010133/1‑100 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  1:1383722/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:568545/1‑100 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATCAGCGGAGCGCATCAGGCATTtttgcctt              >  1:1460948/1‑100 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              >  1:1549970/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  1:1096824/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1527646/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:183407/100‑1 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  2:2681053/100‑1 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  2:2510357/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:1121902/1‑100 (MQ=255)
                                                                                            naaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:2380731/2‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:507124/100‑1 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  2:2435044/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:2648685/100‑1 (MQ=255)
                                                                                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGGGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  2:2577710/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  2:787981/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  2:2021464/1‑100 (MQ=255)
                                                                                                   ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  >  2:1544962/1‑100 (MQ=255)
                                                                                                   ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:787976/100‑1 (MQ=255)
                                                                                                   |                                                                                                         
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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