Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,328 A→C 8.6% intergenic (+367/+247) 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,3280AC 8.6% 1560.9 / 23.8 694intergenic (+367/+247)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (229/344);  new base C (9/58);  total (239/455)
Fisher's exact test for biased strand distribution p-value = 7.90e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 698 total aligned reads displayed.

TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTG  >  CP009273/3263261‑3263395
                                                                   |                                                                   
ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGa                                                                     <  1:3224252/68‑1 (MQ=35)
  ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                                    >  1:43980634/1‑67 (MQ=35)
     ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                <  1:2825105/68‑1 (MQ=35)
      tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                               <  1:22584655/68‑1 (MQ=35)
       gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                              <  1:47491616/68‑1 (MQ=35)
       gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                              <  1:11572631/68‑1 (MQ=35)
       gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                              <  1:18216963/68‑1 (MQ=35)
         gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                            >  1:18942399/1‑68 (MQ=35)
         gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                            >  1:31134523/1‑68 (MQ=35)
         gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCCGGgtt                                                            >  1:3318711/1‑68 (MQ=11)
           ataGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGttttt                                                          >  1:8285111/1‑68 (MQ=255)
             aatgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:55648601/66‑1 (MQ=25)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:39134599/68‑1 (MQ=35)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:3798862/68‑1 (MQ=35)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:5622146/68‑1 (MQ=35)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:13139419/68‑1 (MQ=35)
              gtgtCGTGGACAGTCATTGATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                       >  1:44846596/1‑68 (MQ=255)
               tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGCAGCGGGTTTTTACGt                                                      >  1:56885879/1‑68 (MQ=255)
                gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                     >  1:39480096/1‑68 (MQ=35)
                  cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                   <  1:10024304/68‑1 (MQ=35)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                >  1:41821069/1‑68 (MQ=40)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGGAGAACa                                                >  1:33661305/1‑68 (MQ=255)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGATTTTTACGTAAAACa                                                >  1:19291580/1‑68 (MQ=255)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCACCGAAGCGGGTTTTTACGTtaaaaa                                                >  1:36466561/1‑66 (MQ=11)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                               <  1:52938750/67‑1 (MQ=40)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                               <  1:22044984/67‑1 (MQ=40)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                               <  1:5979534/67‑1 (MQ=40)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                              >  1:57975819/1‑68 (MQ=25)
                        cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGt                                             >  1:18433755/1‑68 (MQ=255)
                            cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                          <  1:5656191/67‑1 (MQ=40)
                            cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                          <  1:4350441/67‑1 (MQ=40)
                              ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                          >  1:23838254/1‑65 (MQ=255)
                                 atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                     >  1:24940276/1‑67 (MQ=255)
                                 atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                     >  1:30409755/1‑67 (MQ=255)
                                 atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                     >  1:4882178/1‑67 (MQ=255)
                                 atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                     >  1:53362632/1‑67 (MQ=255)
                                 atCTTTCTGACCATCCAAAAGTAAAAACCCGCATAAGCGGGTTTTTACGAAAAACACGTCAAACTGAc                                    >  1:24165441/1‑68 (MQ=14)
                                     ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                >  1:25333520/1‑68 (MQ=40)
                                       cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAACATGACCGaga                               >  1:32032330/1‑65 (MQ=14)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTATAACAGGTGAAACTGCACGAAAAg                             >  1:34481718/1‑67 (MQ=14)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                             >  1:6529227/1‑67 (MQ=255)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                             >  1:57286048/1‑67 (MQ=255)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                             >  1:37616203/1‑67 (MQ=255)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                             >  1:50017873/1‑67 (MQ=255)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                            <  1:30133993/68‑1 (MQ=255)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                            <  1:40925123/68‑1 (MQ=255)
                                         gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                            <  1:20295187/68‑1 (MQ=255)
                                                 aaaaGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                    <  1:35667434/68‑1 (MQ=255)
                                                 aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                    <  1:51801076/68‑1 (MQ=255)
                                                 aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                    <  1:55130963/68‑1 (MQ=255)
                                                 aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                    <  1:46518046/68‑1 (MQ=255)
                                                 aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                    <  1:44282022/68‑1 (MQ=255)
                                                   atgTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                  <  1:12039174/66‑1 (MQ=255)
                                                   aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                  <  1:49373876/68‑1 (MQ=40)
                                                   aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                  <  1:7749010/68‑1 (MQ=40)
                                                   aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                  <  1:42421435/68‑1 (MQ=40)
                                                    aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCtttt                 >  1:31729504/1‑68 (MQ=255)
                                                      tataaaCCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg               <  1:35089693/65‑1 (MQ=11)
                                                      tataaaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               <  1:46119004/65‑1 (MQ=255)
                                                      tAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               <  1:14434847/68‑1 (MQ=255)
                                                      cataaaCCCGCCGCATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               <  1:25797989/65‑1 (MQ=11)
                                                       acaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTggg              <  1:16993781/66‑1 (MQ=255)
                                                       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTggg              <  1:29228606/68‑1 (MQ=255)
                                                       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg              <  1:9634467/68‑1 (MQ=35)
                                                       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg              <  1:17331794/68‑1 (MQ=35)
                                                       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               <  1:42787141/67‑1 (MQ=35)
                                                        taaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:36102471/67‑1 (MQ=25)
                                                        caaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:45530529/67‑1 (MQ=255)
                                                        aaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:15424593/68‑1 (MQ=25)
                                                        aaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:15913169/68‑1 (MQ=25)
                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:12575095/68‑1 (MQ=35)
                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:23099117/68‑1 (MQ=35)
                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:48431971/68‑1 (MQ=35)
                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             <  1:55938108/68‑1 (MQ=35)
                                                          aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat           <  1:21589568/68‑1 (MQ=35)
                                                            cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGTTTCTTTTGGGTATAg         >  1:50542778/1‑68 (MQ=255)
                                                            cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg         >  1:11116261/1‑68 (MQ=35)
                                                             ccgccgAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt        <  1:47005441/68‑1 (MQ=255)
                                                             ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt        <  1:33950627/68‑1 (MQ=14)
                                                             ccgccgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt        <  1:2734454/68‑1 (MQ=255)
                                                             ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt        <  1:14693352/68‑1 (MQ=37)
                                                               gccgTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:13804440/68‑1 (MQ=14)
                                                               gccgCATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:52555733/68‑1 (MQ=14)
                                                               gccgCATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:40589236/68‑1 (MQ=14)
                                                               gccgCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:547460/68‑1 (MQ=255)
                                                               gccgCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:20842572/68‑1 (MQ=255)
                                                               gccgCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:6879783/68‑1 (MQ=255)
                                                               gccgCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:19962869/68‑1 (MQ=255)
                                                               gccgCACCGGTTTTTTACGTAAAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:38399868/68‑1 (MQ=11)
                                                               gccgCACCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:50981381/68‑1 (MQ=11)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:32699595/68‑1 (MQ=255)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAgtgt      <  1:21296137/68‑1 (MQ=255)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:9068556/68‑1 (MQ=40)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:106444/68‑1 (MQ=40)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:16191103/68‑1 (MQ=40)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:17834377/68‑1 (MQ=40)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:28941167/68‑1 (MQ=40)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGACGCTTTCTTTTGGGTATAgtgt      <  1:22810717/68‑1 (MQ=255)
                                                               gccgAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:41411392/68‑1 (MQ=255)
                                                               gccgAACCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt      <  1:26880011/68‑1 (MQ=255)
                                                                   |                                                                   
TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTG  >  CP009273/3263261‑3263395

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