Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,315 T→A 8.5% 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,3150TA 8.5% 1074.4 / 12.0 440intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (175/216);  new base A (43/4);  total (218/222)
Fisher's exact test for biased strand distribution p-value = 1.88e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 461 total aligned reads displayed.

CCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGG  >  CP009273/3263248‑3263382
                                                                   |                                                                   
ccGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGTCAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:50143631/68‑1 (MQ=255)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:40769148/68‑1 (MQ=31)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:20628864/68‑1 (MQ=31)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:12655901/68‑1 (MQ=31)
 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                    <  1:14636700/68‑1 (MQ=33)
               ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                       >  1:43980634/1‑67 (MQ=35)
                  ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:8423001/50‑1 (MQ=32)
                   tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                  <  1:20976506/68‑1 (MQ=35)
                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                               >  1:47308463/1‑68 (MQ=35)
                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGttt                                               >  1:53318597/1‑68 (MQ=25)
                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCCGGgtt                                               >  1:3318711/1‑68 (MQ=11)
                      gTATAGTGTCGTGGACAGTCATTCATCATTTTTCCAATACAAAACTAAAAACCCCCCGGACCggggat                                               >  1:21437063/1‑65 (MQ=2)
                      gTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                               >  1:49015282/1‑68 (MQ=255)
                        ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGgtttt                                             >  1:17123256/1‑68 (MQ=25)
                          aatgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:55648601/66‑1 (MQ=25)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:56061256/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:3798862/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:39134599/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:43486175/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:11746536/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:43236761/68‑1 (MQ=35)
                           gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAACCGGGTTTTTACg                                          >  1:18314396/1‑68 (MQ=25)
                           gtgtCGTGGACAGTCATTCATCTTTCTGACCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                          >  1:12403596/1‑68 (MQ=255)
                             gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                        >  1:17446364/1‑68 (MQ=35)
                               cGTGGCCAGTCATTCATCTTTCTGCCCCGCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTaaa                                      >  1:9707427/1‑68 (MQ=11)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:10024304/68‑1 (MQ=35)
                                  ggACAGTCATTGATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                   >  1:25577033/1‑68 (MQ=255)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGGAGAACa                                   >  1:33661305/1‑68 (MQ=255)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCACCGAAGCGGGTTTTTACGTtaaaaa                                   >  1:36466561/1‑66 (MQ=11)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACa                                   >  1:48048765/1‑68 (MQ=255)
                                   gACAGTCATTCATCTTTCTGCCCCTCCTAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:35982157/68‑1 (MQ=255)
                                   gACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  >  1:18569045/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:28315990/67‑1 (MQ=255)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:24473373/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:6753458/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:1572179/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:22044984/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:46612149/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:2713812/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:29099837/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:52938750/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:57638667/1‑68 (MQ=25)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:49458764/1‑68 (MQ=25)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:37223035/1‑68 (MQ=25)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:4059200/1‑68 (MQ=25)
                                    aCAGTCATTCATCTTTCTGACCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTATTTACGTAAAACTgg                                 >  1:7075922/1‑68 (MQ=11)
                                    aCAGTCATTCATCTATCTGACAATCAAAAAGAAAAAACCCCCCGAAGCGTGTTTTTACGTAAAAAAgg                                 >  1:37760165/1‑68 (MQ=2)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 <  1:20086621/63‑1 (MQ=255)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 <  1:23469805/63‑1 (MQ=255)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                <  1:13044744/64‑1 (MQ=40)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                             >  1:1899498/1‑67 (MQ=40)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                             <  1:31268842/67‑1 (MQ=40)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                              <  1:22905310/66‑1 (MQ=255)
                                         cattcatCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                             >  1:5963954/1‑67 (MQ=25)
                                            tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTg                         >  1:31870747/1‑68 (MQ=40)
                                             catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        <  1:16444788/68‑1 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTGAAACAGGTGAAACTGa                        >  1:56351310/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:24940276/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:45650243/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:21705597/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:19459246/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:30409755/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:8679854/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:10480416/1‑67 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:42150084/68‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:16099623/68‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAca                      >  1:22288250/1‑67 (MQ=255)
                                                cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                     <  1:11245331/68‑1 (MQ=255)
                                                cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                     >  1:789041/1‑68 (MQ=255)
                                                cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                     >  1:55007055/1‑68 (MQ=255)
                                                  ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                   >  1:35489901/1‑68 (MQ=40)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTATAACAGGTGAAACTGCACGAAAAg                >  1:34481718/1‑67 (MQ=14)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTTAGCGATAAg                >  1:24014312/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:42653956/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:11022436/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:57107962/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:32374566/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:46396722/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc               <  1:30133993/68‑1 (MQ=255)
                                                      gCCCCTCCAAAAGAAAAAACCCGCCGAAGCGCGTTTTTACGTAAAACAGGTGAAACAGACCCATAAg                >  1:3995519/1‑67 (MQ=11)
                                                      gCCCCTCAAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc               <  1:41681002/68‑1 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:51209038/68‑1 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:13674589/68‑1 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:38383719/68‑1 (MQ=255)
                                                         acTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:6335221/67‑1 (MQ=255)
                                                          cTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGct           <  1:25268944/68‑1 (MQ=255)
                                                           tCCAAAATGAAAACCCCGCGGCGCCTGGTTTATACGTAACACAGGTTAAACTGACCGATAAGCCgggt          <  1:5672650/68‑4 (MQ=0)
                                                              aaaaGTATAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttct       <  1:29775203/68‑1 (MQ=14)
                                                              aaaaGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:551798/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:34908690/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:44282022/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:44757612/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:22624203/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:26450465/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:51986239/68‑1 (MQ=255)
                                                                aattAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt     <  1:28116455/65‑1 (MQ=255)
                                                                aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:18098245/68‑1 (MQ=40)
                                                                 aGTAAAAACCCGCCGAAGCTGGGTGTTAAGAAAAACACGTGGAACTCGCCGCTTAGCCGGTTTCtgtt    >  1:50575841/1‑68 (MQ=2)
                                                                  ttAAAAACCCGCCGTACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTg   <  1:33847126/67‑1 (MQ=11)
                                                                  gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGgtttctttgg   >  1:32115377/1‑66 (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: