Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,441 A→T 7.4% intergenic (+480/+134) 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,4410AT 7.4% 1168.7 / 12.9 470intergenic (+480/+134)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (255/142);  new base T (27/10);  total (312/158)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 476 total aligned reads displayed.

TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTG  >  CP009273/3263374‑3263508
                                                                   |                                                                   
 tctttTGGGTATAGTGTCGTGGACAGTCATTAATCTTTATGCCCCTCCAAAAGCAAAAACCCCCCGaa                                                                    >  1:16073211/1‑68 (MQ=11)
  ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGaa                                                                    >  1:29289623/1‑67 (MQ=25)
      tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                               <  1:1064321/68‑1 (MQ=35)
       gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                              <  1:33516576/68‑1 (MQ=35)
       gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                              <  1:43216186/68‑1 (MQ=35)
         gTATAGTGTCGTGGACAGTCATTCATCTTTGTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGgtt                                                            >  1:13714077/1‑68 (MQ=255)
         gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                             <  1:4228752/67‑1 (MQ=35)
         gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAACCAAAAACACGCCGAAGAGGGtgt                                                            >  1:17212407/1‑66 (MQ=11)
           ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCCGGgtttt                                                          >  1:8455729/1‑68 (MQ=11)
           ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                          >  1:23343507/1‑68 (MQ=35)
           ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCATACAAAAGCAAAAACCCGCCGCAGCGGGttttt                                                          >  1:39445433/1‑68 (MQ=11)
            taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                         <  1:8807377/68‑1 (MQ=35)
             aGTGTCTTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:35072711/68‑1 (MQ=255)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:37168398/68‑1 (MQ=35)
             aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                        <  1:14319399/68‑1 (MQ=35)
              gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                       >  1:5230945/1‑68 (MQ=33)
              gtgtCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTaag                                                       >  1:10556727/1‑66 (MQ=11)
              gtgtCGTGGACAGTCATTCATCTTTCTCCCCCAACAAAAGCAAAAACCCGCCGCAGCGGGGTTTTACg                                                       >  1:15626018/1‑68 (MQ=11)
               tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAATAGCCGCCGAAGCGGGTTTTTAc                                                        <  1:19842557/66‑1 (MQ=255)
               tgtCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGt                                                      >  1:28787190/1‑68 (MQ=255)
                 agaaagcggcttatcggtcagtttcacctggttttcgtTAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                    <  1:16580718‑M1/30‑1 (MQ=255)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACTTAAAcgt                                                >  1:37541549/1‑66 (MQ=9)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACTTAAAcgt                                                >  1:6098301/1‑66 (MQ=9)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                 <  1:30215212/67‑1 (MQ=33)
                     ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGATGGTTTTTACGTAAACCa                                                >  1:9772495/1‑68 (MQ=255)
                     ggACAGTCATTCATCTTTCTGACCATACAAAAGCAAAAACACGCAGAAGCGCCTTTTTACGTAAACCa                                                >  1:21939634/1‑68 (MQ=11)
                     ggACAGTCATTCATCTTTCTCCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                                                >  1:30904984/1‑68 (MQ=25)
                     ggACAGTCATTCATCTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAACCCa                                                >  1:12438738/1‑68 (MQ=255)
                     ggACAGTCATTCATCTTTCAGCCCCTACAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAACCa                                                >  1:21144731/1‑68 (MQ=11)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                               <  1:40587247/67‑1 (MQ=255)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGACGCGGGTTTTTACGTAAACCAgg                                              >  1:32756976/1‑68 (MQ=25)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                              >  1:39003149/1‑68 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                              >  1:7883881/1‑68 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                              >  1:25223266/1‑68 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                               <  1:36604105/67‑1 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                               <  1:33080605/67‑1 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                               <  1:23178883/67‑1 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                               <  1:33786594/67‑1 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                               <  1:24410587/67‑1 (MQ=35)
                       aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAgg                                              >  1:1919239/1‑68 (MQ=25)
                       aCAGTCATTCATCTTTCTCCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                              >  1:40861645/1‑68 (MQ=25)
                          gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa                                           >  1:42891626/1‑68 (MQ=35)
                            cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa                                          <  1:11749011/67‑1 (MQ=35)
                             attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa                                          >  1:26410529/1‑66 (MQ=35)
                               tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTg                                      <  1:38180116/68‑1 (MQ=35)
                                catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     <  1:38487906/68‑1 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:281557/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:37832663/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:31880317/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:34529362/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:29980810/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:11612392/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:11962345/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:21587734/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                     >  1:6460182/1‑67 (MQ=35)
                                 atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                    <  1:12622023/68‑1 (MQ=35)
                                  tcagtttcacctggtttacgtAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                   <  1:5736738‑M1/47‑1 (MQ=255)
                                    tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                 >  1:40219446/1‑68 (MQ=35)
                                    tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                 >  1:27416753/1‑68 (MQ=35)
                                     ttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAt                                >  1:42979999/1‑68 (MQ=35)
                                      tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                               >  1:26797687/1‑68 (MQ=35)
                                         gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                             >  1:26231144/1‑67 (MQ=35)
                                         gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                             >  1:13808725/1‑67 (MQ=35)
                                         gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                             >  1:11049186/1‑67 (MQ=35)
                                         gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                             >  1:18788872/1‑67 (MQ=35)
                                         gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                             >  1:35530254/1‑67 (MQ=35)
                                          ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                           >  1:17956468/1‑68 (MQ=35)
                                           cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                           >  1:15355142/1‑67 (MQ=35)
                                           cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                          >  1:19578134/1‑68 (MQ=35)
                                            acTCAAAAATCAAAAACCCGCCGTACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                         <  1:22979172/67‑1 (MQ=11)
                                                cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                     <  1:10013234/68‑1 (MQ=40)
                                                cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                       >  1:27005235/1‑66 (MQ=40)
                                                aaaaaTCATAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGctttc                     <  1:29476112/67‑1 (MQ=11)
                                                 aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                    <  1:15000124/68‑1 (MQ=255)
                                                 aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                    <  1:25626207/68‑1 (MQ=255)
                                                    aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGgtttctttg                 >  1:13332493/1‑67 (MQ=255)
                                                    aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                 >  1:20329935/1‑68 (MQ=40)
                                                     gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                >  1:8977880/1‑68 (MQ=39)
                                                     gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTg                >  1:20056823/1‑68 (MQ=255)
                                                     gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGACGCTTTCTTATg                >  1:4339306/1‑68 (MQ=255)
                                                      cttAAACCCGCCGTACCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:23848911‑M1/4‑68 (MQ=255)
                                                      ctcAAACCCGCCGTAGCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:34132467‑M1/4‑68 (MQ=255)
                                                      catAAACCCGCCGCATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:22648033‑M1/4‑68 (MQ=255)
                                                      cacAAACCCGCCGAACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:32443562‑M1/4‑68 (MQ=255)
                                                       tgAAACCCGCCGCATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg              >  1:4661193‑M1/3‑68 (MQ=255)
                                                       atAAACCCGCCGAATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:36163092‑M1/3‑67 (MQ=255)
                                                       AAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg              >  1:2057917‑M1/1‑68 (MQ=255)
                                                       AAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:17627991‑M1/1‑67 (MQ=255)
                                                       AAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg               >  1:7509361‑M1/1‑67 (MQ=255)
                                                        tAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:12821352‑M1/2‑68 (MQ=255)
                                                        tAAACCCGCCGTATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:16828780‑M1/2‑68 (MQ=255)
                                                        tAAACCCGCCGTAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:39868678‑M1/2‑68 (MQ=255)
                                                        tAAACCCGCCGAACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:2203764‑M1/2‑68 (MQ=255)
                                                        cAAACCCGCCGTATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:25063857‑M1/2‑68 (MQ=255)
                                                        AAAACCCGCCGTACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:37053095‑M1/1‑68 (MQ=255)
                                                        AAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt             >  1:31325350‑M1/1‑68 (MQ=255)
                                                          AACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat           <  1:41772180‑M1/68‑1 (MQ=255)
                                                             CCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc        >  1:24735908‑M1/1‑68 (MQ=255)
                                                             CCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc        >  1:15765492‑M1/1‑68 (MQ=255)
                                                               GCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGt      >  1:32185601‑M1/1‑68 (MQ=255)
                                                                   |                                                                   
TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTG  >  CP009273/3263374‑3263508

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: