Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,466 G→T 17.5% intergenic (+505/+109) 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,4660GT17.5% 967.4 / 37.2 416intergenic (+505/+109)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (209/108);  new base T (26/49);  total (251/165)
Fisher's exact test for biased strand distribution p-value = 9.58e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

AGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCT  >  CP009273/3263399‑3263533
                                                                   |                                                                   
aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTg                                                                     >  1:43613024/1‑68 (MQ=35)
 gTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGCAGCGCCTTTATACGTAAACCAGGTgt                                                                    >  1:38310503/1‑67 (MQ=11)
  ttatcggtcagtttcacctgatttacgtaaaaacccgcttcgGCGGGTTTTTACGTAAATCAGGTGaa                                                                   >  1:2570392‑M1/43‑68 (MQ=255)
   cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                  <  1:44940/68‑1 (MQ=35)
   cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa                                                                   <  1:40028615/67‑1 (MQ=35)
   cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa                                                                   <  1:17783837/67‑1 (MQ=35)
   cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa                                                                   >  1:56814493/1‑67 (MQ=35)
       catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              <  1:19217346/68‑1 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              >  1:29730089/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              >  1:23529060/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              >  1:44873461/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              >  1:57761425/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              >  1:43136697/1‑67 (MQ=35)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACAAGGTGAAACTGa                                                              >  1:55938553/1‑67 (MQ=25)
        atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGTAATCCACGTCAAACAGa                                                              >  1:18418187/1‑67 (MQ=11)
        atCTTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                              >  1:41266014/1‑67 (MQ=25)
         tCTTTCTGCCCCTCCACAAGCAAAACCCCGCCGAAGAGGGTTTTTAAGTAAACCAGGTGAAACTGAcc                                                            <  1:45378962/68‑1 (MQ=11)
         tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                            >  1:14397874/1‑68 (MQ=35)
         tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                            <  1:30607687/68‑1 (MQ=35)
            ttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAt                                                         >  1:6340006/1‑68 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                      >  1:4991593/1‑67 (MQ=35)
                gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAg                                                      >  1:40832203/1‑67 (MQ=35)
                 ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGCCCGATAAGcc                                                    >  1:10422905/1‑68 (MQ=25)
                 ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                    >  1:46702868/1‑68 (MQ=35)
                 ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                    >  1:27532599/1‑68 (MQ=35)
                  cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                   >  1:26287815/1‑68 (MQ=35)
                  cccTCCAAAAGCAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                    >  1:50979837/1‑66 (MQ=255)
                   ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                    >  1:42736156/1‑66 (MQ=35)
                   ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                  >  1:30125425/1‑68 (MQ=37)
                   ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                  <  1:34193175/68‑1 (MQ=37)
                       cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                >  1:18765995/1‑66 (MQ=40)
                       aaaaaTCACAAACCCGCCGGACCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttc                                              <  1:7089504/67‑1 (MQ=14)
                        aaaaTCATAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGctttct                                             <  1:38863518/68‑1 (MQ=11)
                        aaaaTCACAAACCCGCCGAACCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATACGCCGctttct                                             <  1:57267439/68‑1 (MQ=11)
                        aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                             <  1:20358162/68‑1 (MQ=255)
                        aaaaGCAAAAACCCGCCGAACCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                             <  1:1691793/68‑1 (MQ=255)
                         aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttctt                                            >  1:16199866/1‑68 (MQ=255)
                          aatcAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                           <  1:41840374/65‑1 (MQ=255)
                          aaGCATAAACCCGCCGCAGCGGGTTATTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttcttt                                           <  1:38163388/68‑1 (MQ=14)
                          aaGCACAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttcttt                                           <  1:22534636/68‑1 (MQ=14)
                          aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGcgttcttt                                           >  1:55267020/1‑68 (MQ=255)
                           aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                          >  1:52633316/1‑68 (MQ=40)
                            gtaaaaaCCCGCCGAAGCGGGTTTTTACGTAATCCAGGGGAAACTGACCGATAAGCCGGTGTGTTTTg                                         >  1:47508907/3‑68 (MQ=11)
                            gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                         >  1:23274616/1‑68 (MQ=39)
                            gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                         >  1:39063830/1‑68 (MQ=39)
                             cataaaCCCGCCGCACCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTgg                                        <  1:40480453/65‑1 (MQ=11)
                             cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                        >  1:48063015/1‑68 (MQ=37)
                               taaaCCCGCCGCATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                      <  1:11571747/67‑1 (MQ=11)
                               caaaCCCGCCGAACCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                      <  1:49584433/67‑1 (MQ=11)
                               aaaaCCCGCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                      <  1:44993455/68‑1 (MQ=25)
                               aaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                      <  1:9375276/68‑1 (MQ=35)
                               aaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                      <  1:28334022/68‑1 (MQ=35)
                                 aaCCCGCCGCAGCGGGTTATTACGTAAACCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                    <  1:32605139/68‑1 (MQ=11)
                                 aaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                    <  1:21221340/68‑1 (MQ=35)
                                  aCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                   >  1:18169546/1‑68 (MQ=35)
                                  aCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                   >  1:42390566/1‑68 (MQ=35)
                                  aCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGtata                                   >  1:24115421/1‑68 (MQ=255)
                                   cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                   >  1:55006593/1‑67 (MQ=35)
                                   cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                  >  1:13675717/1‑68 (MQ=35)
                                   cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGtata                                   >  1:6736402/1‑67 (MQ=25)
                                   cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGtata                                   >  1:32203126/1‑67 (MQ=255)
                                    ccgccgTAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc                                 <  1:35392125/68‑1 (MQ=255)
                                    ccgccgCATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGc                                 <  1:54571302/68‑1 (MQ=255)
                                      gccgTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:15212111/68‑1 (MQ=255)
                                      gccgCAGCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:26753842/68‑1 (MQ=14)
                                      gccgATGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:12583521/68‑1 (MQ=255)
                                      gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:45826339/68‑1 (MQ=40)
                                      gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:4881649/68‑1 (MQ=40)
                                      gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:15637363/68‑1 (MQ=40)
                                      gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:8729918/68‑1 (MQ=40)
                                      gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:31310828/68‑1 (MQ=40)
                                      gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                               <  1:1209448/68‑1 (MQ=40)
                                       ccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                >  1:6978472/1‑66 (MQ=39)
                                             CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGgaa                        <  1:48446390‑M1/68‑2 (MQ=255)
                                              GGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                       >  1:39720756‑M1/1‑68 (MQ=255)
                                                 ttattACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                    <  1:13019119/65‑1 (MQ=1)
                                                 TTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                     >  1:19928197‑M1/1‑67 (MQ=255)
                                                 TTTTTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:14082981‑M1/1‑68 (MQ=255)
                                                 TTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                     >  1:28578587‑M1/1‑67 (MQ=255)
                                                 TTTTTACGTAAACCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGt                     >  1:1063421‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:37415532‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCa                   >  1:53076956‑M1/1‑68 (MQ=255)
                                                  TTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:29666562‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:33166381‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCa                   >  1:33561966‑M1/1‑68 (MQ=255)
                                                  TTTTACGTAAACCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:7627609‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:51987438‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:10854478‑M1/1‑67 (MQ=255)
                                                  TTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTc                    >  1:14760867‑M1/1‑67 (MQ=255)
                                                      aCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGgcattca               >  1:53814837/1‑68 (MQ=255)
                                                          aaaTCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCt            <  1:20761180/67‑1 (MQ=255)
                                                          aaaCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCt            <  1:50470846/67‑1 (MQ=40)
                                                             ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATcat           >  1:41380191/1‑63 (MQ=37)
                                                             ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATGtt           >  1:35926631/1‑65 (MQ=255)
                                                             ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTAGTGGAAAGTCATTAATATTTca        >  1:24707348/1‑67 (MQ=11)
                                                                 gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGCGTCGTGGCCAGTCATTCATATTTCTTcc     >  1:3618614/1‑67 (MQ=11)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGTCATTCATCTTTCAGccc    >  1:37188724/1‑67 (MQ=255)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATGTTTCAGcccc   >  1:27792720/1‑68 (MQ=255)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGcc     >  1:35590809/1‑66 (MQ=35)
                                                                  tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCCTTCATCTTTAAGcccc   >  1:10015371/1‑68 (MQ=11)
                                                                   |                                                                   
AGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCT  >  CP009273/3263399‑3263533

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: