Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,188,519 A→T 9.2% intergenic (+115/‑219) PP_3684 → / → PP_5611 hypothetical protein/transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,188,5190AT 9.2% 141.5 / 2.6 51intergenic (+115/‑219)PP_3684/PP_5611hypothetical protein/transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base T (2/3);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAGTGGTGATGTGTGCAAAA  >  NC_002947/4188374‑4188667
                                                                                                                                                 |                                                                                                                                                    
cccGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACa                                                                                                                                                   <  5:271570/149‑1 (MQ=255)
   ggCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGtt                                                                                                                                                >  1:220662/1‑149 (MQ=255)
   ggCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTAGTGGTTTCAGATATCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTATGTGCACTTTGCGCTCGACAACTCCGTATAGACATGAAATGAGATATGTCGCAAAAAGtt                                                                                                                                                >  1:220667/1‑149 (MQ=255)
        aTCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCg                                                                                                                                           >  4:33853/1‑149 (MQ=255)
            gAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCATCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTAGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAg                                                                                                                                       >  8:201657/1‑149 (MQ=255)
                        aaCCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAAt                                                                                                                            >  3:187773/1‑148 (MQ=255)
                         aCCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGa                                                                                                                          <  2:150897/149‑1 (MQ=255)
                          ccTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGaa                                                                                                                         <  4:187773/149‑1 (MQ=255)
                            tGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAata                                                                                                                       >  7:59019/1‑149 (MQ=255)
                                 gctgctTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATatcatc                                                                                                                  <  8:263556/149‑1 (MQ=255)
                                           gTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCTAACAGTTGTCCGCTAGGCTCAAGCAATGAATATCATCGCAGACCAt                                                                                                         <  3:33853/148‑1 (MQ=255)
                                              tCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAAt                                                                                                      >  2:318586/1‑148 (MQ=255)
                                                 gATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAgg                                                                                                  <  7:99875/149‑1 (MQ=255)
                                                  aTCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAgg                                                                                                  >  7:194311/1‑148 (MQ=255)
                                                  aTCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACAATAATGAgg                                                                                                  >  5:159051/1‑148 (MQ=255)
                                                       aaGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACaaa                                                                                            <  7:213705/149‑1 (MQ=255)
                                                                   cTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACt                                                                                <  7:48323/149‑1 (MQ=255)
                                                                       aGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATatca                                                                                                                    >  3:66446/1‑109 (MQ=255)
                                                                       aGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATatca                                                                                                                    <  4:66446/109‑1 (MQ=255)
                                                                       aGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGttt                                                                            >  8:159595/1‑149 (MQ=255)
                                                                         cGACTTGAACTGCGTTTTTTGTGCCCTGTGCGCTCTACAACTCCGTTTATAGAAGTTCTGTGATATGGCGCATATTGTTGTCCGCTAGCCTCAAGCAATGAATATCATCTCAGACCATTATGAGGCACAAACTCTTCAAAACTGTTTgc                                                                          <  2:220667/149‑1 (MQ=255)
                                                                         cGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTgc                                                                          <  2:220662/149‑1 (MQ=255)
                                                                             ttCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCaaaa                                                                                  >  7:109725/1‑137 (MQ=255)
                                                                             ttCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCaaaa                                                                                  <  8:109725/137‑1 (MQ=255)
                                                                             ttCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTgagcaa                                                                      >  3:10134/1‑147 (MQ=255)
                                                                                       tATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTATTTGCGCCACGCTAATCaa                                                            >  3:77167/1‑149 (MQ=255)
                                                                                        aTTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTTAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCaa                                                            >  4:216619/1‑148 (MQ=255)
                                                                                            tgttCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATAcc                                                        <  1:83391/148‑1 (MQ=255)
                                                                                            tgtgCCCTTTGCGCTCGACAACTCCGTTTAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCATCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATAcc                                                        <  6:182094/148‑1 (MQ=255)
                                                                                            tgtgCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCAACGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATAcc                                                        <  1:6630/148‑1 (MQ=255)
                                                                                               gCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGa                                                     <  7:159595/148‑1 (MQ=255)
                                                                                                      gagatCGACAACTCCGTATAGATAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTa                                              <  8:59019/144‑1 (MQ=255)
                                                                                                         cTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGATGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGc                                          <  7:201657/149‑1 (MQ=255)
                                                                                                          tCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGc                                          >  5:102556/1‑148 (MQ=255)
                                                                                                              cAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAg                                     >  3:253262/1‑149 (MQ=255)
                                                                                                              cAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAg                                     >  3:248931/1‑149 (MQ=255)
                                                                                                                 cTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTAtt                                            <  1:255059/139‑1 (MQ=255)
                                                                                                                 cTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTAtt                                            >  2:255059/1‑139 (MQ=255)
                                                                                                                  tCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTTAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTa                                 <  3:216619/149‑1 (MQ=255)
                                                                                                                            cAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACtt                        <  1:318586/148‑1 (MQ=255)
                                                                                                                             aaGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAATCTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTg                       >  4:208891/1‑148 (MQ=255)
                                                                                                                               gTACTGAGATATGGCGCATACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTaa                    >  2:7718/1‑149 (MQ=255)
                                                                                                                               gTACTGAGATATGGCGCATACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTAATGCGAAAGTTTACAAATACTTGTaa                    >  2:2458/1‑149 (MQ=255)
                                                                                                                                 aCTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACAATAATTAGGCACAAAATCTTCAAAAATGTTTGCGCCAAGATAATCAATAACGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAgt                  >  4:22873/1‑149 (MQ=255)
                                                                                                                                  cTGAGATATGGCGCATACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAgtg                 <  1:7718/149‑1 (MQ=255)
                                                                                                                                  cTGAGATATGGCGCATACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAgtg                 <  1:2458/149‑1 (MQ=255)
                                                                                                                                  cTGAGATATGGAGCAAACAGTTGTCCGATAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAgtg                 <  4:77167/149‑1 (MQ=255)
                                                                                                                                         atGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGa                                                                                                                          >  7:6907/1‑37 (MQ=255)
                                                                                                                                         atGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGa                                                                                                                          <  8:6907/37‑1 (MQ=255)
                                                                                                                                            gcgcAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTg                                           >  2:210974/1‑113 (MQ=255)
                                                                                                                                            gcgcAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAGTGGTGAtgtgtg       >  5:75463/1‑149 (MQ=255)
                                                                                                                                                aaaCAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAGTGGTGATGTGTGCaa    >  5:137161/1‑148 (MQ=255)
                                                                                                                                                 aaCAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAGTGGTGATGTGTGCaaaa  >  1:26153/1‑149 (MQ=255)
                                                                                                                                                 aaCAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTAAAAACTTTTTGCGCCAATCTAAACAATACAGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAGTTTTGATCTGTGAaaaa  >  3:248476/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAGTACTGAGATATGGCGCAAACAGTTGTCCGCTAGCCTCAAGCAATGAATATCATCGCAGACCATAATGAGGCACAAACTCTTCAAAACTGTTTGCGCCACGCTAATCAATACCGGATAATTTATTGCGAAAGTTTACAAATACTTGTAAGTGGTGATGTGTGCAAAA  >  NC_002947/4188374‑4188667

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑