Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,375,213 T→G 5.3% P76P (CCA→CCC)  dnaJ ← chaperone protein DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,375,2130TG 5.3% 238.7 / 4.3 76P76P (CCA→CCC) dnaJchaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base T (39/33);  new base G (2/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GCTCCAGGGTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATC  >  NC_002947/5375065‑5375354
                                                                                                                                                    |                                                                                                                                             
gCTCCAGGGTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCtt                                                                                                                                               <  6:355315/149‑1 (MQ=255)
gCTCCAGGGTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCtt                                                                                                                                               <  1:369755/149‑1 (MQ=255)
   ccAGGGTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTggg                                                                                                                                            >  2:261639/1‑149 (MQ=255)
    cAGGGTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACGGAAGAAGTCGCTGAACACATGGCCGAAGAGCTCGGAGAAAATGGCACACCCAACGCCGGCACCCCCGCCCCCCAGCCTTgg                                                                                                                                             >  7:256941/1‑147 (MQ=255)
    cAGGGTGGAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGt                                                                                                                                           <  3:224092/149‑1 (MQ=255)
       ggTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGa                                                                                                                                        >  4:334398/1‑149 (MQ=255)
        gtgtAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGa                                                                                                                                        >  1:59973/1‑148 (MQ=255)
          gtAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCCCCCAGGCTTGGGTCCacac                                                                                                                                     >  8:286470/1‑149 (MQ=255)
             gcgcAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCt                                                                                                                                   <  3:334398/148‑1 (MQ=255)
             gcgcAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTg                                                                                                                                  >  3:129522/1‑149 (MQ=255)
              cgcAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGc                                                                                                                                 <  7:229819/149‑1 (MQ=255)
                 ccggcGCTGACGCGCTTGGTACCGCCTCGACCACCACCGCGACCAGCGCCTCCGAACAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATg                                                                                                                              <  7:168626/145‑1 (MQ=255)
                   gTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCAACAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGc                                                                                                                            >  2:223612/1‑149 (MQ=255)
                    tCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAACAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCAGGTGTCGACCCAAGCATgggt                                                                                                                           <  1:310239/149‑3 (MQ=255)
                     cGCTGCCGCGCTGGGCACCGCCTCGACCACCACAGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCa                                                                                                                          >  3:110384/1‑149 (MQ=255)
                        tGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                       >  7:229282/1‑149 (MQ=255)
                         gCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACt                                                                                                                      >  2:296544/1‑149 (MQ=255)
                          ccGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACt                                                                                                                      <  7:286470/148‑1 (MQ=255)
                                 gggCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAg                                                                                                              >  4:179401/1‑149 (MQ=255)
                                  ggCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAgg                                                                                                             <  2:229740/149‑1 (MQ=255)
                                   gCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGc                                                                                                            >  3:181435/1‑149 (MQ=255)
                                    cACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGcc                                                                                                           >  3:332045/1‑149 (MQ=255)
                                       cGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCa                                                                                                        <  4:110384/149‑1 (MQ=255)
                                        gCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCAc                                                                                                       <  5:336892/149‑1 (MQ=255)
                                           tCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCt                                                                                                    >  5:85522/1‑149 (MQ=255)
                                             gACCACCACCGCGACCACCGCCACCGAACAAGGCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCTTGGCCATACTGGTCGAAGGCCGCACGCt                                                                                                    <  5:91984/147‑1 (MQ=255)
                                                caccaccGCGACCACCGCCACCGAAGAAGTCGATGAACACATCTCCGAAGATGTCGGAGAAATTGGCACCAAAAAAGCCGGCACCACCGCCAACCATGCTTGGCTCGACACCAGCATGGCCATACTGGTCGAATGCCGCACGCTAGCTc                                                                                               >  6:381005/1‑149 (MQ=255)
                                                 accaccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCAGGTGTCGACCCAAGCATgggtg                                                                                                                          >  2:186959/1‑118 (MQ=255)
                                                  ccaccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTcgc                                                                                             >  2:338848/1‑149 (MQ=255)
                                                     ccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTc                                                                                          >  7:321443/1‑149 (MQ=255)
                                                     ccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTc                                                                                          <  7:137596/149‑1 (MQ=255)
                                                        cgACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                       <  3:240696/149‑1 (MQ=255)
                                                        cgACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                       <  1:296544/149‑1 (MQ=255)
                                                        cgACCACCGCCACCGAAGAAGTCGCTGAACACACCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                       >  8:220305/1‑149 (MQ=255)
                                                          accaccGACACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCTCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACa                                                                                     <  5:381005/149‑1 (MQ=255)
                                                                   accgAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCg                                                                             >  1:155357/1‑148 (MQ=255)
                                                                    cctaagaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTa                                                                           <  8:90565/146‑1 (MQ=255)
                                                                    ccgAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTa                                                                           <  8:321443/149‑1 (MQ=255)
                                                                     cgAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAg                                                                          <  4:129522/149‑1 (MQ=255)
                                                                      gaagaaGTCGCTGAACACATCGCCGAGGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAgg                                                                         >  1:202952/1‑149 (MQ=255)
                                                                         gaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGTACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCt                                                                                         >  3:225728/1‑130 (MQ=255)
                                                                         gaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGTACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCt                                                                                         <  4:225728/130‑1 (MQ=255)
                                                                         gaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCt                                                                      <  2:59973/149‑1 (MQ=255)
                                                                          aagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCCCCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTc                                                                     >  4:157154/1‑149 (MQ=255)
                                                                             tCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGtt                                                                  >  2:174193/1‑149 (MQ=255)
                                                                                 tGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGcc                                                              >  1:277223/1‑149 (MQ=255)
                                                                                       caTCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTg                                                        <  1:261639/149‑1 (MQ=255)
                                                                                          cGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAAc                                                     <  7:220305/149‑1 (MQ=255)
                                                                                           gCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAAc                                                     <  2:202952/148‑1 (MQ=255)
                                                                                           gCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAAc                                                     <  2:155357/148‑1 (MQ=255)
                                                                                            ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                       >  2:375573/1‑81 (MQ=255)
                                                                                            ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                       <  1:375573/81‑1 (MQ=255)
                                                                                            ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACtg                                                   >  8:220828/1‑148 (MQ=255)
                                                                                               aaGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGAc                                                                                      >  6:382732/1‑111 (MQ=255)
                                                                                               aaGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGAc                                                                                      <  5:382732/111‑1 (MQ=255)
                                                                                                  aTGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTc                                             <  6:85522/149‑1 (MQ=255)
                                                                                                   tGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTc                                             <  1:223612/148‑1 (MQ=255)
                                                                                                    gTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTcc                                            >  8:141004/1‑148 (MQ=255)
                                                                                                       ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGt                                                 >  4:141548/1‑140 (MQ=255)
                                                                                                       ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGt                                                 <  3:141548/140‑1 (MQ=255)
                                                                                                       ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGAtt                                        <  6:111126/149‑1 (MQ=255)
                                                                                                       ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGAtt                                        >  1:87078/1‑149 (MQ=255)
                                                                                                       ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTC‑AAGGCCGCACGCATGCGCGCGGAGGACAGCACTTCGTAAGCGTAGTTGGCCGCCGTGAACTTGTCGTACGATTc                                       >  8:347502/1‑149 (MQ=255)
                                                                                                        gagaAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTc                                       >  4:52273/1‑149 (MQ=255)
                                                                                                         agaAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCg                                                                    >  4:340632/1‑119 (MQ=255)
                                                                                                         agaAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCg                                                                    <  3:340632/119‑1 (MQ=255)
                                                                                                            aaTTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTa                                   <  4:181435/149‑1 (MQ=255)
                                                                                                              ttGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATc                                 >  7:219947/1‑149 (MQ=255)
                                                                                                                   accaccAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGcctc                                                            <  6:178229/117‑1 (MQ=255)
                                                                                                                   accaccAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGcctc                                                            >  5:178229/1‑117 (MQ=255)
                                                                                                                        caAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTg                       >  6:158813/1‑149 (MQ=255)
                                                                                                                        caAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACAGGTCGAAGGCCGCACGCTTGCTGGCGTGTGACAGCACTTCGTAGGCCTCGTTGGCCTCCGTGAACTTGGCTTCCGATTCTTTTTAGCCCGGGTTg                       >  1:273516/1‑149 (MQ=255)
                                                                                                                          aaGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCg                     >  4:215771/1‑149 (MQ=255)
                                                                                                                                       gtatggccATGCAGGTGTCGACCCAAGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTATGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGggtg            <  3:48515/139‑1 (MQ=255)
                                                                                                                                       gtatggccATGCAGGTGTCGACCCAAGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTATGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGggtg            >  4:48515/7‑145 (MQ=255)
                                                                                                                                           aCCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCa    <  7:141004/149‑1 (MQ=255)
                                                                                                                                             ccATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATc  >  3:96385/1‑149 (MQ=255)
                                                                                                                                             ccATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATc  <  1:174193/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCTCCAGGGTGTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATC  >  NC_002947/5375065‑5375354

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: ‑