Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,285,744:1 +G 100% coding (47/1014 nt) rssB → PcnB‑degradosome interaction factor; response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,285,7441.G90.9% 57.2 / 2.1 22coding (47/1014 nt)rssBPcnB‑degradosome interaction factor; response regulator
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (7/13);  total (8/14)
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.41e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATTCAACCTGTTTGCCCGCAAATATCTACGCTTGATTTCCATCGCGCGCACGCTGCCATTGCGGCCGGACAGCTGGCAGTGGAAAGGAAAATGGACGAACTTTTGCCGTTGGTACGCACCAACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAA  >  CP009273/1285477‑1285984
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
gATTCAACCTGTTTGCCCGCAAATATCTACGCTTGATTTCCATCGCGAGCACGCTGCCATTGCGGCCGGACAGCTGGCAGTGGAAAGGAAAATGGACGAACTTTTGCCGTTGGTACGCACCAACATTTGACCAGACTTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCCCTATTGAGTAAAGCCAGTCAGGGGAGAGAAACTGCAGCAGCCATTGGTCGGAAAACAGATTCACATTGTTGAAAATGAggcaagt                                                                                                                                                                                                                                           >  1:129155/1‑269 (MQ=255)
                                                  gcGCTGCCATTGCGGCCGGACAGCTGGCAGTGGAAAGGAAAATGGAAGAACTTTTGCCGTTGGTACGCACAAACATTTGACCAGAAGTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACg                                                                                                                                                                                           <  1:252980/273‑1 (MQ=255)
                                                                           ttcAGTCGAAAGGAAAAGTGTCAGACTTGTGCCTTTGGTAAGCACCAACATTTTTCCAATATTTTTATCTACACTTAAGGTAATTCTTACAGGTGCAGGTGGCAATAGCAGGCCACTATTGAGTAAAGCCAGTCAAGGGAGAGAACATGACGCCGCCATTGGTCGGAAAACAGATTCTCATTGTTGCAGATGAGGCAGGTATGTCGGTCGCGTCGGGCATAATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGa                                                                                                                                                                    <  2:129149/270‑1 (MQ=255)
                                                                                                                  aCGCACCAACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGATCGGAGAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTgatgat                                                                                                                          <  1:434758/275‑1 (MQ=255)
                                                                                                                           cATTTGACCAGAATTTTTATCTACACTTAAGTTACTTCTGACAGTCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAACCAGATTCTCATTGTTGAAGATGCGGCCGGTATTtcgctcgc                                                                                                                                                                                                                                 >  2:341099/1‑163 (MQ=255)
                                                                                                                               tGACCAGAATGTTTAGCTACACTTAAGTTAATTATGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGCATATGAGGCAGGTATTTCGCTCGCTTCt                                                                                                                                                                                                                             <  1:341111/163‑1 (MQ=255)
                                                                                                                                             aTCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCNTTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGaa                                                                                               <  1:378163/275‑1 (MQ=255)
                                                                                                                                                                 tGACAGGCGCATGTGGCAATAGCATGCCACTATTGAGTAAATCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCTACAACGGTACTGGCGGCTGATGTGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGCCCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGc                                                                           <  1:316744/275‑1 (MQ=255)
                                                                                                                                                                                          gCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGACCGGGCTTAAACTGCTGGAGCATATACGTAAcctggcgcttatac                                                   >  1:267315/1‑261 (MQ=255)
                                                                                                                                                                                             aCTATTGAGTAAAGCCAGTCAGTGGAGAGAACATGACGCAGCCATTGGTAGGAAAAGAGAATCTCATTGTTGAAGATGAGGCAGGTATGTCGCTCGGTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGGTGGGAGGTGTCACTCAAGACCTGATGATATGTGATATCGGGATGCCACGAATGAACGGGCTTAAActgc                                                                                 <  2:482744/241‑1 (MQ=255)
                                                                                                                                                                                              cTATTGAGTAAAGCCAGTCAGTGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATCTCGCGATGcc                                                                                                        <  1:479053/217‑1 (MQ=255)
                                                                                                                                                                                               tATTGAGTAAAGCCAGTCCGGGTAGAGAACATGACGCAGCCATTGGTCGGAAAGCAGATTCTCATTGTTGAAGATGAGGCAGGTAttt                                                                                                                                                                                                                                        >  1:311563/1‑88 (MQ=255)
                                                                                                                                                                                                   gAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTtcgct                                                                                                                                                                                                                                    <  2:311554/88‑1 (MQ=255)
                                                                                                                                                                                                           ccAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGatat                                 >  1:288683/1‑275 (MQ=255)
                                                                                                                                                                                                              gTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTAttt                                                                                                                                                                                                                                        >  1:263998/1‑73 (MQ=255)
                                                                                                                                                                                                                  ggggAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTtcgct                                                                                                                                                                                                                                    <  2:263990/73‑1 (MQ=255)
                                                                                                                                                                                                                               gACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGCGGTTTCAc                                                                                                                                         >  1:306222/1‑151 (MQ=255)
                                                                                                                                                                                                                               gACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTACTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGACTGAACGGGCTTACACTGCTGTAGCATATAGGTAACAGCGGCGAACCAAACCAAGTTCTGGTGATAGATGCCACTGAAAATATgg               >  2:230372/1‑273 (MQ=255)
                                                                                                                                                                                                                                   cAGCAATTGGTCGGACAACAGAATCTCATTGTTGAAGATGGGGCCGGTATTTCGCGCGCTTCTGGATTCATGGTTTTCCTCATTCGGAGCGACAACGGTACTGGCGGATGATGGGGTGGATGTCCTTGAGTTGCTGGGAGGTTTCACTCCa                                                                                                                                     <  2:306213/151‑1 (MQ=255)
                                                                                                                                                                                                                                    aGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATgggg                                                                                                                                                                        >  1:99482/1‑115 (MQ=255)
                                                                                                                                                                                                                                        aTTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGa                                                                                                                                                                    <  2:99480/115‑1 (MQ=255)
                                                                                                                                                                                                                                          ttgttggAAAACTGATTCTCATTTTTGAAGATGAGGTAGGTATTTCGCTCGCTTCTGGAGTCATGGGTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCaa  <  1:230379/270‑1 (MQ=255)
                                                                                                                                                                                                                                             tCGGAAAACAGATTCTCATTGTTGAAGATGAGGCAGGTATTTCNCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTTAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGACCATATGCCTAACAGAGGCGACCAGAcccc                                               <  1:352162/227‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
GATTCAACCTGTTTGCCCGCAAATATCTACGCTTGATTTCCATCGCGCGCACGCTGCCATTGCGGCCGGACAGCTGGCAGTGGAAAGGAAAATGGACGAACTTTTGCCGTTGGTACGCACCAACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAA  >  CP009273/1285477‑1285984

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: