Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,328 A→C 8.8% 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.8% 1040.3 / 13.7 475intergenic (+367/+247)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (154/228);  new base C (10/36);  total (165/310)
Fisher's exact test for biased strand distribution p-value = 1.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 475 total aligned reads displayed.

CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCG  >  CP009273/3263263‑3263393
                                                                 |                                                                 
    tGGGTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                             <  1:13762890/68‑1 (MQ=25)
     gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                            <  1:8806816/68‑1 (MQ=35)
     gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                            <  1:13420708/68‑1 (MQ=35)
       gTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                          >  1:31853320/1‑68 (MQ=255)
           aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                      <  1:9934349/68‑1 (MQ=35)
           aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                      <  1:4047674/68‑1 (MQ=35)
           aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                      <  1:31020592/68‑1 (MQ=35)
           aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                      <  1:22270747/68‑1 (MQ=35)
           aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                      <  1:19903304/68‑1 (MQ=35)
            gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                     >  1:24751388/1‑68 (MQ=35)
              gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                   >  1:10339054/1‑68 (MQ=35)
              gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                   <  1:28184400/68‑1 (MQ=35)
              gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                   >  1:27998988/1‑68 (MQ=35)
              gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCCCCGAAGCGGGTTTTTACGTa                                                   >  1:16877694/1‑68 (MQ=25)
                cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                 <  1:33246746/68‑1 (MQ=35)
                cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                 <  1:6845216/68‑1 (MQ=35)
                cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                 <  1:23986131/68‑1 (MQ=35)
                cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                 <  1:30559926/68‑1 (MQ=35)
                   ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                >  1:10584127/1‑66 (MQ=37)
                   ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                              <  1:28465664/68‑1 (MQ=40)
                   ggACAGTCATTCATCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCCAGGCGGGTTTTTACGTAAAACa                                              >  1:33536061/1‑68 (MQ=14)
                   ggACAGTCATTCATCTTTCTGCCCAAACAAAAGAAAAAACCCGCCGAAGCGGGGTTTTACGTaaaaac                                              >  1:10694583/1‑66 (MQ=11)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:5107838/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:20135521/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:33945540/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:12423126/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:21617963/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:11940645/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:25119846/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                             <  1:2431219/67‑1 (MQ=40)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCAGCCGAAGCGGGTTTTTACGTAAAACAgg                                            >  1:3088584/1‑68 (MQ=255)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAgg                                            >  1:21157491/1‑68 (MQ=255)
                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                            >  1:9323263/1‑68 (MQ=25)
                     aCAGTCATTCATCTTTCTGCCCCTACAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                            >  1:17408795/1‑68 (MQ=255)
                     aCAGTCATTCATCTTTCTCACCCTCCAAAAGAAAAAACCCGCCGCATCCCGGTTTTAAGTAAAACacg                                            >  1:15662060/1‑66 (MQ=2)
                          ccatcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                            <  1:17591524/60‑1 (MQ=255)
                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                            <  1:7196652/63‑1 (MQ=255)
                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                            <  1:29362689/63‑1 (MQ=255)
                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                           <  1:34484792/64‑1 (MQ=40)
                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                           <  1:13139639/64‑1 (MQ=40)
                          cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                         <  1:4629974/66‑1 (MQ=255)
                          catgcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                        >  1:16154137/1‑67 (MQ=255)
                            ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                     <  1:17759377/68‑1 (MQ=40)
                              cttCTTTCTGGCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                   <  1:12212706/66‑1 (MQ=14)
                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGTGGGTTTTTACGTAAAACAGGTGAAACTGa                                   >  1:19377392/1‑67 (MQ=255)
                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                   >  1:20647159/1‑67 (MQ=255)
                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                   >  1:26342161/1‑67 (MQ=255)
                               atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                   >  1:2303395/1‑67 (MQ=255)
                                tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                 <  1:25281244/68‑1 (MQ=255)
                                tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                 >  1:4425933/1‑68 (MQ=40)
                                       gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                           >  1:9599421/1‑67 (MQ=255)
                                       gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                           >  1:8323588/1‑67 (MQ=255)
                                         cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCg                        >  1:17212391/1‑68 (MQ=40)
                                          acTCCAAAAGTTAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                       <  1:32126225/67‑1 (MQ=255)
                                          acTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                       <  1:32877846/67‑1 (MQ=255)
                                            tcaaaaaGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctt                     <  1:33453022/65‑1 (MQ=255)
                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcttt                    <  1:2444756/68‑1 (MQ=255)
                                              cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc                   <  1:10946559/68‑1 (MQ=40)
                                               aaaaTTATAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                  <  1:2574447/68‑1 (MQ=255)
                                               aaaaTTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                  <  1:22557282/68‑1 (MQ=255)
                                               aaaaGTTCAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttct                  <  1:22741835/68‑1 (MQ=14)
                                               aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                  <  1:26604335/68‑1 (MQ=255)
                                               aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                  <  1:2345883/68‑1 (MQ=255)
                                                 aattAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt                <  1:14330753/65‑1 (MQ=255)
                                                 aattAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                <  1:7507761/65‑1 (MQ=255)
                                                 aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                <  1:6538437/68‑1 (MQ=40)
                                                  aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctttg               >  1:8050591/1‑67 (MQ=255)
                                                  aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt               >  1:29971980/1‑68 (MQ=40)
                                                    tataaaCCCGCCGCATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg             <  1:1485168/65‑1 (MQ=11)
                                                    tAAAAACCCGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg             <  1:27479244/68‑1 (MQ=11)
                                                    cataaaCCCGCCGCACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTgg             <  1:11201206/65‑1 (MQ=11)
                                                    attaaaCCCGCCGGATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg             <  1:14024218/65‑1 (MQ=11)
                                                     acaaaCCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg            <  1:5472035/66‑1 (MQ=255)
                                                     aaaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTggg            <  1:3859270/68‑1 (MQ=11)
                                                     aaaaaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg            <  1:26888990/68‑1 (MQ=25)
                                                     aaaaaCCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg             <  1:10019611/67‑1 (MQ=25)
                                                      taaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:12772907/67‑1 (MQ=25)
                                                      taaaCCCGCCGAACCGGTTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:16613005/67‑1 (MQ=11)
                                                      aaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:31544908/68‑1 (MQ=25)
                                                      aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:8731282/68‑1 (MQ=35)
                                                      aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:18339244/68‑1 (MQ=35)
                                                      aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:18626078/68‑1 (MQ=35)
                                                      aaaaCCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt           <  1:19751538/68‑1 (MQ=25)
                                                        aaCCCGCCGAACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGgat         <  1:27087141/68‑1 (MQ=11)
                                                         aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata        >  1:15225610/1‑68 (MQ=35)
                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTGTGGGGTagag       >  1:14547894/1‑68 (MQ=11)
                                                           ccgccgTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt      <  1:29157560/68‑1 (MQ=11)
                                                           ccgccgCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt      <  1:6164077/68‑1 (MQ=255)
                                                           ccgccgAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt      <  1:4872569/68‑1 (MQ=255)
                                                           ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt      <  1:778133/68‑1 (MQ=14)
                                                             gccgTATCGGTTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:12590033/68‑1 (MQ=14)
                                                             gccgTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:16409183/68‑1 (MQ=11)
                                                             gccgTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:21941096/68‑1 (MQ=11)
                                                             gccgTAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:28850338/68‑1 (MQ=255)
                                                             gccgTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:11679947/68‑1 (MQ=255)
                                                             gccgCAGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:3255558/68‑1 (MQ=255)
                                                             gccgCACCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:20816257/68‑1 (MQ=11)
                                                             gccgAAGCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:2035408/68‑1 (MQ=255)
                                                             gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:17903439/68‑1 (MQ=40)
                                                             gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt    <  1:19504921/68‑1 (MQ=40)
                                                                 |                                                                 
CTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCG  >  CP009273/3263263‑3263393

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: