Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,330 G→T 6.8% intergenic (+369/+245) 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,3300GT 6.8% 1215.6 / 3.6 539intergenic (+369/+245)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (178/272);  new base T (4/49);  total (186/353)
Fisher's exact test for biased strand distribution p-value = 6.95e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 539 total aligned reads displayed.

TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGG  >  CP009273/3263265‑3263396
                                                                 |                                                                  
 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                 <  1:36919199/68‑1 (MQ=35)
     gTATAGTGTCGTGGGCAGTCATTCATCTTTGTGCCCCAACAAAAGTAAAAACCCGCCGAACCGGGttt                                                             >  1:11901075/1‑68 (MQ=11)
     gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                             >  1:13933568/1‑68 (MQ=35)
     gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCCGGttt                                                             >  1:37595829/1‑68 (MQ=25)
     gTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                             >  1:25242938/1‑68 (MQ=255)
       ataGTGTCGTGGACAGTCATTCATCATTAGTCCGATACAAAAGTAAAAACCCGCCGAAGCGGGgtttt                                                           >  1:38091638/1‑68 (MQ=11)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:27695372/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:2041723/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:7557510/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:22390626/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:22736697/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         >  1:43332764/1‑68 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:15504523/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:39117234/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:13180501/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:10914076/68‑1 (MQ=35)
          gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                        >  1:6202918/1‑68 (MQ=35)
          gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                        >  1:3494285/1‑68 (MQ=35)
          gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                        >  1:11480055/1‑68 (MQ=35)
            gtCGTGGACAGTCATTCATCTTTGTGCCCCTACAAAAGTAAAAAAACTACGAATCGGGTTTTTAAGTa                                                      >  1:27037563/1‑68 (MQ=11)
            gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                      >  1:6811422/1‑68 (MQ=35)
            gtCGTGGACAGTCATTCATCGTTCTGCCCCTACAAAACTAAAAACCCCCCCAACCGGGGTTTTACGtg                                                      >  1:5580833/1‑67 (MQ=11)
              cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                    >  1:8472277/1‑68 (MQ=35)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                 >  1:2793058/1‑68 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:35287693/1‑68 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:31204061/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:1456117/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:26705047/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCAGAAGCGCGTTTTTACGTAAAACAgg                                               >  1:1241522/1‑68 (MQ=11)
                   aCAGTCATTCATCTTTCTGCCCCTACAAAAGAAAAAACCCGCCGAAGAGGGTTTTTACGTAAAAAAgg                                               >  1:16482212/1‑68 (MQ=11)
                   aCAGTCATTCATCTTTCTGCCACTAAAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:26172195/1‑68 (MQ=11)
                   aCAGTCATTCATCTTTCTCCGCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:30524923/1‑68 (MQ=11)
                      gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                            >  1:41677755/1‑68 (MQ=255)
                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                           <  1:7201348/67‑1 (MQ=40)
                        cattcatCTCTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                           >  1:42650384/1‑67 (MQ=255)
                          ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                        <  1:18516566/68‑1 (MQ=40)
                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                      >  1:23544345/1‑67 (MQ=255)
                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                      >  1:37895429/1‑67 (MQ=255)
                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                      >  1:22219317/1‑67 (MQ=255)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                    <  1:36492666/68‑1 (MQ=255)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                    <  1:31731769/68‑1 (MQ=255)
                                 ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                 <  1:15817554/68‑1 (MQ=255)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                              >  1:33725390/1‑67 (MQ=255)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                              >  1:19791274/1‑67 (MQ=255)
                                       cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCg                           >  1:16149299/1‑68 (MQ=40)
                                        ccTCCAAAAGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                          <  1:32565465/68‑1 (MQ=255)
                                        ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGcc                            >  1:21515637/1‑66 (MQ=255)
                                             aaaaTTATAAACCCGCCGAATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttct                     <  1:42358041/68‑1 (MQ=14)
                                             aaaaTTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:17344489/68‑1 (MQ=255)
                                             aaaaGTTAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGctttct                     <  1:30893767/68‑1 (MQ=11)
                                             aaaaGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:39401052/68‑1 (MQ=255)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:28792369/68‑1 (MQ=40)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:36207146/68‑1 (MQ=40)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:1362480/68‑1 (MQ=40)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:4441160/68‑1 (MQ=40)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:21831168/68‑1 (MQ=255)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:9305409/68‑1 (MQ=37)
                                                  ctcaaaCCCGCCGCATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:20119147/65‑1 (MQ=11)
                                                   ataaaGCCGCCGAATCGGGTTTTTAAGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:15152157/66‑1 (MQ=11)
                                                   acaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:28564822/65‑1 (MQ=255)
                                                   acaaaCCCGCCGGATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTg                 <  1:21016910/64‑1 (MQ=11)
                                                   aaaaaCCCGCCGTAGCGGGTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:19255131/67‑1 (MQ=11)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:38682904/68‑1 (MQ=35)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:11358576/68‑1 (MQ=35)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                 <  1:40923469/66‑1 (MQ=35)
                                                   aaaaaCCCGCCGAACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:10432344/68‑1 (MQ=11)
                                                    taaaCCCGCCGTATCTGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:35756577/67‑1 (MQ=11)
                                                    taaaCCCGCCGTACCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:40297763/67‑1 (MQ=11)
                                                    aaaaCCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:24495032/68‑1 (MQ=11)
                                                    aaaaCCCGCCGCATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:33235880/68‑1 (MQ=11)
                                                    aaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:29770560/68‑1 (MQ=255)
                                                    aaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:32287105/68‑1 (MQ=255)
                                                    aaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:34142992/68‑1 (MQ=25)
                                                    aaaaCCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:34666570/68‑1 (MQ=25)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:35480926/68‑1 (MQ=25)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:33498320/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:1818376/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:38367404/67‑1 (MQ=11)
                                                      aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat            >  1:27389048/1‑68 (MQ=35)
                                                      aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat            <  1:5319883/68‑1 (MQ=35)
                                                        cccGCCGAATCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg          <  1:17086173/68‑1 (MQ=255)
                                                        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGTTTCTTTTGTGTATAg          >  1:30024873/1‑68 (MQ=255)
                                                        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg          >  1:25858411/1‑68 (MQ=35)
                                                        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAATGACCGATAAGCCGCATTCTTTTGGGTAAAg          >  1:7841870/1‑68 (MQ=11)
                                                         ccgccgTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt         <  1:13445550/68‑1 (MQ=11)
                                                           gccgCAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:9409003/68‑1 (MQ=11)
                                                           gccgAATCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:10055398/68‑1 (MQ=255)
                                                           gccgAATCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:43022642/68‑1 (MQ=255)
                                                           gccgAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:30243355/68‑1 (MQ=255)
                                                           gccgAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:41149122/68‑1 (MQ=255)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:24061123/68‑1 (MQ=255)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:22946516/68‑1 (MQ=40)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:17738261/68‑1 (MQ=40)
                                                           gccgAACCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:18894615/68‑1 (MQ=255)
                                                           gccgAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:8794836/68‑1 (MQ=255)
                                                           gccgAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:40469081/68‑1 (MQ=255)
                                                              gtaGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt    <  1:29193021/66‑1 (MQ=255)
                                                              gcaGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt    <  1:24854848/66‑1 (MQ=255)
                                                              gAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt    <  1:6438359/68‑1 (MQ=255)
                                                                aGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTgg  <  1:25600060/68‑1 (MQ=40)
                                                                 |                                                                  
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGG  >  CP009273/3263265‑3263396

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: