Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,330 G→T 6.7% 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.7% 1527.0 / 6.0 696intergenic (+369/+245)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (232/358);  new base T (5/65);  total (240/456)
Fisher's exact test for biased strand distribution p-value = 9.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 697 total aligned reads displayed.

TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGG  >  CP009273/3263265‑3263396
                                                                 |                                                                  
  tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                <  1:22584655/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)
     gTATAGTGTCGTGGACAGTCATTAATATTTATGCCCCAACAAAATTAAAAACCCGCCCCAGCGGGttt                                                             >  1:43170837/1‑68 (MQ=11)
       ataGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGttttt                                                           >  1:8285111/1‑68 (MQ=255)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:23765607/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:13139419/68‑1 (MQ=35)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                         <  1:5622146/68‑1 (MQ=35)
           tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAACCGCCGCAGCGGGTTTTTACGt                                                       >  1:56885879/1‑68 (MQ=255)
            gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                      >  1:46234138/1‑68 (MQ=35)
              cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                    <  1:10024304/68‑1 (MQ=35)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGATTTTTACGTAAAACa                                                 >  1:19291580/1‑68 (MQ=255)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCACCGAAGCGGGTTTTTACGTtaaaaa                                                 >  1:36466561/1‑66 (MQ=11)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGCGTTTTTACGTAAAACa                                                 >  1:4843475/1‑68 (MQ=255)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCTGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:33874850/67‑1 (MQ=255)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:46612149/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:5979534/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:4257475/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:29174740/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                <  1:22044984/67‑1 (MQ=40)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAgg                                               >  1:35001984/1‑68 (MQ=255)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:49458764/1‑68 (MQ=25)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:4059200/1‑68 (MQ=25)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:50100883/1‑68 (MQ=25)
                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                               >  1:57975819/1‑68 (MQ=25)
                   aCAGTCATTCATCTTTCTGAAAATCCAAAAGAAAAAACCCGCCGAAGAGGGTTTTTACGTAAAACAgg                                               >  1:27248203/1‑68 (MQ=11)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGt                                              >  1:18433755/1‑68 (MQ=255)
                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                           <  1:5656191/67‑1 (MQ=40)
                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                            <  1:40963830/66‑1 (MQ=255)
                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                            <  1:47239990/66‑1 (MQ=255)
                        cattcatCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAATACAGGTGaa                                           >  1:33038226/1‑67 (MQ=11)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                    <  1:44101346/68‑1 (MQ=255)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                    <  1:52993568/68‑1 (MQ=255)
                                 ttCTGCCACTCAAAAAGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                 <  1:38084870/68‑1 (MQ=14)
                                    tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                              <  1:25791395/68‑1 (MQ=255)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                              >  1:6529227/1‑67 (MQ=255)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                              >  1:50547622/1‑67 (MQ=255)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                              >  1:57286048/1‑67 (MQ=255)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                             <  1:20295187/68‑1 (MQ=255)
                                        ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                          >  1:51041620/1‑68 (MQ=255)
                                             aaaaTTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:28564252/68‑1 (MQ=255)
                                             aaaaGTTAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:45539224/68‑1 (MQ=255)
                                             aaaaGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:551798/68‑1 (MQ=255)
                                             aaaaGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:35667434/68‑1 (MQ=255)
                                             aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:31955346/68‑1 (MQ=255)
                                             aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:51834422/68‑1 (MQ=255)
                                             aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                     <  1:26450465/68‑1 (MQ=255)
                                               atgTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:12039174/66‑1 (MQ=255)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGCCCGATAAGCCGctttctgt                   >  1:29967262/1‑66 (MQ=255)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:28176007/68‑1 (MQ=40)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                   <  1:7749010/68‑1 (MQ=40)
                                                aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCtttt                  >  1:25293361/1‑68 (MQ=255)
                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:43632614/68‑1 (MQ=37)
                                                  tAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:14434847/68‑1 (MQ=255)
                                                  cacaaaCCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:30341082/65‑1 (MQ=11)
                                                   ataaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:44868127/65‑1 (MQ=11)
                                                   ataaaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:52610386/65‑1 (MQ=255)
                                                   acaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:16993781/66‑1 (MQ=255)
                                                   aaaatCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:33545968/67‑1 (MQ=255)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:9634467/68‑1 (MQ=35)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg               <  1:17331794/68‑1 (MQ=35)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:5478471/67‑1 (MQ=35)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                <  1:32211071/67‑1 (MQ=35)
                                                    taaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:37646982/67‑1 (MQ=11)
                                                    taaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:36102471/67‑1 (MQ=25)
                                                    aaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:41588028/68‑1 (MQ=25)
                                                    aaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:15424593/68‑1 (MQ=25)
                                                    aaaaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:53453713/68‑1 (MQ=25)
                                                    aaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:38453135/68‑1 (MQ=25)
                                                    aaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:15913169/68‑1 (MQ=25)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:12575095/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:32614180/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:31626956/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:47530182/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:24759132/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:23099117/68‑1 (MQ=35)
                                                    aaaaCCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt              <  1:39172067/68‑1 (MQ=25)
                                                      aaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat            <  1:34433902/68‑1 (MQ=25)
                                                      aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat            <  1:21589568/68‑1 (MQ=35)
                                                        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGCCCGATAAGCCGCTTTCGTTTGGGTATAg          >  1:5570713/1‑68 (MQ=255)
                                                        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg          >  1:11116261/1‑68 (MQ=35)
                                                         ccgccgAATCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt         <  1:48888239/68‑1 (MQ=255)
                                                         ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt         <  1:14693352/68‑1 (MQ=37)
                                                           gccgTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:4182555/68‑1 (MQ=14)
                                                           gccgTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:13804440/68‑1 (MQ=14)
                                                           gccgGAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:24102094/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)
                                                           gccgCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:3113282/68‑1 (MQ=255)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAgtgt       <  1:21296137/68‑1 (MQ=255)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:17834377/68‑1 (MQ=40)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:16191103/68‑1 (MQ=40)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:39716608/68‑1 (MQ=40)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:9068556/68‑1 (MQ=40)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:106444/68‑1 (MQ=40)
                                                           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGACGCTTTCTTTTGGGTATAgtgt       <  1:22810717/68‑1 (MQ=255)
                                                           gccgAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:55954192/68‑1 (MQ=255)
                                                           gccgAACCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt       <  1:44315300/68‑1 (MQ=255)
                                                                 |                                                                  
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGG  >  CP009273/3263265‑3263396

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: