Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,825,237 A→T 7.9% L79Q (CTG→CAG)  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,2370AT 7.9% 257.8 / 8.0 87L79Q (CTG→CAG) trmBtRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base A (43/37);  new base T (5/2);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCATG  >  NC_002947/5825090‑5825383
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                    <  2:363742/149‑1 (MQ=255)
                               cccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAg                                                                                                                    <  5:48994/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:380791/1‑149 (MQ=255)
                                                             cGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTgg                                                                                      <  2:225471/149‑1 (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                                                                >  2:380163/1‑149 (MQ=255)
                                                                                   acCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAg                                                                <  8:43375/149‑1 (MQ=255)
                                                                                    cccGGGCGGTGCACTTCAACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGTTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTAtgatc                                                               <  2:372483/149‑5 (MQ=255)
                                                                                     ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                               >  5:12363/1‑148 (MQ=255)
                                                                                     ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                               >  8:272580/1‑148 (MQ=255)
                                                                                     ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                               >  3:282947/1‑148 (MQ=255)
                                                                                     ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                               >  7:97577/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        <  6:312663/149‑1 (MQ=255)
                                                                                                                                           tttCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCg        <  4:213183/149‑1 (MQ=255)
                                                                                                                                             tCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg       <  5:295249/148‑1 (MQ=255)
                                                                                                                                             tCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg       <  1:380791/148‑1 (MQ=255)
                                                                                                                                               cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGTCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGTCCTTCGGTc     <  3:369553/148‑1 (MQ=255)
                                                                                                                                               cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGtt                   >  4:357206/1‑134 (MQ=255)
                                                                                                                                               cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGtt                   <  3:357206/134‑1 (MQ=255)
                                                                                                                                               cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCa    >  7:17975/1‑149 (MQ=255)
                                                                                                                                                 gcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCATg  <  6:115838/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCATG  >  NC_002947/5825090‑5825383

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: