Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,754,612 T→G 26.4% K85N (AAA→AAC U712_19355 ← Putative ribosomal RNA large subunit methyltransferase YwbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,754,6120TG26.4% 83.9 / 8.0 54K85N (AAA→AACU712_19355Putative ribosomal RNA large subunit methyltransferase YwbD
Reads supporting (aligned to +/- strand):  ref base T (12/27);  new base G (12/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 5.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGTACCATTGAATCAGCAGGTAGCCGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCCGCGCGAGAAACTCGCCGCCGCTCTCAGACACGATATC  >  CP006881/3754476‑3754749
                                                                                                                                        |                                                                                                                                           
tGCTGTACCATTGAATCAGCAGGTAGCCGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCTTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCttt                                                                                                                                           >  5:370582/1‑139 (MQ=255)
 gCTGTACCATTGAATCAGCAGGTAGCCGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCttt                                                                                                                                           <  5:336246/138‑1 (MQ=255)
       ccATTGAATCAGCAGGTAGCCGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTTTTTAt                                                                                                                                   >  2:313264/1‑139 (MQ=255)
                         gCCGTCATATTAGTCGATGGTGACGCCGCCGACACCGTCACCCTGACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCGTGGGCCGCTTTTGTTAATTTGGTGAGAAAAAATGCt                                                                                                                <  1:33225/139‑1 (MQ=255)
                         gCCGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCt                                                                                                                <  3:116621/139‑1 (MQ=255)
                         agcGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTAGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCt                                                                                                                <  3:241827/137‑1 (MQ=255)
                             tcctaataGTCGCTGGTGACGCCGCCGACACCGGCCCCCTCACCGTTGAACCGGCGGACAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCGGGGGCCGCGTTTGTTAATTTGCTGAGAAAACATGCTTGGt                                                                                                            <  8:238279/136‑1 (MQ=255)
                             tCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTCTAACTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGGTTATTTTCTGAGAAAAAATTCTTggg                                                                                                            >  7:57211/1‑138 (MQ=255)
                               ataataGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCa                                                                                                          <  1:51329/139‑1 (MQ=255)
                                  ataGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAAt                                                                                                        <  4:210954/138‑1 (MQ=255)
                                  ataGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATc                                                                                                       <  8:335706/139‑1 (MQ=255)
                                      tCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGtt                                                                                                   <  7:276577/139‑1 (MQ=255)
                                      tCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGtt                                                                                                   <  3:192709/139‑1 (MQ=255)
                                      tCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTTATTTTGTGGTAAAAAATGCTTGGGGAATATGtt                                                                                                   >  4:345984/1‑139 (MQ=255)
                                      tCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTGGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGtt                                                                                                   <  3:283028/139‑1 (MQ=255)
                                      tCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTAGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTTTTTATTTTCTGAGAAAAAATGCTTTGTCAATCTgct                                                                                                   >  4:241827/1‑137 (MQ=255)
                                        gATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTTTTTAATTGGTGACAAAAAATTATTGGTGAATCTGttga                                                                                                 >  6:124155/1‑137 (MQ=255)
                                             tGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTGATGttt                                                                                            >  5:104631/1‑139 (MQ=255)
                                              gACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTa                                                                                           <  4:314669/139‑1 (MQ=255)
                                                cgccgcCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTAtt                                                                                         <  6:299672/139‑1 (MQ=255)
                                                   cgccgAGACCGTCACCCTCACCGTTGAACAGGCGGATAGCTGTTGTATCTTGAGCCTCAAATAGTTTTGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCg                                                                                      <  6:181728/139‑1 (MQ=255)
                                                    gccgACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTTATTTGCTGAGAAAAAATGCTTGGTGCATCTGTTGATGTTTTTTGCgg                                                                                     >  3:283124/1‑139 (MQ=255)
                                                       gACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTa                                                                                  <  2:314483/139‑1 (MQ=255)
                                                       gACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTa                                                                                  <  5:57703/139‑1 (MQ=255)
                                                             gTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTTATTTTATGAGGAAAAAAAGATTGTTCAACAGTTCATATATTTTGAGGGTAAgggaa                                                                            >  4:268671/1‑135 (MQ=255)
                                                                      aCCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTAtt                                                                                         >  8:272682/1‑117 (MQ=255)
                                                                      aCCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTAtt                                                                                         <  7:272682/117‑1 (MQ=255)
                                                                      aCCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCaa                                                                     <  4:353325/137‑1 (MQ=255)
                                                                      aCCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTTATTTGGTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTTTTGCGGGTAAGCGTCCAGCCCAcc                                                                   >  1:58913/1‑139 (MQ=255)
                                                                       ccGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAAccc                                                                  <  2:42493/139‑1 (MQ=255)
                                                                         gTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCt                                                                <  2:349866/139‑1 (MQ=255)
                                                                            gAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTTCTGTGAAAAAATGCTTGGTCAATATGTTAATGTTTATTTCGGGTAAGCGTCCAGCCACCCCCGTTg                                                             >  8:217979/1‑139 (MQ=255)
                                                                            gAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCCTTTTTGAATTTGATGAGAAAAAAAGCATCGTGAAAATGTTAATGTTAATTTCGGGTAAGTGTACAGCCAACCACTTTg                                                             >  8:72863/1‑139 (MQ=255)
                                                                                 ggCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTg                                                        <  7:217979/139‑1 (MQ=255)
                                                                                  gCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGGTTGTTTATTTGGTGAGAAAAAATTGGTGGGGAATCTGTTGATGTTTATTTCGGGTTAGCGGCCCGCCACCCCCGTTGgttttc                                                       >  4:12777/1‑137 (MQ=255)
                                                                                   cGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCGTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCa                                                      <  1:281132/139‑1 (MQ=255)
                                                                                   cGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCa                                                      <  2:202443/139‑1 (MQ=255)
                                                                                      aaaGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGc                                                   <  2:1682/139‑1 (MQ=255)
                                                                                      aaaGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGc                                                   <  6:363732/139‑1 (MQ=255)
                                                                                      aaaGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGCGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGc                                                   <  1:255148/139‑1 (MQ=255)
                                                                                            gttgtATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTTATTTTCTGAGAAAAAATGCTTGGTGAATCTGTTGATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGGAGGGCGTaa                                              >  3:38356/1‑138 (MQ=255)
                                                                                            gttgtATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTTATTTGGTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTTCGGGTAAGCGTCCAGCCAACCCCATGGTTTTGCAGGCCGTaa                                              >  6:182660/1‑138 (MQ=255)
                                                                                                  tCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTTAATTGCTGAGAAAAAATGCTTTGTGAATCTGTTCATGTTTATTTCGGGTAAGCGTCCAGCCAACCCCGTTGTTTTTCAGGCCGTAATAGcccc                                       >  3:330445/1‑139 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTTATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTTCGGGTTAGCGTCCAGCCAACCCCGTTGTTTTGCAGGCCTTAATAGCCCcgcgcg                                  >  7:226718/1‑139 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCa                                                      <  2:37669/119‑1 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCcgcgcg                                  >  3:81297/1‑139 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGGCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCa                                                      >  1:37669/1‑119 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTTATTTTGGGTGGAAAAAATGGTGGGGAATCTGTTGATGTTTTTTTAGGGTTAGCGGGCAGCCAAACCCGTTGTTTTTGAGGCCGTGATTGCCCcccgcc                                  >  7:326451/1‑138 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGTTTGTTTATTTTCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTTCGGGGAAGCGTCCAGCCAACCCCGTTGTTTTTCAGGCCGTAATATCCCcccgcg                                  >  5:225443/1‑139 (MQ=255)
                                                                                                       aGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCGGTTTTTTATTTTTTGAGAAAAAAATGGTGGTGGATAAGTTGATTTTTTTTTGGGGTAAGCGGCCAGGCCACCCCATTTTTTTGGAGGCCCTAAAAGCCCcccgcg                                  >  2:182572/1‑139 (MQ=255)
                                                                                                            cAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTGAATCTTTTCATTTTTATTTCGGGTAAGCGTCCAGCCCACCCCGTTGATTTTCAGGCCGTAAAAGCCCCCCCCGAGaa                              >  8:104918/1‑138 (MQ=255)
                                                                                                                      ggCGCGCGCGTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCCGCGCGAGAAACTcgccgccg                   <  2:197723/139‑1 (MQ=255)
                                                                                                                            cgcTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCCGCGCGAGAAACTCGCCGCCGctct               >  7:303311/1‑137 (MQ=255)
                                                                                                                                 gggCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCCGCGCGAGAAACTCGCCGCCGCTCTCAGaca         <  3:177653/138‑1 (MQ=255)
                                                                                                                                 gggCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCCGCGCGAGAAACTCGCCGCCGCTCTCAGaca         <  8:303311/138‑1 (MQ=255)
                                                                                                                                       cTTTTGTTAATTTGGTGAGTATAAATGCGTGGTCAATATGTTCAGGTTTATTGGGGGTTAGCGTGCAGCCCCCCCCGTTGTTTTGGAGGCCGTAATAGCCCCGCGCGCGGGCCTCGGCGCCGCTGTCAGCCACgcgatc  >  6:115860/1‑134 (MQ=255)
                                                                                                                                        |                                                                                                                                           
TGCTGTACCATTGAATCAGCAGGTAGCCGTCATAATAGTCGATGGTGACGCCGCCGACACCGTCACCCTCACCGTTGAACAGGCGGAAAGCTGTTGTATCTTGAGCCTCAAATAGTTTGGCGCGCGCTTGGGCCGCTTTTGTTAATTTGCTGAGAAAAAATGCTTGGTCAATCTGTTCATGTTTATTGCGGGTAAGCGTCCAGCCAACCCCTTTGTTTTGCAGGCCGTAATAGCCCCGCGCGAGAAACTCGCCGCCGCTCTCAGACACGATATC  >  CP006881/3754476‑3754749

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