Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,020,867 A→G 57.6% S159P (TCC→CCC)  sulA ← SOS cell division inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,020,8670AG57.6% 15.2 / 37.8 33S159P (TCC→CCC) sulASOS cell division inhibitor
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base G (6/13);  total (11/22)
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.56e-01

AGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCA  >  W3110S.gb/1020806‑1020934
                                                             |                                                                   
aGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAg                                                                  >  1:3640513/1‑65 (MQ=255)
aGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAg                                                                  >  1:504044/1‑65 (MQ=255)
       ggATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGtctc                                                           <  1:1108750/65‑1 (MQ=255)
               ccTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGgaa                                                                   <  1:3340549/49‑1 (MQ=255)
               ccTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGagag                                               <  1:2664757/69‑1 (MQ=255)
                cTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAgg                                              >  1:1815488/1‑69 (MQ=255)
                  aaaTTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGAt                                            >  1:3022234/1‑69 (MQ=255)
                  aaaTTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGAt                                            >  1:590380/1‑69 (MQ=255)
                     tttACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATgc                                          <  1:3192736/68‑1 (MQ=255)
                      ttACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATgc                                          <  1:876881/67‑1 (MQ=255)
                       tACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCt                                       <  1:2565356/69‑1 (MQ=255)
                          ttAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTa                                     <  1:631668/68‑1 (MQ=255)
                          ttAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTa                                     <  1:1756810/68‑1 (MQ=255)
                                  aCAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTgag                                                 <  1:3408742/48‑1 (MQ=255)
                                   cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATgcg                                         >  1:2394840/1‑55 (MQ=255)
                                   cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATgcg                                         >  1:2370913/1‑55 (MQ=255)
                                   cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATgcg                                         >  1:2137628/1‑55 (MQ=255)
                                   cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATgcg                                         <  1:681825/55‑1 (MQ=255)
                                   cAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATgcg                                         >  1:3337261/1‑55 (MQ=255)
                                       ttAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       <  1:2215232/69‑1 (MQ=255)
                                       ttAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       >  1:3469694/1‑69 (MQ=255)
                                       ttAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       <  1:3092586/69‑1 (MQ=255)
                                       ttAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       <  1:1703324/69‑1 (MQ=255)
                                            gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAccc                  <  1:2904250/69‑1 (MQ=255)
                                            gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACAGGACGCATAATAAAccc                  <  1:2448250/69‑1 (MQ=255)
                                               aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       <  1:2547224/61‑1 (MQ=255)
                                               aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       <  1:2327838/61‑1 (MQ=255)
                                               aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                       >  1:2056131/1‑61 (MQ=255)
                                                 tttttAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg             <  1:2685977/69‑1 (MQ=255)
                                                  ttttAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg             <  1:429739/68‑1 (MQ=255)
                                                   tttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCa                            <  1:816015/52‑1 (MQ=255)
                                                   tttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCa                            <  1:1185355/52‑1 (MQ=255)
                                                   attAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCa                            <  1:2561247/51‑1 (MQ=255)
                                                            gAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCa  <  1:3630328/69‑1 (MQ=255)
                                                             |                                                                   
AGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTTACCTTCA  >  W3110S.gb/1020806‑1020934

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