Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,315 T→A 11.8% intergenic (+354/+260) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,3150TA11.8% 943.9 / 28.3 396intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (184/156);  new base A (35/18);  total (221/175)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 404 total aligned reads displayed.

CCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGG  >  CP009273/3263248‑3263382
                                                                   |                                                                   
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGT                                                                     <  1:21783158‑M1/68‑1 (MQ=255)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGT                                                                     <  1:7124943‑M1/68‑1 (MQ=255)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGT                                                                     <  1:15563103‑M1/68‑1 (MQ=255)
       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAAC                                                               >  1:38962413‑M1/1‑67 (MQ=255)
         cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCG                                                            <  1:18874499‑M1/68‑1 (MQ=255)
             ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGa                                                        <  1:7744542/68‑1 (MQ=35)
                 tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCGGCCCCTCCAAAAGAAAAAACCCGCCGAAgcg                                                    <  1:32535211‑M2/68‑4 (MQ=255)
                  ttGGGTATAGTGTCTTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgg                                                   >  1:36336843‑M2/1‑64 (MQ=255)
                   tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcggg                                                  >  1:8556685‑M2/1‑63 (MQ=255)
                    gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAACACCCGCCGAAgcgggt                                                 >  1:2977949‑M2/1‑62 (MQ=255)
                    gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggt                                                 >  1:25509811‑M2/1‑62 (MQ=255)
                    gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggt                                                 <  1:27155587‑M2/68‑7 (MQ=255)
                     ggTATAGTGTCGTGGACAGTCATTCATCTTTGTGCCCCTACAAAACAAAAAACCCGCCGAAgcgggtt                                                <  1:16692554‑M2/68‑8 (MQ=255)
                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggttt                                               <  1:36540530‑M2/68‑9 (MQ=255)
                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggttt                                               <  1:13773792‑M2/68‑9 (MQ=255)
                        ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAgcgggttttt                                             <  1:24640556‑M2/68‑11 (MQ=255)
                        ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAgcgggttttt                                             <  1:15974087‑M2/68‑11 (MQ=255)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggtttttac                                           >  1:37917368‑M2/1‑56 (MQ=255)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggtttttac                                           >  1:38623069‑M2/1‑56 (MQ=255)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggtttttac                                           <  1:4427188‑M2/68‑13 (MQ=255)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggtttttac                                           >  1:5142246‑M2/1‑56 (MQ=255)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAgcgggtttttac                                           >  1:9591503‑M2/1‑56 (MQ=255)
                           gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                          >  1:30758216/1‑68 (MQ=35)
                           gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACg                                          >  1:10125892/1‑68 (MQ=25)
                           gtgtCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                          >  1:23159900/1‑68 (MQ=25)
                           gtgtCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAAAAAAAAACCCGCCGAAgcggggttttaag                                          <  1:33519458‑M2/68‑14 (MQ=255)
                            tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                         >  1:14984572/1‑68 (MQ=35)
                             gtCGTGGACAGTCATTCATCTTTCTGACCCAACAAAAGTAAAAACCCGCCGCAGCGGCGTTTTAAGTa                                        >  1:22681853/1‑68 (MQ=11)
                             gtCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                        >  1:10575967/1‑68 (MQ=255)
                             gtCCTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                        >  1:2103157/1‑68 (MQ=25)
                               cGTGGTCAGAGATTCTTCTTTCTGCCCGTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:28250093/68‑1 (MQ=11)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:31425714/68‑1 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:17679075/68‑1 (MQ=35)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                   <  1:22957343/68‑1 (MQ=40)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                   >  1:6137760/1‑68 (MQ=40)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                   >  1:33898735/1‑68 (MQ=25)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACa                                   >  1:14860695/1‑68 (MQ=255)
                                  ggACAGTCATTCATCTTTCTGCCCATCCAAAAGTAAAAACCCGCCGCAGCGCGTTGTTACGTaaatcc                                   >  1:32320816/1‑65 (MQ=11)
                                    aCAGTCATTCATCTTTCTGTAAATCCAAAAGAAAAAACCAGCCGAAAACTCATTTTACGTAAAACAgg                                 >  1:27958341/1‑68 (MQ=2)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:25186082/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:25328292/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:14003414/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:37661175/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:5623276/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:12850822/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:12555861/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAgg                                 >  1:14325086/1‑68 (MQ=255)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCCGGTTTTTCCGTAAAACAgg                                 >  1:23965689/1‑68 (MQ=11)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCAGCCGAAGCCGGTTTTTACGTAAAACAgg                                 >  1:37150970/1‑68 (MQ=11)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAACAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:35762823/1‑68 (MQ=255)
                                    aCAGTCATTCATCTTTCTGCCCCACAAAAAGTAAAAAACCGCCGAACCGGCTTTTTACGTAAAAGAgg                                 >  1:4938940/1‑68 (MQ=14)
                                    aCAGTCATTCATCTCTCTGCCCCTCCAAAAGTAAACACCCGCCGCAGCGGGTTTTTACGTAAAACAg                                  <  1:6826242/67‑1 (MQ=14)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 <  1:28331275/63‑1 (MQ=255)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 <  1:19287342/63‑1 (MQ=255)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                              <  1:36702006/66‑1 (MQ=255)
                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                           <  1:35315617/68‑1 (MQ=255)
                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                           <  1:19695756/68‑1 (MQ=255)
                                          attcatCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                             >  1:8871069/1‑66 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                >  1:11899966/1‑59 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTgc                        >  1:26682733/1‑66 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:24234262/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:34253553/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:32846758/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:29177671/1‑67 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:11266559/68‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:33161958/68‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:17368610/68‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      >  1:34526981/1‑68 (MQ=40)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      >  1:6280867/1‑68 (MQ=40)
                                                 tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                    >  1:28649022/1‑68 (MQ=255)
                                                  ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                   >  1:9123435/1‑68 (MQ=40)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:31001960/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:18172831/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACAGATacg                >  1:18653805/1‑65 (MQ=255)
                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGca              >  1:16241106/1‑66 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:4404673/68‑1 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:13544067/68‑1 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:11409349/68‑1 (MQ=255)
                                                         ccTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:31885711/68‑1 (MQ=255)
                                                             cAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc        <  1:25949249/68‑1 (MQ=255)
                                                             cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttc        <  1:33428165/68‑1 (MQ=255)
                                                             aaaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc        <  1:7464067/67‑1 (MQ=255)
                                                              aaaaTTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:9294528/68‑1 (MQ=255)
                                                              aaaaTTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:11637613/68‑1 (MQ=255)
                                                              aaaaGTATAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttct       <  1:24098459/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:15039510/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:12232423/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:20226931/68‑1 (MQ=255)
                                                                aattAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:35478252/65‑1 (MQ=255)
                                                                aaGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:22208253/68‑1 (MQ=255)
                                                                aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:26186148/68‑1 (MQ=40)
                                                                aaGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:16921545/68‑1 (MQ=25)
                                                                 cttaaAAACCCGCttcggcgggtttttacttttggaggggcagaaagatgaatgactgtccacgacac    <  1:21515216‑M1/66‑56 (MQ=255)
                                                                 aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt    >  1:17069101/1‑68 (MQ=40)
                                                                  gTAAAAACCCGCttcggcgggtttttacttttggaggggcagaaagatgaatgactgtccacgacact   >  1:13223964‑M1/1‑12 (MQ=255)
                                                                  gTAAAAACCCGCttcggcgggtttttacttttggaggggcagaaagatgaatgactgtccacgacact   >  1:30038569‑M1/1‑12 (MQ=255)
                                                                  gTAAAAACCCGCttcggcgggtttttacttttggaggggcagaaagatgaatgactgtccacgacact   >  1:8133066‑M1/1‑12 (MQ=255)
                                                                  gTAAAAACCCGCttcggcgggtttttacttttggaggggcagaaagatgaatgactgtccacgacact   >  1:23364593‑M1/1‑12 (MQ=255)
                                                                  gTAAAAACCCGCttcggcgggtttttacttttggaggggcagaaagatgaatgactgtccacgacact   >  1:23093351‑M1/1‑12 (MQ=255)
                                                                   tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGCCCGATAAGCCGCTTTCTTTTgg  >  1:23690008/1‑68 (MQ=255)
                                                                   |                                                                   
CCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGG  >  CP009273/3263248‑3263382

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