Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,479,939:1 +G 8.8% coding (463/993 nt) prmA → ribosomal protein L11 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,479,9391.G 8.8% 253.8 / 16.3 68coding (463/993 nt)prmAribosomal protein L11 methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (30/32);  new base G (3/3);  total (33/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 = 8.33e-01

GTGTTCGCCCACCTGCGCCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTGGACTTTGG  >  NC_002947/5479795‑5480084
                                                                                                                                                 |                                                                                                                                                  
gtgtTCGCCCACCTGCGCCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCa                                                                                                                                                <  6:108081/149‑1 (MQ=255)
   ttCGCCCACCTGCGCCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGGTGGAGGACCAGGGCGGGGAACGGAGGTGGATGGACCACTTCCAGCCGATGCGCTTCGGCCGGCGCCTGTGGATCGTGCCCAGCTGGGACGa                                                                                                                                             >  1:83935/1‑149 (MQ=255)
      gCCCACCTGCGCCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGCATCGTGCCAGCTGGGCACGAtcc                                                                                                                                          <  4:16238/149‑3 (MQ=255)
               cgcCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGaaaa                                                                                                                                 <  4:265147/149‑1 (MQ=255)
               cgcCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGaaaa                                                                                                                                 <  3:106252/149‑1 (MQ=255)
                          cGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTg                                                                                                                      >  2:172259/1‑148 (MQ=255)
                            gcgcCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAAc                                                                                                                   <  3:324189/149‑1 (MQ=255)
                               ccGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACctg                                                                                                                >  6:178932/1‑149 (MQ=255)
                                cGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACctgc                                                                                                               >  4:305032/1‑149 (MQ=255)
                                cGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACctgc                                                                                                               <  5:111303/149‑1 (MQ=255)
                                       cccGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTgctgct                                                                                                           <  4:199790/146‑1 (MQ=255)
                                       cccGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTgctgct                                                                                                           >  3:199790/1‑146 (MQ=255)
                                         cGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGAcc                                                                                                      <  7:36921/149‑1 (MQ=255)
                                              aCCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGAt                                                                                                                                            <  4:238439/107‑1 (MQ=255)
                                              aCCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGAt                                                                                                                                            >  3:238439/1‑107 (MQ=255)
                                              aCCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGggc                                                                                                 >  4:152913/1‑149 (MQ=255)
                                                cAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGggcct                                                                                               >  1:17595/1‑149 (MQ=255)
                                                cAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGggcct                                                                                               >  5:42244/1‑149 (MQ=255)
                                                 aGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGggcctg                                                                                              <  5:178932/149‑1 (MQ=255)
                                                       aGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCttt                                                                                        <  6:186041/149‑1 (MQ=255)
                                                        ggTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCACCGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGTCCTGGACTTTg                                                                                       >  6:326491/1‑149 (MQ=255)
                                                         gTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTg                                                                                       <  2:359324/148‑1 (MQ=255)
                                                          tGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTg                                                                                       <  4:92550/147‑1 (MQ=255)
                                                            aTCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCa                                                                                    <  6:42244/148‑1 (MQ=255)
                                                              cGAGGACCAGGACTGGGAACGCAGGTGTATGGACAACTTCCAGCCGATGCGGTTCGGCCGCCGCCTGTGGATCGTGCCGAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGggcctg                                                                                              <  2:252742/136‑1 (MQ=255)
                                                              cGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCAGGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCg                                                                                 >  4:124948/1‑149 (MQ=255)
                                                              cGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAACCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGggcctg                                                                                              >  1:252742/1‑136 (MQ=255)
                                                                aGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGACGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCgg                                                                                >  5:219733/1‑148 (MQ=255)
                                                                        gACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAATGAc                                                                                                                            <  3:272975/96‑1 (MQ=255)
                                                                        gACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGAc                                                                                                                            >  4:272975/1‑96 (MQ=255)
                                                                           tGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCccac                                                                    >  2:206353/1‑149 (MQ=255)
                                                                            gggAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCccacc                                                                   >  3:216560/1‑149 (MQ=255)
                                                                             ggAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCcgacca                                                                  >  7:283828/1‑149 (MQ=255)
                                                                                 cGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGc                                                              <  4:216560/149‑1 (MQ=255)
                                                                                    aGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCt                                                           <  1:172259/149‑1 (MQ=255)
                                                                                    aGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCt                                                           >  8:116158/1‑149 (MQ=255)
                                                                                        ggatggaCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGc                                                       >  4:196221/1‑149 (MQ=255)
                                                                                         gatggaCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGTTGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGcc                                                      >  4:135515/1‑149 (MQ=255)
                                                                                         gatggaCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGTTGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGcc                                                      >  4:130367/1‑149 (MQ=255)
                                                                                         gatggaCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGcc                                                      >  3:287147/1‑149 (MQ=255)
                                                                                            ggaCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTg                                                    >  1:117625/1‑148 (MQ=255)
                                                                                              aCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCGACCACGGCCCTGTGCCTGGaa                                                 <  8:283828/149‑1 (MQ=255)
                                                                                              aCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGaa                                                 <  3:124948/149‑1 (MQ=255)
                                                                                              aCAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGaa                                                 >  8:279489/1‑149 (MQ=255)
                                                                                                aaCTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTTGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATg                                               >  4:256995/1‑149 (MQ=255)
                                                                                                   ttCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCt                                            >  8:115860/1‑149 (MQ=255)
                                                                                                     ccAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCt                                            <  3:305032/147‑1 (MQ=255)
                                                                                                             tGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGTTGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCcag                                  <  3:135515/149‑1 (MQ=255)
                                                                                                             tGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGTTGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCcag                                  <  3:130367/149‑1 (MQ=255)
                                                                                                                      gccgccgccTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTTAACCTGCTTCTCGACCCGGGCATGGCCTTTGGCACCGGCACTCACACCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGaa                         <  5:326491/149‑1 (MQ=255)
                                                                                                                        cgccgccTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAgg                       <  4:347231/149‑1 (MQ=255)
                                                                                                                         gccgccTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAgg                       >  4:81325/1‑148 (MQ=255)
                                                                                                                           cgccTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCAc                    <  3:152913/149‑1 (MQ=255)
                                                                                                                               tgtgGATCGTTCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCAGCCAg                <  3:256995/149‑1 (MQ=255)
                                                                                                                               tgtgGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAg                >  2:154243/1‑149 (MQ=255)
                                                                                                                               ggggcaTCGTGCCAGCTGGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAg                        >  2:70339/6‑141 (MQ=255)
                                                                                                                               ggggcaTCGTGCCAGCTGGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAg                        <  1:70339/136‑1 (MQ=255)
                                                                                                                                 tgGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAgg               <  2:83935/148‑1 (MQ=255)
                                                                                                                                 tgGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGt              <  7:116158/149‑1 (MQ=255)
                                                                                                                                 ggcaTCGTGCCAGCTGGGCACGATCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTg             >  7:184526/4‑149 (MQ=255)
                                                                                                                                   gATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGc            <  7:115860/149‑1 (MQ=255)
                                                                                                                                    aTCGTGCCAGCTGGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCt           >  1:82373/1‑149 (MQ=255)
                                                                                                                                     tCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTg          <  3:196221/149‑1 (MQ=255)
                                                                                                                                        tGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTGGAc       >  2:89137/1‑149 (MQ=255)
                                                                                                                                        gtgccAGCTGGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTgg         <  2:82373/145‑1 (MQ=255)
                                                                                                                                           ccAGCTGGCACGATGCCCCGGAAAAGGAGCGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGTACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTGGACtt     >  1:204829/1‑149 (MQ=255)
                                                                                                                                             aGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTGGACTTTgg  <  6:219733/149‑1 (MQ=255)
                                                                                                                                             aGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCGCGACCCGGGCCTGGCCTTTGGTACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTGGACTTTgg  <  2:204829/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTGTTCGCCCACCTGCGCCTGCTGACCGGCGCCGAGCTGCCCGAGCACCAGGCCGAGGTGATCGAGGACCAGGACTGGGAACGCAGCTGGATGGACAACTTCCAGCCGATGCGCTTCGGCCGCCGCCTGTGGATCGTGCCCAGCTGGCACGATGCCCCGGAAAAGGACGCCGTGAACCTGCTGCTCGACCCGGGCCTGGCCTTTGGCACCGGCACTCACCCCACCACGGCCCTGTGCCTGGAATGGCTCGACGGCCAGCAGCTTGAAGGCACCCAGGTGCTGGACTTTGG  >  NC_002947/5479795‑5480084

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: