Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,850,295 G→A 6.2% R279H (CGC→CAC)  rlmD → 23S rRNA (uracil(1939)‑C(5))‑methyltransferase RlmD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,850,2950GA 6.2% 190.3 / 2.1 64R279H (CGC→CAC) rlmD23S rRNA (uracil(1939)‑C(5))‑methyltransferase RlmD
Reads supporting (aligned to +/- strand):  ref base G (36/24);  new base A (2/2);  total (38/26)
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.45e-01

CCGGCTGAGGACCAGGCCAGGTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGCCCGCCAGGCG  >  NC_002947/1850147‑1850443
                                                                                                                                                    |                                                                                                                                                     
ccGGCTGAGGACCAGGCCAGGTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGgcg                                                                                                                                                       >  7:6119/1‑149 (MQ=255)
         gACCAGGCCAGGTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCAGGGCGa                                                                                                                                              >  7:176555/1‑149 (MQ=255)
            caggccaggTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAt‑tt                                                                                                                                          >  1:37353/1‑149 (MQ=255)
                    gTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAgg                                                                                                                                  >  6:137518/1‑149 (MQ=255)
                         aaGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAAt                                                                                                                             <  8:31695/149‑1 (MQ=255)
                           gCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGc                                                                                                                           >  1:9469/1‑149 (MQ=255)
                           gCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGc                                                                                                                           >  2:328062/1‑149 (MQ=255)
                            cGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGcc                                                                                                                          <  1:328062/149‑1 (MQ=255)
                              ttCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGcc                                                                                                                          >  6:305260/1‑147 (MQ=255)
                                  gCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGt                                                                                                                     >  6:175320/1‑148 (MQ=255)
                                    gAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCTTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCACCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCaa                                                                                                                  >  3:231721/1‑149 (MQ=255)
                                    gAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCTTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCACCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCaa                                                                                                                  >  3:236471/1‑149 (MQ=255)
                                              aCGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCgatg                                                                                                        <  8:352777/149‑1 (MQ=255)
                                               cGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCgatg                                                                                                        >  3:295484/1‑148 (MQ=255)
                                                   cAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCgat                                                                                                         >  7:204501/1‑143 (MQ=255)
                                                   cAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCgat                                                                                                         <  8:204501/143‑1 (MQ=255)
                                                            cTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCttt                                                                                          <  2:37353/149‑1 (MQ=255)
                                                             tGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTg                                                                                         >  3:166706/1‑149 (MQ=255)
                                                             tGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTg                                                                                         >  8:109436/1‑149 (MQ=255)
                                                               catggCGAAGGTGAGGCGGCGCGCGTGGTCGCAGGGGCGCAGTTGGGCTTTGCCGTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCATGCTTTgg                                                                                        <  7:159704/145‑1 (MQ=255)
                                                               cAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTg                                                                                         <  5:137518/147‑1 (MQ=255)
                                                               cAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGc                                                                                       <  8:6119/149‑1 (MQ=255)
                                                                  ggCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTg                                                                                    >  1:207592/1‑149 (MQ=255)
                                                                   gCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTgg                                                                                   <  6:9653/149‑1 (MQ=255)
                                                                    cGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgc                                                                                  <  2:187006/149‑1 (MQ=255)
                                                                     gAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgct                                                                                 >  6:93937/1‑149 (MQ=255)
                                                                         gTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgctcgct                                                                             >  3:37437/1‑149 (MQ=255)
                                                                         gTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgctcgct                                                                             >  6:44850/1‑149 (MQ=255)
                                                                               cGGCTCCCGTGGTCCCAGCTGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTTCAGGTTAATGCCGAGGTGAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACaa                                                                       <  8:176555/149‑1 (MQ=255)
                                                                                 gcgcCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCCCAAGc                                                                     >  2:57722/1‑149 (MQ=255)
                                                                                 gcgcCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGc                                                                     <  5:44850/149‑1 (MQ=255)
                                                                                  cgcCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGcc                                                                    >  8:130870/1‑149 (MQ=255)
                                                                                          gACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGc                                                             <  4:166706/148‑1 (MQ=255)
                                                                                             ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCt                                                                                     <  8:301303/121‑1 (MQ=255)
                                                                                             ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCt                                                                                     >  7:301303/1‑121 (MQ=255)
                                                                                             ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCgtgt                                                         >  1:1568/1‑149 (MQ=255)
                                                                                                   gcgcAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGa                                                   >  8:285135/1‑149 (MQ=255)
                                                                                                    cgcAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGAc                                                  <  5:175320/149‑1 (MQ=255)
                                                                                                       aGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTg                                               <  1:147289/149‑1 (MQ=255)
                                                                                                             gCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGc                                         <  2:9469/149‑1 (MQ=255)
                                                                                                              cTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTg                                          >  6:140293/1‑147 (MQ=255)
                                                                                                              cTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTg                                          <  5:140293/147‑1 (MQ=255)
                                                                                                               tttGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGTGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCgg                                       >  2:101477/1‑149 (MQ=255)
                                                                                                                  gCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCt                                    >  5:183738/1‑149 (MQ=255)
                                                                                                                  gCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCt                                    >  4:343339/1‑149 (MQ=255)
                                                                                                                   cccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTg                                   >  6:106223/1‑149 (MQ=255)
                                                                                                                   cccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGCGCTGGACTTGTTCTGCGGTCTg                                   >  1:107622/1‑149 (MQ=255)
                                                                                                                    ccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTgg                                  <  8:168864/149‑1 (MQ=255)
                                                                                                                       ggCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGc                                >  7:111067/1‑148 (MQ=255)
                                                                                                                          ccccTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATATTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGCGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGACCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAAc                             >  1:66598/1‑149 (MQ=255)
                                                                                                                          ccccTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACt                            >  5:220168/1‑149 (MQ=255)
                                                                                                                           cccTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACtt                           >  6:24578/1‑149 (MQ=255)
                                                                                                                            ccTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACttt                          <  2:207592/149‑1 (MQ=255)
                                                                                                                             cTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTTTACTTGTTCTGCGGTGTGGGCATCTTTg                         <  2:359092/149‑1 (MQ=255)
                                                                                                                                gCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGCGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGcc                       <  2:107622/148‑1 (MQ=255)
                                                                                                                                   gctggagctggCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAAAGTTGcccaga                   <  7:109436/149‑4 (MQ=255)
                                                                                                                                     tggagctggCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCg                 >  3:12394/1‑149 (MQ=255)
                                                                                                                                       gagctggCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCt               >  7:326580/1‑149 (MQ=255)
                                                                                                                                       gagctggCCTGGCACCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAAAGTTGcccagacc                 <  4:236471/147‑6 (MQ=255)
                                                                                                                                       gagctggCCTGGCACCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAAAGTTGcccagacc                 <  4:231721/147‑6 (MQ=255)
                                                                                                                                           tggcctggcGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGcc           >  5:94208/1‑149 (MQ=255)
                                                                                                                                             gcctggcGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGCCCg         >  7:102005/1‑149 (MQ=255)
                                                                                                                                               ctggcGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGCCCGc        >  2:351936/1‑148 (MQ=255)
                                                                                                                                                 ggcGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGc                                                       >  6:193200/1‑99 (MQ=255)
                                                                                                                                                 ggcGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGc                                                       <  5:193200/99‑1 (MQ=255)
                                                                                                                                                    gcCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGCCCGCCAGgcg  >  5:79601/1‑149 (MQ=255)
                                                                                                                                                    gcCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGCCCGCCAGgcg  >  8:74327/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGCTGAGGACCAGGCCAGGTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGAT‑TTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTGCCGCTGGCCCGCCAGGCG  >  NC_002947/1850147‑1850443

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