Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,825,234 A→T 6.9% L80Q (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,2340AT 6.9% 261.3 / 6.7 87L80Q (CTG→CAG) trmBtRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base A (44/37);  new base T (4/2);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

TTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCA  >  NC_002947/5825090‑5825381
                                                                                                                                                |                                                                                                                                                   
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)
                                                        cAGCACGCCA‑TGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAacac                                                                                        >  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                                                             >  8:272580/1‑148 (MQ=255)
                                                                                     ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                             >  5:12363/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                 <  3:357206/134‑1 (MQ=255)
                                                                                                                                               cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGtt                 >  4:357206/1‑134 (MQ=255)
                                                                                                                                               cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCa  >  7:17975/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTCA  >  NC_002947/5825090‑5825381

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