Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,096,084 G→A 87.5% E160E (GAG→GAA hisB → fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,096,0840GA87.5% 96.1 / 7.6 40E160E (GAG→GAAhisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (15/20);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

AAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGA  >  W3110S.gb/2096020‑2096146
                                                                |                                                              
aaaaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACaa                                                             <  1:2400736/68‑1 (MQ=255)
aaaaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACaa                                                             <  1:2267174/68‑1 (MQ=255)
  aaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCAc                                                        >  1:2435853/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCg                                            >  1:371920/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCg                                            >  1:1943881/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                            >  1:922849/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                            >  1:1300120/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                            >  1:137487/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                            >  1:455896/1‑71 (MQ=255)
              aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCg                                            >  1:2242865/1‑71 (MQ=255)
               cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGt                                           <  1:1457633/71‑1 (MQ=255)
                  gTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTAc                                        >  1:2449023/1‑71 (MQ=255)
                       cGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCa                                   <  1:2691086/71‑1 (MQ=255)
                       cGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCa                                   >  1:746349/1‑71 (MQ=255)
                          tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACt                                                           <  1:1030340/44‑1 (MQ=255)
                            cGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGt                                <  1:2410929/69‑1 (MQ=255)
                            cGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAg                              <  1:612720/71‑1 (MQ=255)
                             gACCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAg                              <  1:191845/70‑1 (MQ=255)
                              aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGAcc                                             <  1:2102479/54‑1 (MQ=255)
                              aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGAcc                                             <  1:2504839/54‑1 (MQ=255)
                               ccGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTa                                         >  1:1764112/1‑57 (MQ=255)
                                 gcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                         <  1:441692/71‑1 (MQ=255)
                                 gcgAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                         <  1:1875806/71‑1 (MQ=255)
                                 gcgAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                         <  1:2315646/71‑1 (MQ=255)
                                 gcgAAACCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                         <  1:321752/71‑1 (MQ=255)
                                     aaCCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGt                                           <  1:1232486/49‑1 (MQ=255)
                                     aaCCCTGAACTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATa                     <  1:1664182/71‑1 (MQ=255)
                                       cccTGAACTGGCCAATGATTGGCGAACAACTCACCAGACg                                                  >  1:986932/1‑40 (MQ=255)
                                       cccTGAACTGGCCAATGATTGGCGAACAACTCAACAGACg                                                  >  1:1354202/1‑40 (MQ=255)
                                            aaCTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga              <  1:1811105/71‑1 (MQ=255)
                                            aaCTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAaga              <  1:124619/71‑1 (MQ=255)
                                             aCTGGCCAATGATTGGCGAACAACTCACCAGACGTGACCGATACGCTCACGTAGTGCGTaa                       >  1:1420965/1‑61 (MQ=255)
                                                  gcAATGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACg                                 <  1:411760/45‑1 (MQ=255)
                                                  ccAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCa        <  1:469867/71‑1 (MQ=255)
                                                     aTGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCa                  <  1:1398210/58‑1 (MQ=255)
                                                      tGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                       >  1:421952/1‑52 (MQ=255)
                                                      tGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                       >  1:57856/1‑52 (MQ=255)
                                                      tGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt    <  1:672653/71‑1 (MQ=255)
                                                      tGATTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGAtt    >  1:897763/1‑71 (MQ=255)
                                                        aTTGGCGAACAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGa  >  1:2424521/1‑71 (MQ=255)
                                                                |                                                              
AAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGA  >  W3110S.gb/2096020‑2096146

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: