Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,293 C→T 6.5% intergenic (‑588/+323) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,2930CT 6.5% 286.5 / 6.8 94intergenic (‑588/+323)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base C (45/40);  new base T (3/3);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

GCCAGCCATAAGTAACGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTT  >  NC_000913/3708148‑3708441
                                                                                                                                                 |                                                                                                                                                    
gccagccaTAAGTAACGAAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGcgcc                                                                                                                                                   <  8:352044/149‑1 (MQ=255)
     ccaTAAGTAACGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATc                                                                                                                                              >  3:12235/1‑149 (MQ=255)
        tAAGTAACGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAGGGCGGTCGCGCCGCATCCGa                                                                                                                                           >  8:213888/1‑149 (MQ=255)
           gTAACGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATTAGGAGTTCGAGAAGCATCCGACaa                                                                                                                                        >  7:12453/1‑149 (MQ=255)
            tAACGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAAt                                                                                                                                       <  7:280369/149‑1 (MQ=255)
              aCGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATgg                                                                                                                                     <  7:213888/149‑1 (MQ=255)
                gCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGcc                                                                                                                                   >  4:214231/1‑149 (MQ=255)
                      aaaGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGcacc                                                                                                                             >  7:255186/1‑148 (MQ=255)
                      aaaGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGcaca                                                                                                                             >  7:249011/1‑149 (MQ=255)
                      aaaGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGcaca                                                                                                                             >  3:354329/1‑149 (MQ=255)
                         gATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACAttt                                                                                                                          <  7:10773/149‑1 (MQ=255)
                           tttAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGGGTGTCGTGTCGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGc                                                                                                                        <  5:128902/149‑1 (MQ=255)
                             tAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCt                                                                                                                      >  1:116957/1‑149 (MQ=255)
                             tAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCt                                                                                                                      >  7:292587/1‑149 (MQ=255)
                              aaGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCt                                                                                                                      <  1:36659/148‑1 (MQ=255)
                               aGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGATCGAATAAGTCGCTAGCGCATAATACAAAAATGGCCGGAACATTTGCCTGa                                                                                                                    >  8:169297/1‑149 (MQ=255)
                                gCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGAt                                                                                                                   <  7:395653/149‑1 (MQ=255)
                                gCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGAt                                                                                                                   >  4:226101/1‑149 (MQ=255)
                                 cAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATg                                                                                                                  >  7:403248/1‑149 (MQ=255)
                                 cAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATg                                                                                                                  <  5:145033/149‑1 (MQ=255)
                                 cAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCAGTTACTGGATATTGAGATATGGGGCTCAAATTTGTAAGGGAAATTTGTTTAGTTTAGGTAGAATACGGCGGTCGAGCAGCATCCGACAATGGCCGTCACATTTGCCTGATg                                                                                                                  >  7:288472/1‑149 (MQ=255)
                                     tataAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGAc                                                                                                              <  7:231790/149‑1 (MQ=255)
                                           aaaaGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGAc                                                                                                        >  1:414757/1‑149 (MQ=255)
                                           aaaaGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGAc                                                                                                        >  5:103901/1‑149 (MQ=255)
                                             aaGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGAcgc                                                                                                      <  7:245395/149‑1 (MQ=255)
                                               gACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGt                                                                                                    >  8:56815/1‑149 (MQ=255)
                                               gACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGt                                                                                                    <  8:400507/149‑1 (MQ=255)
                                                ccAATGGTGTCTGTCACAAATCCGTTAATGGATAGTGATATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTc                                                                                                   <  8:255186/148‑1 (MQ=255)
                                                aCAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTc                                                                                                   <  2:116957/149‑1 (MQ=255)
                                                aCAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGATCTAATAAGTCGCTCTATCATAATAAAACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTc                                                                                                   >  2:135684/1‑149 (MQ=255)
                                                 cAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTc                                                                                                   <  8:249011/148‑1 (MQ=255)
                                                 cAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCt                                                                                                  >  1:373439/1‑149 (MQ=255)
                                                  aaTGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCt                                                                                                  >  1:304668/1‑148 (MQ=255)
                                                    tGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTAt                                                                                               <  4:76082/149‑1 (MQ=255)
                                                         tCTGTCACTTATCGGTTAATGGATAGTGAGATATGGGGCGCAATTTTGGGAGGGATATGTGTGTCGTGTAGGTCGAATATGGCGGTCGCGCCGATTCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                           <  6:362994/148‑1 (MQ=255)
                                                             tCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                       <  4:12235/148‑1 (MQ=255)
                                                             tCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                      <  8:211958/149‑1 (MQ=255)
                                                             tCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                      >  8:413183/1‑149 (MQ=255)
                                                               acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCTCCGCATCCGAAAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACgt                                                                                    >  2:56615/1‑149 (MQ=255)
                                                                 aaaTCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACgtc                                                                                   >  2:299248/1‑148 (MQ=255)
                                                                      cGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACgtcgagtc                                                                              <  8:351622/148‑1 (MQ=255)
                                                                        ttAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGaa                                                                           >  5:238622/1‑149 (MQ=255)
                                                                            tGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTc                                                                       >  1:95712/1‑149 (MQ=255)
                                                                             ggATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCt                                                                                        <  2:229902/131‑1 (MQ=255)
                                                                             ggATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCt                                                                                        >  1:229902/1‑131 (MQ=255)
                                                                              gATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTccc                                                                     <  2:304668/149‑1 (MQ=255)
                                                                              gATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTcc                                                                      <  3:297029/148‑1 (MQ=255)
                                                                                       atatGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCGCCCTGCCCCCg                                                             >  3:281876/1‑148 (MQ=255)
                                                                                        tatGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCg                                                             >  2:101457/1‑147 (MQ=255)
                                                                                        tatGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGa                                                           <  8:403248/149‑1 (MQ=255)
                                                                                        tatGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGa                                                           <  6:103901/149‑1 (MQ=255)
                                                                                          tGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAg                                                          <  3:214231/148‑1 (MQ=255)
                                                                                            gggCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGtt                                                       >  5:185386/1‑149 (MQ=255)
                                                                                              gcgcAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGTCGGACGCTCCGAATAAGAAAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTtct                                                     >  8:141759/1‑149 (MQ=38)
                                                                                              gcgcAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATTGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTtct                                                     >  7:245883/1‑149 (MQ=255)
                                                                                              gcgcAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGACCCTCCCTGCCCCCGGAGGTtct                                                     >  7:157245/1‑149 (MQ=255)
                                                                                              gcgcAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTtct                                                     >  5:69083/1‑149 (MQ=255)
                                                                                                 cAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGaa                                                                           <  5:9205/124‑1 (MQ=255)
                                                                                                 cAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGaa                                                                           >  6:9205/1‑124 (MQ=255)
                                                                                                   aaTTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATcc                                                >  8:61185/1‑149 (MQ=255)
                                                                                                    aTTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCt                                               >  7:68840/1‑149 (MQ=255)
                                                                                                     tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGTGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCt                                               >  8:81716/1‑148 (MQ=255)
                                                                                                     tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGTGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCt                                               >  6:382039/1‑148 (MQ=255)
                                                                                                     tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGTGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCt                                               >  6:388077/1‑148 (MQ=255)
                                                                                                       tGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCa                                             <  3:226101/148‑1 (MQ=255)
                                                                                                            gggAAATGTGTGGCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCACGAGCt                                       <  7:262765/149‑1 (MQ=255)
                                                                                                              gAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGTGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGc                                     <  5:388077/149‑1 (MQ=255)
                                                                                                              gAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGTGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGc                                     <  5:382039/149‑1 (MQ=255)
                                                                                                              gAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGTGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGc                                     <  7:81716/149‑1 (MQ=255)
                                                                                                              gAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTAATCCTCATGAGCTGc                                     >  3:265966/1‑149 (MQ=255)
                                                                                                                 aTGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCAt                                            <  8:374456/139‑1 (MQ=255)
                                                                                                                 aTGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCAt                                            >  7:374456/1‑139 (MQ=255)
                                                                                                                 aTGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTg                                      >  1:160823/1‑145 (MQ=255)
                                                                                                                  tgtgtgTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTATTATCAGTCCTACGTCGATTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAt                                 <  8:12453/149‑1 (MQ=255)
                                                                                                                    tgtgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATg                                >  4:220478/1‑148 (MQ=255)
                                                                                                                    tgtgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATg                                >  1:274049/1‑148 (MQ=255)
                                                                                                                    tgtgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGc                               <  7:61185/149‑1 (MQ=255)
                                                                                                                      tgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAAc                                                                          >  5:150195/1‑104 (MQ=255)
                                                                                                                      tgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAAc                                                                          <  6:150195/104‑1 (MQ=255)
                                                                                                                         cGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCaaaa                           >  5:320407/1‑148 (MQ=255)
                                                                                                                          gtgtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCaaaaaa                         <  2:274049/149‑1 (MQ=255)
                                                                                                                          gtgtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCaaaaaa                         <  8:292587/149‑1 (MQ=255)
                                                                                                                          gggtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCaaaaaa                         <  2:414757/147‑1 (MQ=255)
                                                                                                                            gtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGc                       >  4:187817/1‑149 (MQ=255)
                                                                                                                            gtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGc                       >  6:186806/1‑149 (MQ=255)
                                                                                                                             tAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGc                       <  8:68840/148‑1 (MQ=255)
                                                                                                                               ggTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCaaaaaa                         <  5:300805/144‑1 (MQ=255)
                                                                                                                               ggTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCaaaaaa                         >  6:300805/1‑144 (MQ=255)
                                                                                                                                gTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGc                   <  1:101457/149‑1 (MQ=255)
                                                                                                                                gTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGc                   >  1:273989/1‑149 (MQ=255)
                                                                                                                                     aTAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt              >  2:118263/1‑149 (MQ=255)
                                                                                                                                     aTAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt              <  5:186806/149‑1 (MQ=255)
                                                                                                                                        aGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGa            <  6:238622/148‑1 (MQ=255)
                                                                                                                                         ggCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGc          <  1:299248/149‑1 (MQ=255)
                                                                                                                                           cGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAg        >  5:382837/1‑149 (MQ=255)
                                                                                                                                                 cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt  <  2:273989/149‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt  <  3:220478/149‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt  <  1:135684/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCAGCCATAAGTAACGCAAAAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTT  >  NC_000913/3708148‑3708441

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