Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,980 A→G 5.2% intergenic (‑599/+32) trmF ← / ← trmE tRNA uridine 5‑carboxymethylaminomethyl modification protein/GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,9800AG 5.2% 226.8 / 2.2 77intergenic (‑599/+32)trmF/trmEtRNA uridine 5‑carboxymethylaminomethyl modification protein/GTPase
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base G (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

CCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGAAGC  >  NC_002947/4837‑5126
                                                                                                                                               |                                                                                                                                                   
cccGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgc                                                                                                                                                >  7:48038/1‑149 (MQ=255)
 ccGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgcg                                                                                                                                               <  8:83550/149‑1 (MQ=255)
 ccGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGGCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgcg                                                                                                                                               <  3:327208/149‑1 (MQ=255)
  cgcgAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGACCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgcgg                                                                                                                                              >  6:93815/1‑149 (MQ=255)
   gcgAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACTAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAgcggc                                                                                                                                             <  4:385346/149‑1 (MQ=255)
     gAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggctt                                                                                                                                           <  2:268389/149‑3 (MQ=255)
        aGGATGTGTAAATCCACAGAGATCTGGTGCTTGTAAGAGATTCATCACGGGTCCACAAGCTCCCACAGGTGCTTCTGTCAGAGGCCTGCATTCATCTTCTCAAACGAATTTCAGGCAAAAAAAAGCAGCTCCTTAGGTAgaggcttttt                                                                                                                                        >  4:250819/1‑144 (MQ=255)
         ggAGGGGTACATCCACAGAGATCTGGTGCTTGTCCGCGCTGCATCACGGGTCCACCCGCTCCCACAGGTGCTGGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCTTAGGGAgcggctttttt                                                                                                                                       <  5:193261/149‑7 (MQ=255)
          gaggtgTACATCCACAGAGATCTGTGGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATGCATCTTCTGCAACTAATTTCAGCCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGc                                                                                                                                      <  6:198584/144‑1 (MQ=255)
               gTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGc                                                                                                                                <  8:48038/149‑1 (MQ=255)
                  cATCCACATCGATCTGGTGCTTGTAAGCGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTg                                                                                                                             <  4:289758/149‑1 (MQ=255)
                  cATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGAAAAAAAAAGCCGCTCCGTAAGGAGCGGCGCTTGGCCTTGCGt                                                                                                                              <  6:302220/148‑1 (MQ=255)
                       acaGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgt                                                                                                                        <  3:173795/149‑1 (MQ=255)
                        caGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCtgtg                                                                                                                       <  1:227981/149‑1 (MQ=255)
                             aTCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGAtt                                                                                                                   <  2:201949/148‑1 (MQ=255)
                                tGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACtt                                                                                                               >  8:240258/1‑149 (MQ=255)
                                  gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                             >  4:199384/1‑149 (MQ=255)
                                  gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                             <  4:353641/149‑1 (MQ=255)
                                  gTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGc                                                                                                             <  4:347942/149‑1 (MQ=255)
                                     cTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGa                                                                                                          >  5:346145/1‑149 (MQ=255)
                                     cTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGa                                                                                                          >  7:165458/1‑149 (MQ=255)
                                     cTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGa                                                                                                          >  5:209526/1‑149 (MQ=255)
                                        gTAACATCTTCATCACGGGGCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATg                                                                                                        <  4:73018/148‑1 (MQ=255)
                                         tAAGAGGTTCATCACGGGTCCAACCGCTCCCACAGGTGCTTGTGTCAGCGGCCTGCATTCATCTTCTGGAACGAATTTCACGCAAAAAAAAGCCGCTCCCTAAGGATCGGCGCTTTGCCTTCTCGTGGCTGTGGATTACTTGCCGATg                                                                                                        <  1:4951/148‑1 (MQ=255)
                                           agagCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCGGTGGGTTACTTGCCGAGGCAGa                                                                                                    >  8:257669/1‑149 (MQ=255)
                                                  catcaCGGGACCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTg                                                                                               <  5:93815/147‑1 (MQ=255)
                                                  cataaCGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTgg                                                                                              >  7:151819/1‑148 (MQ=255)
                                                     caCGGGTCCACCCGCTCCCACAGGTGCATGGCTCAGCGGCCTGCATTCAACTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCGAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCGGgaga                                                                                           >  2:282299/1‑148 (MQ=255)
                                                          gTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATcc                                                                                     >  1:73555/1‑149 (MQ=255)
                                                          gTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATcc                                                                                     >  7:51357/1‑149 (MQ=255)
                                                                ccGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCaa                                                                               <  8:51357/149‑1 (MQ=255)
                                                                  gCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAg                                                                              >  8:218981/1‑148 (MQ=255)
                                                                      ccACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAg                                                                           <  7:35568/147‑1 (MQ=255)
                                                                        acaGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCa                                                                       >  7:39148/1‑149 (MQ=255)
                                                                          aGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATc                                                                     <  6:209526/149‑1 (MQ=255)
                                                                            gTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATcc                                                                    >  6:397644/1‑148 (MQ=255)
                                                                                  gggTCCGGGGCCTGCATTCATCTTCTGCAACTAATTTCCGGAAAAAAAACGCGGCTCCCTATGGCGGGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGa                                                             <  2:395111/149‑1 (MQ=255)
                                                                                   ggTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGa                                                             <  3:199384/148‑1 (MQ=255)
                                                                                    gTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAAt                                                           >  8:252316/1‑149 (MQ=255)
                                                                                        gCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTc                                                       <  7:240258/149‑1 (MQ=255)
                                                                                        gCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTc                                                       <  6:371375/149‑1 (MQ=255)
                                                                                           gCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGACCAAGCAGGTCATCCGAGATGAATGCTCCtt                                                    >  6:240793/1‑148 (MQ=255)
                                                                                              tGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGa                                                 >  4:56819/1‑149 (MQ=255)
                                                                                               gCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGa                                                 >  3:261711/1‑148 (MQ=255)
                                                                                               gCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGa                                                 >  4:213447/1‑148 (MQ=255)
                                                                                               gCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAt                                                >  5:149034/1‑149 (MQ=255)
                                                                                                cATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAtc                                               <  5:397644/149‑1 (MQ=255)
                                                                                                   tCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTcc                                            >  3:351021/1‑149 (MQ=255)
                                                                                                   tCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTcc                                            >  3:149035/1‑149 (MQ=255)
                                                                                                          ctGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAgc                                      <  8:165458/148‑1 (MQ=255)
                                                                                                           tGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCt                                                     >  5:268653/1‑132 (MQ=255)
                                                                                                           tGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAgcg                                     >  7:12656/1‑148 (MQ=255)
                                                                                                             cAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCt                                                     >  6:57204/1‑130 (MQ=255)
                                                                                                             cAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCt                                                     <  5:57204/130‑1 (MQ=255)
                                                                                                             cAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCAt                                  >  8:206559/1‑149 (MQ=255)
                                                                                                             cAACGAATTTCAGGCAAAAAAAAGCCGCTCACTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCTATGCAGAAGCTGGAGAAAATCCGACCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCAt                                  >  8:220663/1‑149 (MQ=255)
                                                                                                              aaCGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCt                                                     <  5:291110/129‑1 (MQ=255)
                                                                                                              aaCGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCt                                                     >  6:291110/1‑129 (MQ=255)
                                                                                                              aaCGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATg                                 <  2:73555/149‑1 (MQ=255)
                                                                                                                cGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCt                               >  2:187407/1‑149 (MQ=255)
                                                                                                                        aGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTg                        <  3:56819/148‑1 (MQ=255)
                                                                                                                         ggCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGc                      <  3:213447/149‑1 (MQ=255)
                                                                                                                         ggCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGc                      <  6:149034/149‑1 (MQ=255)
                                                                                                                          gCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAt                                                >  4:391023/1‑122 (MQ=255)
                                                                                                                          gCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGAt                                                <  3:391023/122‑1 (MQ=255)
                                                                                                                          gCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCg                     >  4:257594/1‑149 (MQ=255)
                                                                                                                            aaaaaaaaGCCGCTCCCTAAGGAGCGGCTATTGTACTTTAATGGCTATGTATTACTTGACGATTAAGAACCTGGAGAAAATAAGAACAAGCAGGTCATACAAGCTGAAAGATAATTTGAACTCCCCAAGACCAAGCTTTGcattgatt                    >  6:106469/1‑141 (MQ=255)
                                                                                                                             aaaaaaaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTaa                  <  4:149035/149‑1 (MQ=255)
                                                                                                                               caagcgCCGCTCCTTAGGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGTTGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGt                >  7:84527/6‑149 (MQ=255)
                                                                                                                                aaaaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTc               <  7:252316/149‑1 (MQ=255)
                                                                                                                                 aaaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGAGCATGCTGTGCCTGGCGTAAGTcc              >  3:231131/1‑149 (MQ=255)
                                                                                                                                  aaGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAATCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCt             <  8:151819/149‑1 (MQ=255)
                                                                                                                                    gCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCa           <  8:199875/149‑1 (MQ=255)
                                                                                                                                          accTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAg      <  4:86617/147‑1 (MQ=255)
                                                                                                                                           ccTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAg      >  5:106323/1‑147 (MQ=255)
                                                                                                                                           ccTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGaa    <  3:257594/149‑1 (MQ=255)
                                                                                                                                             tAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGAAGc  >  1:390944/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCCGCGAAAGGATGGGTACATCCACAGAGATCTGGTGCTTGTAAGAGCTTCATCACGGGTCCACCCGCTCCCACAGGTGCTTGGGTCAGCGGCCTGCATTCATCTTCTGCAACGAATTTCAGGCAAAAAAAAGCCGCTCCCTAAGGAGCGGCGCTTGGCCTTGCGTGGCTGTGGATTACTTGCCGATGCAGAAGCTGGAGAAAATCCGCCCAAGCAGGTCATCCGAGCTGAATGCTCCTGTGATCTCCCCCAGCGCATGCTGTGCCTGGCGTAAGTCCTCAGCCAGAAGC  >  NC_002947/4837‑5126

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: