Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,348,085 G→A 25.6% E89K (GAA→AAA)  npr → phosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,348,0850GA25.6% 78.0 / 24.5 43E89K (GAA→AAA) nprphosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr)
Reads supporting (aligned to +/- strand):  ref base G (19/13);  new base A (6/5);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCA  >  W3110S.gb/3348019‑3348150
                                                                  |                                                                 
gCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTtgatg                                                                   >  1:1186854/1‑67 (MQ=255)
   aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGAtt                                                             >  1:1798478/1‑70 (MQ=255)
   aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTa                                                            <  1:1853257/71‑1 (MQ=255)
    ccGGGCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTaa                                                           <  1:436406/71‑1 (MQ=255)
        tCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCtt                                                       <  1:2234275/71‑1 (MQ=255)
           acaGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTtcat                                                    <  1:2761567/71‑1 (MQ=255)
                 agaggaagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGaa                                                                 <  1:1678953/52‑1 (MQ=255)
                       agCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATa                                        <  1:1822254/71‑1 (MQ=255)
                        gCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGAtt                                                             <  1:1783989/49‑1 (MQ=255)
                            tGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc                                       <  1:2175848/67‑1 (MQ=255)
                            tGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTa                                                            <  1:1386574/46‑1 (MQ=255)
                             ggCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTtc                                                      >  1:455986/1‑51 (MQ=255)
                                cgccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGAc                                           >  1:606357/1‑59 (MQ=255)
                                    gTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCa                            >  1:2279338/1‑70 (MQ=255)
                                     ttATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTa                                                            <  1:457721/37‑1 (MQ=255)
                                     ttATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATAc                                       <  1:66175/58‑1 (MQ=255)
                                     ttATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATAc                                       <  1:2237751/58‑1 (MQ=255)
                                     ttATCGCACTCTTTAATTCTGGTTTTGATAAAGATTa                                                            <  1:1292200/37‑1 (MQ=37)
                                       aTCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACt                        >  1:274040/1‑71 (MQ=255)
                                       aTCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACaa                                     >  1:2662207/1‑58 (MQ=255)
                                              tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTc                       <  1:1622907/65‑1 (MQ=255)
                                               ctTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATa                                        >  1:2243322/1‑47 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCAc                                                 >  1:272126/1‑37 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCAc                                                 >  1:2129412/1‑37 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACttt                                              <  1:1232403/40‑1 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc                                       <  1:221756/47‑1 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTccc                     >  1:1644324/1‑65 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTccc                     >  1:2361700/1‑65 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGa                >  1:1227194/1‑70 (MQ=255)
                                                tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGa                >  1:2727912/1‑70 (MQ=255)
                                                  tAATTCTGGTTTTGATGAAGATTAATCTTCATCACtt                                               >  1:1484642/1‑37 (MQ=255)
                                                  tAATTCTGGTTTTGATGAAGATTAATCTTCATCACtt                                               >  1:2675082/1‑37 (MQ=255)
                                                  tAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa             >  1:217902/1‑71 (MQ=255)
                                                  tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGa                                            >  1:895923/1‑40 (MQ=255)
                                                  tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGa                                            >  1:1273801/1‑40 (MQ=255)
                                                  tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa             >  1:1019832/1‑71 (MQ=255)
                                                       cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTccc                                >  1:1560461/1‑47 (MQ=255)
                                                       cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTccc                                >  1:134946/1‑47 (MQ=255)
                                                       cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTccc                                >  1:1032393/1‑47 (MQ=255)
                                                       cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAg         >  1:1192376/1‑70 (MQ=255)
                                                       cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAg         <  1:113842/70‑1 (MQ=255)
                                                        tGGTTTTGATGAAGATTAATCTTCATCACTTTGACa                                          >  1:1913415/1‑36 (MQ=255)
                                                           ttttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaa       <  1:490878/68‑1 (MQ=255)
                                                            tttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAAc                         >  1:1488686/1‑49 (MQ=255)
                                                              tgatgaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCa  <  1:2188290/70‑1 (MQ=255)
                                                                  |                                                                 
GCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCA  >  W3110S.gb/3348019‑3348150

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: