Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,824,714 T→A 8.6% intergenic (+138/+36) PP_5102 → / ← trmB hypothetical protein/tRNA (guanosine(46)‑N7)‑methyltransferase TrmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,824,7140TA 8.6% 196.5 / 5.8 67intergenic (+138/+36)PP_5102/trmBhypothetical protein/tRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base A (3/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGT  >  NC_002947/5824569‑5824854
                                                                                                                                                 |                                                                                                                                            
gcgGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAg                                                                                                                                           <  4:62723/149‑1 (MQ=255)
  ggCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAggg                                                                                                                                         <  6:301221/149‑1 (MQ=255)
          gAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTa                                                                                                                                 >  7:109891/1‑149 (MQ=255)
          gAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTa                                                                                                                                 >  5:257561/1‑149 (MQ=255)
             tGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAcc                                                                                                                               >  4:208990/1‑148 (MQ=255)
                ccccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAACGaa                                                                                                                            >  8:303313/1‑148 (MQ=255)
                 cccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGaa                                                                                                                            >  6:174922/1‑147 (MQ=255)
                  ccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAg                                                                                                                         >  7:129744/1‑149 (MQ=255)
                  ccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAg                                                                                                                         >  4:236766/1‑149 (MQ=255)
                   cTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAgg                                                                                                                        <  6:194513/149‑1 (MQ=255)
                      gTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGcc                                                                                                                      <  1:48255/148‑1 (MQ=255)
                      gTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGc                                                                                                                       <  1:64437/147‑1 (MQ=255)
                      gTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCg                                                                                                                     >  3:288216/1‑149 (MQ=255)
                      gTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGcc                                                                                                                      <  2:202476/147‑1 (MQ=255)
                       tAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCg                                                                                                                     <  6:257561/148‑1 (MQ=255)
                         ggAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTa                                                                                                                   <  7:303313/148‑1 (MQ=255)
                          gAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTAt                                                                                                                  >  6:112333/1‑148 (MQ=255)
                            gctgATTCAACCGCGAAAGGGCAGGCCCTGCTTGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGTTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGtttttt                                                                                                                                <  5:9648/129‑3 (MQ=18)
                            gCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAg                                                                                                               >  3:253022/1‑149 (MQ=255)
                             cGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGAGGTTTTACCGAAAAGGCCGTATCAGc                                                                                                              >  2:104853/1‑149 (MQ=255)
                                    aaCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCaa                                                                                                       >  8:196088/1‑149 (MQ=255)
                                      ccGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATc                                                                                                     <  8:109891/149‑1 (MQ=255)
                                          gAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGtttt                                                                                                                                  >  1:6945/1‑116 (MQ=18)
                                          gAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGtttt                                                                                                                                  <  2:6945/116‑1 (MQ=18)
                                             tgggCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGGTACCGTTTTAGCGTTTAAACCCTCTCCTACAAAGAAAAAACCTCCCTTTTCAGGGCGGTTTTACCGAAAAGGCGGTATCAGCCGGTCAATCCACCTtc                                                                                              <  3:77892/148‑1 (MQ=37)
                                              gggCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTtct                                                                                             <  8:183632/149‑1 (MQ=255)
                                               ggccggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTtctc                                                                                            <  8:129744/149‑1 (MQ=255)
                                                  cggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCg                                                                                                                     >  7:303290/1‑121 (MQ=255)
                                                  cggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCg                                                                                                                     <  8:303290/121‑1 (MQ=255)
                                                   ggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAGGGCCGTATCAGCCGCTCAATCCACCTTCTCGAAc                                                                                        <  8:260489/149‑1 (MQ=255)
                                                               cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTccccc                                                                            >  1:188836/1‑147 (MQ=255)
                                                                        cccGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGc                                                                   >  8:283544/1‑149 (MQ=255)
                                                                         ccGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTccc                                                                              >  1:135042/1‑137 (MQ=255)
                                                                                           gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCgcggc                                                <  5:174922/149‑1 (MQ=255)
                                                                                             cgcTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGt                                              <  4:74580/149‑1 (MQ=255)
                                                                                                  ccGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGaa                                         >  2:67502/1‑149 (MQ=255)
                                                                                                  ccGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGaa                                         <  7:283544/149‑1 (MQ=255)
                                                                                                   cGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTTACCGCGGCGTTCGAAc                                        >  1:250244/1‑149 (MQ=255)
                                                                                                    gTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACt                                       >  4:151190/1‑149 (MQ=255)
                                                                                                       ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATGCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCCGCGGCGCACCGCGGCGTTCGAACTggg                                    >  4:221451/1‑149 (MQ=255)
                                                                                                       ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTTACCGCGGCGTTCGAACTTg                                     <  2:250244/148‑1 (MQ=255)
                                                                                                       ttCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                     <  3:151190/148‑1 (MQ=255)
                                                                                                         cgcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGt                                   >  7:86124/1‑148 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAg                                                              <  1:10828/120‑1 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAg                                                              >  2:10828/1‑120 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTg                                  <  7:196088/148‑1 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGa                                 >  5:306272/1‑149 (MQ=255)
                                                                                                              gTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTACAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAAACCACCTTCACCAACTTCGAGTCCCACACCCCCTGCCCCAGCCGCTAACCGCGGCCTTCGAACTTCGTGATCgg                             >  7:60013/1‑149 (MQ=255)
                                                                                                                 aaCCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGgcg                          <  3:236766/149‑1 (MQ=255)
                                                                                                                       cTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTc                    >  6:51088/1‑149 (MQ=255)
                                                                                                                          cTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCggg                 <  1:104853/149‑1 (MQ=255)
                                                                                                                          cTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCggg                 >  2:233612/1‑149 (MQ=255)
                                                                                                                          cTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCg                   <  4:253022/147‑1 (MQ=255)
                                                                                                                            aCAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgc               >  1:276567/1‑149 (MQ=255)
                                                                                                                            aCAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgc               >  3:162958/1‑149 (MQ=255)
                                                                                                                             cAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgcg              >  5:107252/1‑149 (MQ=255)
                                                                                                                               aaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGc                                                                                                              >  2:102975/1‑51 (MQ=255)
                                                                                                                                aGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGt                                                                     >  5:159588/1‑91 (MQ=255)
                                                                                                                                aGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGt                                                                     <  6:159588/91‑1 (MQ=255)
                                                                                                                                 ggtaaaaCCGCCCTGAAAAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAg                                                                                  >  4:99762/4‑77 (MQ=25)
                                                                                                                                 ggtaaaaCCGCCCTGAAAAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAg                                                                                  <  3:99762/74‑1 (MQ=25)
                                                                                                                                     aaaCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACAt      <  4:162958/149‑1 (MQ=255)
                                                                                                                                        ccGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCg                                                           <  5:80039/93‑1 (MQ=255)
                                                                                                                                        ccGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCg                                                           >  6:80039/1‑93 (MQ=255)
                                                                                                                                         cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGt  <  6:306272/149‑1 (MQ=255)
                                                                                                                                         cGCCCTTTCAGGGCGGTTTTTCCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGcac        <  4:300354/142‑1 (MQ=255)
                                                                                                                                         cGCCCTTTCAGGGCGGTTTTTCCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGcac        >  3:300354/1‑142 (MQ=255)
                                                                                                                                              ttttCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCg                                                           >  1:151812/1‑87 (MQ=255)
                                                                                                                                              ttttCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCg                                                           <  2:151812/87‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GCGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGT  >  NC_002947/5824569‑5824854

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: