Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,825,234 A→T 7.1% 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 7.1% 162.0 / 3.6 56L80Q (CTG→CAG) trmBtRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base A (30/22);  new base T (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTC  >  NC_002947/5825087‑5825380
                                                                                                                                                   |                                                                                                                                                  
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:25733/1‑149 (MQ=255)
              gATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCa                                                                                                                                     <  2:274339/149‑1 (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)
                                     gCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTg                                                                                                              >  6:98641/1‑149 (MQ=255)
                                       aGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAAt                                                                                                                                           <  5:75542/118‑1 (MQ=255)
                                       aGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAAt                                                                                                                                           >  6:75542/1‑118 (MQ=255)
                                          tttttCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgg                                                                                                         <  2:68220/149‑1 (MQ=255)
                                             ttCAGGCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGGCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgggg                                                                                                       <  2:40357/148‑1 (MQ=255)
                                             ttCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCgggg                                                                                                       <  8:181080/148‑1 (MQ=255)
                                                 gCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaa                                                                                                  >  7:308113/1‑149 (MQ=255)
                                                 gCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaa                                                                                                  >  1:219523/1‑149 (MQ=255)
                                                       gggTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGggccgaacggccg                                                                                            <  7:332517/149‑1 (MQ=255)
                                                          tCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGCCGCTGCCgcac                                                                                                                             >  6:369574/1‑110 (MQ=25)
                                                          tCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCATGGGCCATTCCCTGCTGGAAATGGCCGCTGCCgcac                                                                                                                             <  5:369574/113‑4 (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)
                                                                                                                                                  cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg    >  2:368741/1‑146 (MQ=255)
                                                                                                                                                  cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCgg    <  1:368741/146‑1 (MQ=255)
                                                                                                                                                  cagcagGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTc  <  6:144091/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGTTCAGCACTTCGATGGCATCGCAGTCATACACCCGCAGGTTTTTCAGCCCCTGGGTCAGCACGCCATTGAGCAGCGCGCCAACACCCGGGCGGTGCACTTCCACACCGATGAAATCCTGCTCAGGTGCGGCAGCGGCCATTTCCAGCAGGGAATGGCCCATGCCGAAGCCGATCTCCAAGGTGCGCGGGGCCGAACGGCCGAACACCTGGTCGTAGTCCACCGGGCTATCAGCCAGTGGCAGGATGTACAGCGGGCCGCCTTGGTCCAGACCGCGTTGCTGGCCTTCGGTC  >  NC_002947/5825087‑5825380

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