Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,825,231 T→C 7.0% E81G (GAA→GGA)  trmB ← tRNA (guanosine(46)‑N7)‑methyltransferase TrmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,825,2310TC 7.0% 249.3 / 6.6 86E81G (GAA→GGA) trmBtRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base T (45/35);  new base C (4/2);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGG  >  NC_002947/5825083‑5825378
                                                                                                                                                    |                                                                                                                                                   
gCAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCAttt                                                                                                                                                     >  1:293579/1‑149 (MQ=255)
gCAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCAttt                                                                                                                                                     >  6:48994/1‑149 (MQ=255)
 cAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGGCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTc                                                                                                                                                    <  2:338891/149‑1 (MQ=255)
  aGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTc                                                                                                                                                    >  4:188245/1‑148 (MQ=255)
       ttCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCcagcag                                                                                                                                              <  5:234825/149‑1 (MQ=255)
         tatgaCTGCGATGCCATCGAAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAgg                                                                                                                                             >  6:116646/5‑148 (MQ=255)
            cACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAAt                                                                                                                                         <  5:311982/149‑1 (MQ=255)
               ttCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGATCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGc                                                                                                                                      >  6:374545/1‑149 (MQ=255)
                 cGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGc                                                                                                                                      <  8:40376/147‑1 (MQ=255)
                   aTGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCAt                                                                                                                                  <  6:375759/149‑1 (MQ=255)
                       cATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGccgct                                                                                                                                  <  1:144459/145‑4 (MQ=255)
                       cATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGccgct                                                                                                                                  >  2:144459/1‑142 (MQ=255)
                       cATCGCAGTCATACACCCGCAGGTTTATCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGcc                                                                                                                               >  1:9753/1‑148 (MQ=255)
                           gCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGCCGCTGCCgcac                                                                                                                           <  6:125743/148‑4 (MQ=38)
                           gCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAg                                                                                                                           <  3:188245/148‑1 (MQ=255)
                           gCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGc                                                                                                                          <  3:251185/149‑1 (MQ=255)
                             aGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGcc                                                                                                                         >  5:333103/1‑148 (MQ=255)
                                   acacCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGAtctc                                                                                                                   >  4:235218/1‑148 (MQ=255)
                                    cacCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTcc                                                                                                                  <  6:323021/148‑1 (MQ=255)
                                     acCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCaa                                                                                                                <  2:293579/149‑1 (MQ=255)
                                      cccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAg                                                                                                               <  5:48994/149‑1 (MQ=255)
                                      cccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAg                                                                                                               <  2:363742/149‑1 (MQ=255)
                                       ccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAgg                                                                                                              >  6:281052/1‑149 (MQ=255)
                                        cGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGt                                                                                                             >  6:378491/1‑149 (MQ=255)
                                         gCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTg                                                                                                            <  2:9753/149‑1 (MQ=255)
                                           aGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCATCACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTgcg                                                                                                          >  1:279763/1‑149 (MQ=255)
                                              tttttCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgg                                                                                                       >  3:263240/1‑149 (MQ=255)
                                                 ttCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgggg                                                                                                     <  3:372929/148‑1 (MQ=255)
                                                  tCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggc                                                                                                    <  7:226915/148‑1 (MQ=255)
                                                   cAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggcc                                                                                                   <  5:378491/148‑1 (MQ=255)
                                                      ccccTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaac                                                                                               >  5:128353/1‑149 (MQ=255)
                                                        ccTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCTGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGGGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaacgg                                                                                             <  6:1640/149‑1 (MQ=255)
                                                        ccTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaacgg                                                                                             >  8:345417/1‑149 (MQ=255)
                                                         cTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGTGGggccgaacggc                                                                                            <  3:235218/149‑1 (MQ=255)
                                                         cTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaacggc                                                                                            <  5:281052/149‑1 (MQ=255)
                                                               cAGCACGCCATGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAacac                                                                                     >  4:125707/1‑149 (MQ=255)
                                                                aGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAacac                                                                                     <  4:263240/149‑1 (MQ=255)
                                                                    cGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTgg                                                                                 >  6:295249/1‑149 (MQ=255)
                                                                    cGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTgg                                                                                 <  2:225471/149‑1 (MQ=255)
                                                                    cGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTgg                                                                                 >  2:380791/1‑149 (MQ=255)
                                                                     gCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTACACACAGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCAAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCAGAACGGCCGAACACCTGGt                                                                                >  1:372483/1‑149 (MQ=255)
                                                                            agtagcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTc                                                                         >  8:131044/4‑149 (MQ=255)
                                                                            agcagcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCTGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGACCGGCCGAACCCCTGGTCGTAGTc                                                                         >  8:212187/1‑149 (MQ=255)
                                                                            agcagcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTc                                                                         >  6:220954/1‑149 (MQ=255)
                                                                             gcagcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTTCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTc                                                                         >  6:157542/1‑148 (MQ=255)
                                                                             gcagcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTc                                                                         >  6:112645/1‑148 (MQ=255)
                                                                              cagcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCa                                                                       >  4:369553/1‑149 (MQ=255)
                                                                               agcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCAc                                                                      >  5:304499/1‑149 (MQ=255)
                                                                                 cgcgcCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCg                                                                    >  7:43375/1‑149 (MQ=255)
                                                                                  gcgcCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCgg                                                                   >  3:213183/1‑149 (MQ=255)
                                                                                   cgcCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCggg                                                                  <  2:341571/149‑1 (MQ=255)
                                                                                    gcCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGc                                                                 >  8:133539/1‑149 (MQ=255)
                                                                                        acacCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATc                                                             <  7:317394/149‑1 (MQ=255)
                                                                                        acacCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATc                                                             <  3:125707/149‑1 (MQ=255)
                                                                                         cacCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCa                                                            <  2:279763/149‑1 (MQ=255)
                                                                                          acCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAg                                                           <  8:43375/149‑1 (MQ=255)
                                                                                          acCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAg                                                           >  2:380163/1‑149 (MQ=255)
                                                                                           cccGGGCGGTGCACTTCAACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGTTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTAtgatc                                                          <  2:372483/149‑5 (MQ=255)
                                                                                            ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  8:272580/1‑148 (MQ=255)
                                                                                            ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  3:282947/1‑148 (MQ=255)
                                                                                            ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  7:97577/1‑148 (MQ=255)
                                                                                            ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  5:12363/1‑148 (MQ=255)
                                                                                            ccGGGCGGTGCACTTCCACACCGAAGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGCTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  6:211599/1‑148 (MQ=255)
                                                                                             cGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  5:278403/1‑147 (MQ=255)
                                                                                                  ggTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGc                                                   <  5:116646/149‑1 (MQ=255)
                                                                                                    tGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAg                                                 >  5:115838/1‑149 (MQ=255)
                                                                                                       aCTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGAt                                              <  7:345417/149‑1 (MQ=255)
                                                                                                                 cGATGAAATCCTGCTCAGGTGCGGCAGCCGCCATTTCCAGCAGGGAAAGGACCAAGCCGAAGCCGATCTCCAAGGTGCGAGGGGCCGAACGGCCGAACACCTGCACGAAGTCAACCCGGCTAACAGCCAGTGGCAGGATGTACAGCgtg                                    >  7:40441/1‑147 (MQ=255)
                                                                                                                   aTGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGgcc                                  >  7:16812/1‑149 (MQ=255)
                                                                                                                    tGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGgccg                                 >  1:271825/1‑149 (MQ=255)
                                                                                                                       aaTCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGgccgcc                               <  8:97577/148‑1 (MQ=255)
                                                                                                                        aTCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCtt                             >  5:312663/1‑149 (MQ=255)
                                                                                                                         tCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTg                            <  6:12363/149‑1 (MQ=255)
                                                                                                                           cTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTTCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTgg                           <  5:157542/148‑1 (MQ=255)
                                                                                                                               tCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGGGGCAGGATGTACCGCGGGCCGCCGTgggcgag                      >  3:295316/1‑144 (MQ=255)
                                                                                                                                 aGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCAAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGAc                    >  6:274399/1‑149 (MQ=255)
                                                                                                                                    ttcGGCATGGGCCATTCCCTGCTGGAAATGGCCGCTGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACcgc                 >  8:161277/3‑149 (MQ=38)
                                                                                                                                    tGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACcgc                 >  2:379597/1‑149 (MQ=255)
                                                                                                                                       ggCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGtt              <  8:16812/149‑1 (MQ=255)
                                                                                                                                         catgggCCATTCCCTGCTGGAAATGGCCCATGCCGAAGCCGATCTCCAAGGTGAGCGGGGCCGAACGCCCTAACACCTGGTCGTAGTCCACCGGGCTATCAGACAGTGGCAGGATGTACAGCGGGCCCCCTTGATCCAGACCGCGATg             >  8:78477/5‑148 (MQ=255)
                                                                                                                                           tgggCCATTCCCTGCTGGAAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGCACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTg          >  1:156261/3‑149 (MQ=255)
                                                                                                                                           gCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTg          <  4:282947/149‑1 (MQ=255)
                                                                                                                                           gCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTg          <  2:271825/149‑1 (MQ=255)
                                                                                                                                              gCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGcc       <  6:304499/149‑1 (MQ=255)
                                                                                                                                               ccATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCt      >  7:296632/1‑149 (MQ=255)
                                                                                                                                               ccATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTAGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCt      >  6:62240/1‑149 (MQ=255)
                                                                                                                                                cATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCtt     >  3:213269/1‑149 (MQ=255)
                                                                                                                                                 aTTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTc    >  5:244414/1‑149 (MQ=255)
                                                                                                                                                  tttCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCg   <  4:213183/149‑1 (MQ=255)
                                                                                                                                                  tttCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCg   <  6:312663/149‑1 (MQ=255)
                                                                                                                                                    tCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg  <  5:295249/148‑1 (MQ=255)
                                                                                                                                                    tCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg  <  1:380791/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGG  >  NC_002947/5825083‑5825378

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: