Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,750,454 G→T 9.9% W30L (TGG→TTG)  PP_1550 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,750,4540GT 9.9% 155.1 / 5.2 50W30L (TGG→TTG) PP_1550Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base T (2/3);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGCCGTCCAGCTGTTCAG  >  NC_002947/1750310‑1750603
                                                                                                                                                |                                                                                                                                                     
gCTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCgg                                                                                                                                                  >  2:781119/1‑150 (MQ=255)
 cTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGa                                                                                                                                                 <  1:690077/150‑1 (MQ=255)
  ttGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCATCCGCCCATTGGCGGac                                                                                                                                                <  1:428485/150‑2 (MQ=255)
         gAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCAGGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTg                                                                                                                                         >  2:772436/1‑150 (MQ=255)
          aaTATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGa                                                                                                                                        <  1:947188/150‑1 (MQ=255)
          aaTATGCAAGCATGCTTGAATGCCGGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGa                                                                                                                                        <  2:1089272/150‑1 (MQ=255)
                   gCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCa                                                                                                                               >  1:200839/1‑150 (MQ=255)
                     aTGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAGGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGt                                                                                                                             <  1:897007/150‑1 (MQ=255)
                       gCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTAc                                                                                                                           >  2:932763/1‑150 (MQ=255)
                               gCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGa                                                                                                                   <  1:781119/150‑1 (MQ=255)
                                  tGGACGAGCAGAGGAGGTCCCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCg                                                                                                                >  1:343563/1‑150 (MQ=255)
                                       gAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGa                                                                                                           >  1:109643/1‑150 (MQ=255)
                                       gAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGa                                                                                                           <  2:1020362/150‑1 (MQ=255)
                                             aggaggTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAgg                                                                                                     >  1:897465/1‑150 (MQ=255)
                                             aggagcTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAgg                                                                                                     <  1:175953/150‑1 (MQ=255)
                                                  gTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGgcgcg                                                                                                >  1:448041/1‑150 (MQ=255)
                                                     gCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTg                                                                                             >  1:1169245/1‑150 (MQ=255)
                                                       cATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGaa                                                                                           >  2:928242/1‑150 (MQ=255)
                                                       cATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGaa                                                                                           <  2:448041/150‑1 (MQ=255)
                                                       cATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGaa                                                                                           <  2:1098172/150‑1 (MQ=255)
                                                             aaggagCCAGGCACTCCACCATGGCGAGTCTATTTCGGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGCGCGTTTGACCAAGGGCGCGCTGAAGGCAgg                                                                                     <  2:672166/147‑1 (MQ=255)
                                                                  gCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGa                                                                                                                                                 >  2:872113/1‑85 (MQ=255)
                                                                  gCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGa                                                                                                                                                 <  1:872113/85‑1 (MQ=255)
                                                                   ccAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAg                                                                               >  1:897336/1‑150 (MQ=255)
                                                                    cAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCt                                                                                                                                          >  1:72266/1‑90 (MQ=255)
                                                                    cAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCt                                                                                                                                          <  2:72266/90‑1 (MQ=255)
                                                                       gCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGAt                                                                           <  2:109643/150‑1 (MQ=255)
                                                                           tCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCg                                                                       >  1:203359/1‑150 (MQ=255)
                                                                                           tATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAg                                                       >  2:1039798/1‑150 (MQ=255)
                                                                                              ttCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGAAGCCGCTGCGGTACCTTGTAGaaa                                                    <  2:897336/150‑1 (MQ=255)
                                                                                                  gCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAAtttt                                                <  1:932763/150‑1 (MQ=255)
                                                                                                   cccTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTc                                               >  2:1056892/1‑150 (MQ=255)
                                                                                                    ccTGCCGAAGGCGCTTGGGGTGGCCTATGAGGCCGTCCGCCAATGGGCGGATGTTCTTGAGCTACGACACGACGAAATCGAGCGGTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTcg                                              <  2:343563/150‑2 (MQ=255)
                                                                                                        gCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTa                                          >  2:885546/1‑150 (MQ=255)
                                                                                                        gCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTCTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATGTTCACCTa                                          >  1:990795/1‑150 (MQ=255)
                                                                                                          gAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATc                                        <  2:832774/150‑1 (MQ=255)
                                                                                                              gcgcTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGc                                    <  2:521015/150‑1 (MQ=255)
                                                                                                                     gggTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTgcgc                             <  2:283702/150‑1 (MQ=255)
                                                                                                                       gTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCAc                                             <  1:1126022/132‑1 (MQ=255)
                                                                                                                       gTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCAc                                             >  2:1126022/1‑132 (MQ=255)
                                                                                                                        tGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCag                          <  1:1039798/150‑1 (MQ=255)
                                                                                                                            cTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGCATAAATTTTCACCTATCTGGCCTTGCGCCagtagt                      >  2:885590/1‑150 (MQ=255)
                                                                                                                             tATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTa                                                        <  2:504018/115‑1 (MQ=255)
                                                                                                                             tATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTa                                                        >  1:504018/1‑115 (MQ=255)
                                                                                                                                 aGGCCGTCCGCCAATGGGCGGATGTTCGCTAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGc                 >  1:635457/1‑150 (MQ=255)
                                                                                                                                     catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                                >  1:788297/3‑131 (MQ=255)
                                                                                                                                     catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                                <  2:788408/129‑1 (MQ=255)
                                                                                                                                     catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                                <  2:788297/129‑1 (MQ=255)
                                                                                                                                     catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                                >  1:788408/3‑131 (MQ=255)
                                                                                                                                               tGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGCCGTCCAGCTGTTCa   <  2:200839/150‑1 (MQ=255)
                                                                                                                                                gggCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGCCGTCCAGCTGTTCAg  <  2:990795/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGCCGTCCAGCTGTTCAG  >  NC_002947/1750310‑1750603

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