Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,750,451 A→C 10.0% Q29P (CAA→CCA)  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,4510AC10.0% 127.9 / 8.0 50Q29P (CAA→CCA) PP_1550Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base C (2/3);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

AAGCATGCTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGC  >  NC_002947/1750304‑1750588
                                                                                                                                                   |                                                                                                                                         
aaGCATGCTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAAt                                                                                                                                         >  1:1089272/1‑150 (MQ=255)
      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:1098172/150‑1 (MQ=255)
                                                             cATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGaa                                                                            <  2:448041/150‑1 (MQ=255)
                                                                   aaggagCCAGGCACTCCACCATGGCGAGTCTATTTCGGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGCGCGTTTGACCAAGGGCGCGCTGAAGGCAgg                                                                      <  2:672166/147‑1 (MQ=255)
                                                                        gCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGa                                                                                                                                  <  1:872113/85‑1 (MQ=255)
                                                                        gCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGa                                                                                                                                  >  2:872113/1‑85 (MQ=255)
                                                                         ccAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAg                                                                >  1:897336/1‑150 (MQ=255)
                                                                          cAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCt                                                                                                                           <  2:72266/90‑1 (MQ=255)
                                                                          cAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCt                                                                                                                           >  1:72266/1‑90 (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                              >  2:1126022/1‑132 (MQ=255)
                                                                                                                             gTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCAc                              <  1:1126022/132‑1 (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                 <  2:788408/129‑1 (MQ=255)
                                                                                                                                           catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                 <  2:788297/129‑1 (MQ=255)
                                                                                                                                           catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                 >  1:788408/3‑131 (MQ=255)
                                                                                                                                           catCCGCCCATTGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTg                 >  1:788297/3‑131 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
AAGCATGCTTGTTTTGAATATGCAAGCATGCTTGAATGCCTGGACGAGCAGAGGAGGTCGCCATGACCAAGAGCCAGGCAATCCAACATTTCGGGTCTATTTCAGCCCTGGCGAAGGCGCTAGGGGTGACCTATGAGGCCGTCCGCCAATGGGCGGATGTTCCTGAGCTTCGCCAGTACCAAATCGAGCGTTTGACCAAGGGCGCGCTGAAGGCAGGGAAAAGGGATGCCGCTGCGTGACCTTGTAGAAATTTTCACCTATCTGGCCTTGCGCCAGTAGTGGGGC  >  NC_002947/1750304‑1750588

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: ‑