Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 703,662 T→C 8.9% K124K (AAA→AAG nagB ← glucosamine‑6‑phosphate deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb703,6620TC 8.9% 118.3 / 4.1 45K124K (AAA→AAGnagBglucosamine‑6‑phosphate deaminase
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/2);  total (23/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.33e-01

AAATGCAATATGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGTTTTCT  >  W3110S.gb/703596‑703731
                                                                  |                                                                     
aaaTGCAATATGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAttttt                                                                    <  1:2572259/70‑1 (MQ=255)
        tatGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttt                                                                     <  1:4243091/61‑1 (MQ=255)
        tatGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttt                                                                     <  1:4449367/61‑1 (MQ=255)
        tatGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATAc                                                           <  1:1005232/71‑1 (MQ=255)
         atGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATAc                                                           >  1:50870/1‑70 (MQ=255)
          tGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTg                                                         >  1:1106349/1‑71 (MQ=255)
          tGACCGTCGTTACCTACACCGCCCATAAACAGATCAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTg                                                         <  1:4124940/71‑1 (MQ=255)
           gACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATCTTttc                                                                  >  1:2576694/1‑61 (MQ=255)
                tcgtTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa                                                              <  1:241604/60‑1 (MQ=255)
                 cgtTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa                                                              >  1:1063760/1‑59 (MQ=255)
                  gtTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACt                                                 >  1:2313007/1‑71 (MQ=255)
                    tACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCTGCACTc                                                <  1:9345/70‑1 (MQ=255)
                     aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttt                                                                     >  1:2410687/1‑48 (MQ=255)
                     aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa                                                            >  1:2109303/1‑57 (MQ=255)
                     aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCgg                                              >  1:1516886/1‑71 (MQ=255)
                     aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCgg                                              >  1:3597836/1‑71 (MQ=255)
                     aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATCTTTTCTTCATACTGGCGGCACTCgg                                              >  1:3552935/1‑71 (MQ=255)
                       cTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACt                                                          <  1:1530668/57‑1 (MQ=255)
                        tACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa                                                              <  1:2302147/52‑1 (MQ=255)
                         acacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttttt                                                                   >  1:3083516/1‑46 (MQ=255)
                         acacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTttct                                                                 >  1:3450918/1‑48 (MQ=255)
                         acacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa                                                            >  1:3961957/1‑53 (MQ=255)
                          cacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa                                                              >  1:1055703/1‑50 (MQ=255)
                          cacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa                                                              >  1:1423521/1‑50 (MQ=255)
                           acCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGa                                        <  1:789977/71‑1 (MQ=255)
                              gCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATAc                                                           >  1:864063/1‑49 (MQ=255)
                                ccATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc                                    <  1:2405670/70‑1 (MQ=255)
                                ccATAAACAGATGAATTTTTCCGTAAGAACGGATCTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc                                    <  1:2629184/70‑1 (MQ=255)
                                   tAAACAGATGAATTTTTCCGTAAGAACGGAtttttt                                                                   <  1:4538330/36‑1 (MQ=255)
                                   tAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc                                    <  1:45351/67‑1 (MQ=255)
                                   tAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCggg                                >  1:2836006/1‑71 (MQ=255)
                                   tAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCgg                                 <  1:902098/70‑1 (MQ=255)
                                    aaaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa                                                            >  1:867047/1‑42 (MQ=255)
                                    aaaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCgggg                               >  1:2767919/1‑71 (MQ=255)
                                     aaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGc                              >  1:4114768/1‑71 (MQ=255)
                                     aaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGAGGCACTCGGCGt                                           >  1:1314094/1‑58 (MQ=255)
                                      aCAGATGAATTTTTCCGTAAGAACGGATTTTttct                                                                 <  1:3334032/35‑1 (MQ=255)
                                        aGATGAATTTTTCCGTAAGAACGGATTTTttcttc                                                               >  1:4501037/1‑35 (MQ=255)
                                        aGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGtt                           <  1:19808/71‑1 (MQ=255)
                                         gATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa                                                            <  1:1576736/37‑1 (MQ=255)
                                          aTGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGa                                        <  1:1747911/56‑1 (MQ=255)
                                                 tttCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGata                                      <  1:1328191/51‑1 (MQ=255)
                                                   tCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGa                                        >  1:1310002/1‑47 (MQ=255)
                                                   tCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc                                    <  1:183590/51‑1 (MQ=255)
                                                        aaGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGtt            >  1:521637/1‑70 (MQ=255)
                                                          gAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGtt            <  1:4528665/68‑1 (MQ=255)
                                                             cGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGtt      <  1:739123/71‑1 (MQ=255)
                                                             cGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGtt      <  1:3779968/71‑1 (MQ=255)
                                                             cGGATCTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGtt      <  1:2152143/71‑1 (MQ=255)
                                                              ggATTTTTTCTTCATACTGGCGGCACTCGGCGTCGata                                      <  1:2533647/38‑1 (MQ=255)
                                                                aTTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGt                            >  1:340320/1‑46 (MQ=255)
                                                                 ttttttCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGTTTTct  >  1:4339549/1‑71 (MQ=255)
                                                                  |                                                                     
AAATGCAATATGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGTTTTCT  >  W3110S.gb/703596‑703731

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