Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,375,225 T→A 6.0% A72A (GCA→GCT)  dnaJ ← chaperone protein DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,375,2250TA 6.0% 230.2 / 2.7 77A72A (GCA→GCT) dnaJchaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (3/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAACGGCGGT  >  NC_002947/5375078‑5375367
                                                                                                                                                   |                                                                                                                                              
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                                                                                                    <  1:296544/149‑1 (MQ=255)
                                           cgACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                                    <  3:240696/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                                                                                                                                    <  1:375573/81‑1 (MQ=255)
                                                                               ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                                    >  2:375573/1‑81 (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)
                                                                                                                                         ggTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcgg      <  3:179401/149‑1 (MQ=255)
                                                                                                                                          gTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcggc     <  3:157154/149‑1 (MQ=255)
                                                                                                                                           tCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcggcg    >  7:15713/1‑149 (MQ=255)
                                                                                                                                           tCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCACCAAAcggcg    <  8:229282/149‑1 (MQ=255)
                                                                                                                                             gACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcggcgg   >  2:17325/1‑148 (MQ=255)
                                                                                                                                              acacCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAACGGCGGt  >  4:168134/1‑148 (MQ=255)
                                                                                                                                                  cTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTc                                                          >  2:324463/1‑88 (MQ=255)
                                                                                                                                                  cTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTc                                                          <  1:324463/88‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAACGGCGGT  >  NC_002947/5375078‑5375367

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