Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,159 C→T 7.7% S126L (TCG→TTG)  PP_3045 → ATP‑dependent protease ClpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,430,1590CT 7.7% 187.4 / 4.5 65S126L (TCG→TTG) PP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (2/3);  total (30/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 = 1.16e-01

ACAAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGA  >  NC_002947/3430014‑3430306
                                                                                                                                                 |                                                                                                                                                   
acaAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGaa                                                                                                                                                  >  5:224803/1‑149 (MQ=255)
   aaCGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGc                                                                                                                                               >  8:326415/1‑149 (MQ=255)
    aCGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGcc                                                                                                                                              <  7:344324/149‑1 (MQ=255)
       gAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATg                                                                                                                                           <  3:187207/149‑1 (MQ=255)
        aGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATga                                                                                                                                          <  7:337062/149‑2 (MQ=255)
        aGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATga                                                                                                                                          <  5:181359/149‑2 (MQ=255)
           tGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCtga                                                                                                                                       <  4:181848/149‑1 (MQ=255)
               cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTgatgat                                                                                                                                   >  2:346095/1‑149 (MQ=255)
               cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTgatgat                                                                                                                                   >  4:227254/1‑149 (MQ=255)
                tGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATc                                                                                                                                  >  3:313582/1‑149 (MQ=255)
                 gTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATc                                                                                                                                  >  1:23868/1‑148 (MQ=255)
                  ttCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATcc                                                                                                                                 <  5:67044/148‑1 (MQ=255)
                  ttCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATc                                                                                                                                  <  1:346095/147‑1 (MQ=255)
                     gAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATaa                                                                                                                             >  3:177327/1‑149 (MQ=255)
                         ggCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCt                                                                                                                         >  2:183042/1‑149 (MQ=255)
                          gCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCtt                                                                                                                        <  4:284119/149‑1 (MQ=255)
                          gCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCt                                                                                                                         >  5:363608/1‑148 (MQ=255)
                                cATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCt                                                                                                                  >  1:305855/1‑149 (MQ=255)
                                 aTCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTg                                                                                                                 >  2:1085/1‑149 (MQ=255)
                                 aTCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCACAACCCTTACACCTg                                                                                                                 <  7:318858/149‑1 (MQ=255)
                                  tCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCTTGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTgg                                                                                                                <  1:76169/149‑1 (MQ=255)
                                  tCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTgg                                                                                                                <  3:300099/149‑1 (MQ=255)
                                    cacaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTggg                                                                                                               <  4:174072/148‑1 (MQ=255)
                                     acaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATc                                                                                                                                                     <  7:143838/109‑1 (MQ=255)
                                     acaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATc                                                                                                                                                     >  8:143838/1‑109 (MQ=255)
                                     acaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCAGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGtt                                                                                                             >  6:335297/1‑149 (MQ=255)
                                        aCGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGAGGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGGGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGc                                                                                                          <  3:227254/149‑1 (MQ=255)
                                         cgcgCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCg                                                                                                         <  8:133299/149‑1 (MQ=255)
                                           cgcTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGt                                                                                                       <  8:4247/149‑1 (MQ=255)
                                              tGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGAt                                                                                                    >  6:297829/1‑149 (MQ=255)
                                              tGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGAt                                                                                                    >  5:99115/1‑149 (MQ=255)
                                                   cGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGa                                                                                               <  7:326415/149‑1 (MQ=255)
                                                   cGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGa                                                                                               <  4:313582/149‑1 (MQ=255)
                                                    gCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGACATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAg                                                                                              >  3:134905/1‑149 (MQ=255)
                                                        gggTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCTAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGAcc                                                                                          >  1:236533/1‑149 (MQ=255)
                                                        gggTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGa                                                                                            <  6:99115/147‑1 (MQ=255)
                                                                   gCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGt                                                                                <  2:305855/148‑1 (MQ=255)
                                                                    ctccGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGTCGGTAATGACGAGGACCTTCGTCGGGTg                                                                               <  4:177327/146‑1 (MQ=255)
                                                                    cACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTgg                                                                              >  1:23730/1‑149 (MQ=255)
                                                                     aCCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGAc                                                                                           <  2:104049/135‑1 (MQ=255)
                                                                     aCCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGAc                                                                                           >  1:104049/1‑135 (MQ=255)
                                                                              gTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCt                                                                    >  3:159109/1‑149 (MQ=255)
                                                                              gTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCt                                                                    <  3:23514/149‑1 (MQ=255)
                                                                                  aTGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGAt                                                                >  5:299647/1‑149 (MQ=255)
                                                                                               aaGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAgg                                                   <  1:148266/149‑1 (MQ=255)
                                                                                                 gTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCa                                                 <  1:183042/149‑1 (MQ=255)
                                                                                                       gcagCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGa                                                                                                     <  7:119967/91‑1 (MQ=255)
                                                                                                       gcagCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGa                                                                                                     >  8:119967/1‑91 (MQ=255)
                                                                                                       gcagCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATc                                           <  5:238771/149‑1 (MQ=255)
                                                                                                       gcagCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATc                                           >  1:271991/1‑149 (MQ=255)
                                                                                                       gcagCGTGGCAGTCTGCGGCGCTCACCGTGAGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATc                                           >  1:29769/1‑149 (MQ=255)
                                                                                                            gTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCTTCGGGTTGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGc                                      <  2:23868/149‑1 (MQ=255)
                                                                                                                        ggcgCTCACCGTGTGGTCATGGCATTGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCaa                          <  5:297829/149‑1 (MQ=255)
                                                                                                                        ggcgCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGcc                            >  3:247160/1‑147 (MQ=255)
                                                                                                                        ggcgCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCaa                          <  4:247160/149‑1 (MQ=255)
                                                                                                                        ggcgCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCaa                          <  5:335297/149‑1 (MQ=255)
                                                                                                                           gcTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGc                       <  6:363608/149‑1 (MQ=255)
                                                                                                                            cTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGc                                                                     <  2:28658/102‑1 (MQ=255)
                                                                                                                            cTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGc                                                                     >  1:28658/1‑102 (MQ=255)
                                                                                                                             tCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCta                     >  3:286250/1‑148 (MQ=255)
                                                                                                                                       gTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGa           >  3:301550/1‑149 (MQ=255)
                                                                                                                                       gTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGa           >  2:243611/1‑149 (MQ=255)
                                                                                                                                       agcATGGCATTCGATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGa           >  1:59841/3‑149 (MQ=255)
                                                                                                                                         cATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGAt          >  7:137015/1‑148 (MQ=255)
                                                                                                                                             gCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGAcgc     <  2:23730/149‑1 (MQ=255)
                                                                                                                                              cATCGAATGCCATGCTGATGATCCATAACCCTTACTCCTGGGTTGGCGTTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCATTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGAcgcg    <  4:134905/149‑1 (MQ=255)
                                                                                                                                                tCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGa  <  1:1085/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACAAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGA  >  NC_002947/3430014‑3430306

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