Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,196,188 G→A 10.0% A145T (GCC→ACC)  PP_5618 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,196,1880GA10.0% 159.0 / 4.7 57A145T (GCC→ACC) PP_5618Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (2/4);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAACTCGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGA  >  NC_002947/4196043‑4196334
                                                                                                                                                 |                                                                                                                                                  
gTTAACTCGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAg                                                                                                                                                <  1:37074/150‑1 (MQ=255)
       cGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTgg                                                                                                                                         <  1:1070272/150‑1 (MQ=255)
              ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAggg                                                                                                                                              <  2:818618/138‑1 (MQ=255)
              ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAggg                                                                                                                                              >  1:818618/1‑138 (MQ=255)
              ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGt                                                                                                                                  >  2:1053127/1‑150 (MQ=255)
              ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGt                                                                                                                                  <  2:608216/150‑1 (MQ=255)
                    tGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCa                                                                                                                            >  2:736674/1‑150 (MQ=255)
                     gggAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCAc                                                                                                                           <  1:1102979/150‑1 (MQ=255)
                     gggAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGTagacttt                                                                                                                           <  2:1133373/150‑7 (MQ=255)
                           aCAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAGGCCACTGCTCg                                                                                                                     >  2:290190/1‑150 (MQ=255)
                             aGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCt                                                                                                                   >  1:419631/1‑150 (MQ=255)
                               cAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTc                                                                                                                 <  2:286723/150‑1 (MQ=255)
                               cAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTc                                                                                                                 <  2:712563/150‑1 (MQ=255)
                               cAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTc                                                                                                                 >  1:712563/1‑150 (MQ=255)
                                           gCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCCAACTGTTTCGTTGCCATTGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCga                                                                                                     <  1:1053127/150‑1 (MQ=255)
                                             tGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGt                                                                                                                                        >  2:1186207/1‑113 (MQ=255)
                                             tGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGt                                                                                                                                        <  1:1186207/113‑1 (MQ=255)
                                                 aGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAggg                                                                                               >  1:633464/1‑150 (MQ=255)
                                                   ggaggaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTc                                                                                                                      >  2:972582/1‑125 (MQ=255)
                                                   ggaggaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTc                                                                                                                      <  1:972582/125‑1 (MQ=255)
                                                      ggaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATg                                                                                          <  2:663983/150‑1 (MQ=255)
                                                       gaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTggcaacg                                                                                                                                    <  2:365436/107‑6 (MQ=255)
                                                       gaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTggcaacg                                                                                                                                    >  1:365436/1‑102 (MQ=255)
                                                         tCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGcc                                                                                       <  1:546162/150‑1 (MQ=255)
                                                            aaaaGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGa                                                                                    >  1:443526/1‑150 (MQ=255)
                                                                   gCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTc                                                                             <  2:419631/150‑1 (MQ=255)
                                                                      cttCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCa                                                                          >  2:15832/2‑150 (MQ=255)
                                                                       ttCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCAt                                                                         <  1:86424/150‑1 (MQ=255)
                                                                          aaGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATgag                                                                      <  2:332694/150‑1 (MQ=255)
                                                                           aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAc                                                                                                         >  2:684611/1‑114 (MQ=255)
                                                                           aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAc                                                                                                         <  1:684611/114‑1 (MQ=255)
                                                                           aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGTa                                                                                                                                 >  2:97650/1‑90 (MQ=38)
                                                                           aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGTa                                                                                                                                 <  1:97650/90‑1 (MQ=37)
                                                                              aTGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAcgc                                                                                                       >  1:365237/1‑113 (MQ=255)
                                                                              aTGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAcgc                                                                                                       <  2:365237/113‑1 (MQ=255)
                                                                               ccgatctcTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTCATCGCCCGACAAGCTCGCATGAGATCGc                                                                 <  1:1068436/143‑1 (MQ=255)
                                                                                gCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGcc                                                                >  2:428875/1‑150 (MQ=255)
                                                                                       tAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTAGCTTCTTCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGccagatcgg                                                         >  2:1068436/1‑143 (MQ=255)
                                                                                             ccTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGc                                                   <  1:512375/150‑1 (MQ=255)
                                                                                                 aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                               <  1:905884/150‑1 (MQ=255)
                                                                                                 aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                               >  1:896221/1‑150 (MQ=255)
                                                                                                 aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                               <  1:736674/150‑1 (MQ=255)
                                                                                                 aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                               >  2:905884/1‑150 (MQ=255)
                                                                                                         agagaCAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGt                                       >  1:672806/1‑150 (MQ=255)
                                                                                                             aCAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGccc                                   >  2:687323/1‑150 (MQ=255)
                                                                                                               aGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAgg                                              >  2:293315/1‑137 (MQ=255)
                                                                                                               aGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAgg                                              <  1:293315/137‑1 (MQ=255)
                                                                                                                gCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTAt                                <  1:15832/150‑1 (MQ=255)
                                                                                                                 ccGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTa                                                                                                          >  1:1170232/1‑75 (MQ=255)
                                                                                                                 ccGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTa                                                                                                          <  2:1170232/75‑1 (MQ=255)
                                                                                                                  cGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATga                                                                       <  2:842373/109‑1 (MQ=255)
                                                                                                                  cGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATga                                                                       >  1:842373/1‑109 (MQ=255)
                                                                                                                              gcTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCg                  >  2:450426/1‑150 (MQ=255)
                                                                                                                               cTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGt                 <  2:494620/150‑1 (MQ=255)
                                                                                                                                        tCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGc        >  1:216841/1‑150 (MQ=255)
                                                                                                                                          tACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCg      >  1:668604/1‑150 (MQ=255)
                                                                                                                                              gTTGCCAGGGCCTGGTGACGGTAGGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGa  <  1:290190/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTTAACTCGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGA  >  NC_002947/4196043‑4196334

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