Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,620,666 T→G 15.9% L492W (TTG→TGG)  rhsB → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,620,6660TG15.9% 60.1 / 3.5 25L492W (TTG→TGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base G (2/2);  total (12/13)
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.85e-01

TAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACC  >  NC_000913/3620612‑3620908
                                                      |                                                                                                                                                                                                                                                  
gcccgtcaggcccggtggCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGcccc                                                                                                                                                                                                                         >  5:680945‑M2/19‑82 (MQ=255)
gcccgtcaggcccggtggCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGcccc                                                                                                                                                                                                                         <  6:680945‑M2/64‑1 (MQ=255)
 cccgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGCCGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg                                                                                                                                                                             >  5:307187‑M2/18‑125 (MQ=255)
 cccgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg                                                                                                                                                                             <  6:307187‑M2/108‑1 (MQ=255)
   cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTgg                                                                                                                                                                                                                                                  >  8:1007620‑M2/16‑54 (MQ=255)
   cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTgg                                                                                                                                                                                                                                                  <  7:1007620‑M2/39‑1 (MQ=255)
   cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTa                                                                                                                                                                                              >  8:808779‑M2/16‑106 (MQ=255)
   cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTa                                                                                                                                                                                              <  7:808779‑M2/91‑1 (MQ=255)
   cgtcaggcccggtggCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg                                                                                                                                                                             >  8:276045‑M2/16‑123 (MQ=255)
   cgtcaggcccgggggCTGACGGGCGGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg                                                                                                                                                                             <  7:276045‑M2/108‑1 (MQ=255)
                   TGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCGGGAAACTGCCCCTGATGGCGa                                                                                                                                                                                                                >  8:173027‑M2/1‑72 (MQ=255)
                   TGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCGGGAAACTGCCCCTGACGGCGa                                                                                                                                                                                                                <  7:173027‑M2/72‑1 (MQ=255)
                   TGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCCCACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCAGGAAAACCATGCCGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCAACTGTTCCGGTTATGTCACCAGTTATGACCCTGCCCGCTTCTGGCAGCTGACGACTGTCCACCGCgcggaagagcaggacaagt     >  1:122828‑M2/1‑266 (MQ=255)
                    GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGtt                                                                                                                                                                                               <  7:806904‑M2/88‑1 (MQ=255)
                    GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGtt                                                                                                                                                                                               >  8:806904‑M2/1‑88 (MQ=255)
                    GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATaa                                                                                                                                                                                        >  7:262918‑M2/1‑95 (MQ=255)
                    GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATaa                                                                                                                                                                                        <  8:262918‑M2/95‑1 (MQ=255)
                    GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGaaaa                                                                                                                                   >  6:243401‑M2/1‑148 (MQ=255)
                     ACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAAc                                                                                                                                  >  8:270654‑M2/1‑148 (MQ=255)
                      CGGGCTGGAATTGCGCCGGGAATATGATGAATTGGTCCGTCTGATGCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAAcc                                                                                                                                 <  8:236934‑M2/148‑1 (MQ=255)
                      CGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAAcc                                                                                                                                 >  7:552730‑M2/1‑148 (MQ=255)
                       GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGTAACGGAAGATGCCACCGGCAGCCGGAAAACCa                                                                                                                                <  6:227526‑M2/148‑1 (MQ=255)
                       GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCa                                                                                                                                <  7:803460‑M2/148‑1 (MQ=255)
                       GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCa                                                                                                                                <  5:243401‑M2/148‑1 (MQ=255)
                       GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTAcc  <  2:277833‑M2/274‑1 (MQ=255)
                                                      |                                                                                                                                                                                                                                                  
TAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACC  >  NC_000913/3620612‑3620908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: