Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,315 T→A 6.3% 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 6.3% 741.6 / 9.2 312intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (129/158);  new base A (23/2);  total (152/160)
Fisher's exact test for biased strand distribution p-value = 3.16e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

CCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGG  >  CP009273/3263248‑3263382
                                                                   |                                                                   
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:27578693/68‑1 (MQ=31)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:21321115/68‑1 (MQ=31)
ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                     <  1:15673345/68‑1 (MQ=31)
 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                    <  1:18474077/68‑1 (MQ=33)
 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                    <  1:3985434/68‑1 (MQ=33)
 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                    <  1:30073042/68‑1 (MQ=33)
 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                    <  1:26321557/68‑1 (MQ=33)
       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAAGTAAAAAc                                                               >  1:23234330/1‑67 (MQ=25)
                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                               >  1:11856401/1‑68 (MQ=35)
                      gTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAAAACGCCCAAGCGGGttt                                               >  1:11658145/1‑68 (MQ=11)
                        ataGTGTCGTGGACAGTCCTTCATCTTTCGGCCCCTACAAAAGTAAAAACCCGCCGCCGCGGGttttt                                             >  1:17947753/1‑68 (MQ=11)
                        ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAACTAAAAACCCGCCGCAGCGGGttttt                                             >  1:28773037/1‑68 (MQ=11)
                          tgtgtCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:26657508/67‑1 (MQ=255)
                          aGTGTCTTTGACATTCATTCATCTTTCTGCCCCTCCAAATGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:14442090/68‑1 (MQ=11)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:14738689/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:9934349/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:31020592/68‑1 (MQ=35)
                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                           <  1:24952014/68‑1 (MQ=35)
                           gtgtCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                          >  1:13021247/1‑68 (MQ=255)
                            tgtCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                         >  1:13511543/1‑68 (MQ=255)
                             gtCGTGGCCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGGTTTTACGTa                                        >  1:20625597/1‑68 (MQ=11)
                             gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACggt                                        >  1:25828960/1‑66 (MQ=25)
                             gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                        >  1:13185068/1‑68 (MQ=35)
                             gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                        >  1:4859904/1‑68 (MQ=35)
                             gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                        >  1:24134068/1‑68 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:6845216/68‑1 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:7917855/68‑1 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:1994301/68‑1 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      >  1:22306997/1‑68 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:18385774/68‑1 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:17394729/68‑1 (MQ=35)
                               cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                      <  1:33246746/68‑1 (MQ=35)
                                  ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                   >  1:9191594/1‑68 (MQ=40)
                                  ggACAGTCATTCATCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCCAGGCGGGTTTTTACGTAAAACa                                   >  1:33536061/1‑68 (MQ=14)
                                  ggACAGTCATTCATCTTTCTGCCCCTACAAAAGTAAAAACCCGCCGCAGCGGGGTTTTACGTAAAACa                                   >  1:22834541/1‑68 (MQ=14)
                                  ggACAGTCATTCATCTTTCTGCCCCTACAAAAGAAAAAACCCGCCGAAGAGGCTTGTTAAGTAAAACa                                   >  1:25232195/1‑68 (MQ=11)
                                  ggACAGTCATTCATCTTTATGCAACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                   >  1:19470000/1‑68 (MQ=14)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAgg                                 >  1:10643946/1‑68 (MQ=255)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAgg                                 >  1:10790749/1‑68 (MQ=255)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGAGGGTTGTTACGTAAAACAgg                                 >  1:30694709/1‑68 (MQ=14)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACTgg                                 >  1:5908430/1‑68 (MQ=255)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:16473659/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:29578258/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:19517855/1‑68 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:27240399/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:28189205/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:28141013/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:20999650/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                  <  1:12423126/67‑1 (MQ=40)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:8180358/1‑68 (MQ=25)
                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAAACCGAAGAAGCGCGTTTATACGTAAAACAgg                                 >  1:20219645/1‑68 (MQ=11)
                                    aCAGTCATTCATCTTTCGGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                 >  1:11298320/1‑68 (MQ=255)
                                       gTCATTCATCTTTCTGCCCCTACAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                              >  1:2529815/1‑68 (MQ=255)
                                         cattcatCTTTCTGCCCCTCCAAAATTAAAAAGCCGCCGAAGCGTGTTTTTACGTAAAACAGGTGaa                             <  1:24319350/67‑1 (MQ=14)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                <  1:34484792/64‑1 (MQ=40)
                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                             >  1:10194225/1‑67 (MQ=40)
                                         catgcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                             >  1:16154137/1‑67 (MQ=255)
                                             cttCTTTCTGGCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        <  1:12212706/66‑1 (MQ=14)
                                             catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        <  1:5163161/68‑1 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:4302349/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:12591021/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:7213226/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:19352389/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:1979505/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:17824772/1‑67 (MQ=255)
                                              atCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                        >  1:2166683/1‑67 (MQ=255)
                                               tatttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:27002787/66‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:27904740/68‑1 (MQ=255)
                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                      <  1:32991268/68‑1 (MQ=255)
                                                cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                     >  1:9681130/1‑68 (MQ=255)
                                                cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                     >  1:2186437/1‑68 (MQ=255)
                                                 tttCTGCCCCTCCAAAAGTAAAAACCCGCCGCATCGGGTTTTTTCGTAAAACAGGTGAAACTGACCGa                    <  1:9590687/68‑1 (MQ=14)
                                                 tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                    <  1:15067191/68‑1 (MQ=40)
                                                   tCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                  >  1:5528286/1‑68 (MQ=255)
                                                   tCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                  >  1:7589812/1‑68 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:8408529/1‑67 (MQ=255)
                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                >  1:1669458/1‑67 (MQ=255)
                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCg             >  1:17212391/1‑68 (MQ=40)
                                                         ccTCCAAAAGTTAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:24431642/68‑1 (MQ=255)
                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:32085305/68‑1 (MQ=255)
                                                         acTCCAAAAGTTAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc            <  1:32126225/67‑1 (MQ=255)
                                                           tcaaaaaGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctt          <  1:33453022/65‑1 (MQ=255)
                                                              aaaaTTATAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:8781037/68‑1 (MQ=14)
                                                              aaaaTTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttct       <  1:31706739/68‑1 (MQ=14)
                                                              aaaaTTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:22557282/68‑1 (MQ=255)
                                                              aaaaGTAAAACCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:34004205/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:10372030/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAATCGGGTTTTTACGTAAATGAGGTGAAACTGACCTATAAGCCGctttct       <  1:26460525/68‑1 (MQ=11)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:29307709/68‑1 (MQ=255)
                                                              aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct       <  1:2345883/68‑1 (MQ=255)
                                                               aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt      <  1:17554117/68‑1 (MQ=40)
                                                                aaGTATAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt     <  1:16824886/68‑1 (MQ=14)
                                                                aaGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt     <  1:20721719/68‑1 (MQ=14)
                                                                aaGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:13853838/68‑1 (MQ=255)
                                                                aaGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:22726644/68‑1 (MQ=255)
                                                                aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt     <  1:24739647/68‑1 (MQ=255)
                                                                aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt     <  1:6538437/68‑1 (MQ=40)
                                                                 aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGggttctttg    >  1:32506924/1‑67 (MQ=11)
                                                                 aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcgttctgtg    >  1:4521103/1‑65 (MQ=255)
                                                                  gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg   >  1:28448501/1‑68 (MQ=39)
                                                                   |                                                                   
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: