Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,020,867 A→G 68.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,8670AG68.6% 42.4 / 25.4 35S159P (TCC→CCC) sulASOS cell division inhibitor
Reads supporting (aligned to +/- strand):  ref base A (5/6);  new base G (22/2);  total (27/8)
Fisher's exact test for biased strand distribution p-value = 5.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTT  >  W3110S.gb/1020813‑1020927
                                                      |                                                            
ggATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGc                                               >  1:2984401/1‑70 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:1270929/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:3222911/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:2968204/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:2757073/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:735339/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:2135137/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:1822152/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:928097/1‑56 (MQ=255)
     aaTCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTg                                                        >  1:788056/1‑56 (MQ=255)
                      aTGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGa                         <  1:778636/70‑1 (MQ=255)
                      aTGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGa                         <  1:2204706/70‑1 (MQ=255)
                                    agTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAcacc          <  1:1863626/71‑1 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAAcccc          >  1:1824933/1‑70 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:2359054/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:3004882/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:1647052/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:646407/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:739519/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:2148503/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:2061252/1‑71 (MQ=255)
                                     gTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa         >  1:1900445/1‑71 (MQ=255)
                                        aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg      >  1:772814/1‑71 (MQ=255)
                                        aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg      >  1:1054178/1‑71 (MQ=255)
                                        aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg      >  1:930950/1‑71 (MQ=255)
                                        aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg      >  1:2871654/1‑71 (MQ=255)
                                           ttttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATaaa              >  1:3149969/1‑60 (MQ=255)
                                           ttttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGt   >  1:312655/1‑71 (MQ=255)
                                           ttttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGt   >  1:1784598/1‑71 (MQ=255)
                                           ttttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGt   >  1:1695531/1‑71 (MQ=255)
                                             ttAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGtt  <  1:2097026/70‑1 (MQ=255)
                                             ttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGtt  <  1:2030949/70‑1 (MQ=255)
                                             ttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGtt  <  1:1389853/70‑1 (MQ=255)
                                             ttAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCGTACCGGACGCATAATAAACCCCATAGCGtt  <  1:2658110/70‑1 (MQ=255)
                                                 cccGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata                <  1:2636658/52‑1 (MQ=255)
                                                      |                                                            
GGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAGCGTT  >  W3110S.gb/1020813‑1020927

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