Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,825,231 T→C 6.6% 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 6.6% 174.1 / 2.0 60E81G (GAA→GGA) trmBtRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGG  >  NC_002947/5825084‑5825378
                                                                                                                                                   |                                                                                                                                                   
cAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTc                                                                                                                                                    <  7:342982/149‑1 (MQ=255)
   cGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCcag                                                                                                                                                 <  2:68968/149‑1 (MQ=255)
     gTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCcagca                                                                                                                                               <  1:158806/149‑1 (MQ=255)
     gTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCcagca                                                                                                                                               <  4:12560/149‑1 (MQ=255)
        cAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGTGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAgg                                                                                                                                             >  4:241132/1‑148 (MQ=255)
         aGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCATGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTgga                                                                                                                                            >  4:193126/1‑147 (MQ=255)
           cACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAAt                                                                                                                                         >  2:133603/1‑149 (MQ=255)
            aCTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATg                                                                                                                                        >  1:40357/1‑149 (MQ=255)
                 gATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCa                                                                                                                                   <  2:274339/149‑1 (MQ=255)
                 gATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCa                                                                                                                                   >  2:25733/1‑149 (MQ=255)
                 gATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCa                                                                                                                                   >  7:9598/1‑149 (MQ=255)
                  aTGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCAt                                                                                                                                  <  1:25733/149‑1 (MQ=255)
                    ggCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGc                                                                                                                                >  2:31283/1‑149 (MQ=255)
                     gCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCATGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGccgctgc                                                                                                                                <  3:193126/148‑6 (MQ=38)
                        tCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGaa                                                                                                                            <  4:257115/149‑1 (MQ=255)
                              tCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGTGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGAt                                                                                                                      <  3:241132/149‑1 (MQ=255)
                              tCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCGCATGCCGAAGCCGAt                                                                                                                      >  3:255222/1‑149 (MQ=255)
                              tCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGAt                                                                                                                      >  1:339229/1‑149 (MQ=255)
                               cATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGAt                                                                                                                      <  2:171764/148‑1 (MQ=255)
                                 tACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGAtctc                                                                                                                   >  1:68220/1‑149 (MQ=255)
                                   cacCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTcc                                                                                                                  <  7:183557/148‑1 (MQ=255)
                                     cccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAg                                                                                                               <  1:31283/149‑1 (MQ=255)
                                     cccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAg                                                                                                               <  1:53354/149‑1 (MQ=255)
                                      ccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTc                                                                                                                                                    <  5:47986/111‑1 (MQ=255)
                                      ccGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTc                                                                                                                                                    >  6:47986/1‑111 (MQ=255)
                                        gCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTg                                                                                                            >  6:98641/1‑149 (MQ=255)
                                          aGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAAt                                                                                                                                         >  6:75542/1‑118 (MQ=255)
                                          aGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAAt                                                                                                                                         <  5:75542/118‑1 (MQ=255)
                                             tttttCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgg                                                                                                       <  2:68220/149‑1 (MQ=255)
                                                ttCAGGCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGGCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgggg                                                                                                     <  2:40357/148‑1 (MQ=255)
                                                ttCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgggg                                                                                                     <  8:181080/148‑1 (MQ=255)
                                                    gCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaa                                                                                                >  1:219523/1‑149 (MQ=255)
                                                    gCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaa                                                                                                >  7:308113/1‑149 (MQ=255)
                                                          gggTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaacggccg                                                                                          <  7:332517/149‑1 (MQ=255)
                                                             tCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGCCGCTGCCgcac                                                                                                                           <  5:369574/113‑4 (MQ=25)
                                                             tCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGCCGCTGCCgcac                                                                                                                           >  6:369574/1‑110 (MQ=25)
                                                                   cGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTgg                                                                                 <  8:145008/149‑1 (MQ=255)
                                                                       aTTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCg                                                                              >  3:283352/1‑148 (MQ=255)
                                                                              agcGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCAc                                                                      <  1:133603/149‑1 (MQ=255)
                                                                                      aaCACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTAt                                                              >  7:213439/1‑149 (MQ=255)
                                                                                       acacCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTAt                                                              >  3:61661/1‑148 (MQ=255)
                                                                                       acacCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATc                                                             <  2:339229/149‑1 (MQ=255)
                                                                                          cccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATAAGc                                                          >  8:232996/1‑149 (MQ=255)
                                                                                           ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  4:352788/1‑148 (MQ=255)
                                                                                           ccGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGc                                                          >  6:217232/1‑148 (MQ=255)
                                                                                                  gTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCa                                                  <  5:217232/149‑1 (MQ=255)
                                                                                                    gCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAg                                                 >  5:193963/1‑148 (MQ=255)
                                                                                                        ttCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATg                                             >  2:173192/1‑148 (MQ=255)
                                                                                                        ttCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGACGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCatgatg                                             >  5:219735/1‑148 (MQ=255)
                                                                                                            acacCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAg                                        >  5:252870/1‑149 (MQ=255)
                                                                                                               ccGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCAAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCgg                                     <  8:9598/149‑1 (MQ=255)
                                                                                                                  aTGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGgcc                                  >  7:326073/1‑149 (MQ=255)
                                                                                                                        tCCTGCTCAGGTGCGGGAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTg                            <  5:98641/149‑1 (MQ=255)
                                                                                                                           tGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTc                         >  1:238628/1‑149 (MQ=255)
                                                                                                                            gCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTc                         >  5:57831/1‑148 (MQ=255)
                                                                                                                            gCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTc                         >  2:103951/1‑148 (MQ=255)
                                                                                                                                 ggTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGAcc                   <  2:219523/149‑1 (MQ=255)
                                                                                                                                  gTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACcg                  <  4:61661/149‑1 (MQ=255)
                                                                                                                                         aGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCt           <  6:193963/149‑1 (MQ=255)
                                                                                                                                            ggCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGc        >  5:165572/1‑149 (MQ=255)
                                                                                                                                                 tttCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCg   >  5:219632/1‑149 (MQ=255)
                                                                                                                                                  ttCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg  >  5:144091/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGG  >  NC_002947/5825084‑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: ‑