Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,371,273 Δ2 bp 100% intergenic (‑6/+299) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,371,2730A.100.0% 217.2 / NA 54intergenic (‑6/+300)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (25/29);  total (25/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*W3110S.gb4,371,2740A.98.1% 209.8 / ‑1.3 54intergenic (‑7/+299)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (24/29);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTA  >  W3110S.gb/4371203‑4371336
                                                                      ||                                                              
gggTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGa                                                                 <  1:3989105/71‑1 (MQ=255)
gggTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGa                                                                 <  1:3022682/71‑1 (MQ=255)
 ggTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa                                                                <  1:2298413/71‑1 (MQ=255)
   ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                              <  1:1294690/71‑2 (MQ=255)
   ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                              <  1:3163228/71‑2 (MQ=255)
   ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                               >  1:562001/1‑70 (MQ=255)
           cagcagCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                     >  1:1169477/1‑70 (MQ=255)
              cagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                  <  1:1693216/70‑1 (MQ=255)
                gcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa                                                                <  1:4111279/56‑1 (MQ=255)
                 cTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCg                                              <  1:2494004/71‑1 (MQ=255)
                     tctgtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCAt                                          >  1:304575/1‑71 (MQ=255)
                       tgtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc                                                <  1:2770342/63‑1 (MQ=255)
                         tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGAt                                       <  1:1218070/70‑1 (MQ=255)
                         tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGAt                                       <  1:3054353/70‑1 (MQ=255)
                            gatgaACCAGATAGTCCAGACCGCCAGCAACCTGCATTGCAG‑‑AAATAGCGGCGATAACCGCCATGATa                                      >  1:2811867/1‑68 (MQ=255)
                             atgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                     >  1:3671090/1‑68 (MQ=255)
                              tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg                                   >  1:3567417/1‑69 (MQ=255)
                              tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga                                 <  1:2711354/71‑1 (MQ=255)
                              tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga                                 <  1:2549795/71‑1 (MQ=255)
                                   cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa                            <  1:4013627/71‑1 (MQ=255)
                                    aGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                  >  1:2191592/1‑48 (MQ=255)
                                      aTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGTCGATAACCGCCATGATAATGgag                                <  1:812184/64‑1 (MQ=255)
                                      aTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATgg                                  <  1:3698772/62‑1 (MQ=255)
                                      aTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag                                <  1:1418908/64‑1 (MQ=255)
                                       tAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAcc                                               >  1:2896507/1‑48 (MQ=38)
                                        aGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag                                >  1:2105229/1‑62 (MQ=255)
                                         gTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCg                       <  1:2782855/70‑1 (MQ=255)
                                         gTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGa                      <  1:45308/71‑1 (MQ=255)
                                            cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAAc                    <  1:2474549/70‑1 (MQ=255)
                                            cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAAc                    >  1:232745/1‑70 (MQ=255)
                                            cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACg                   <  1:2075935/71‑1 (MQ=255)
                                            cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACg                   >  1:649568/1‑71 (MQ=255)
                                             aGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                  >  1:2010407/1‑71 (MQ=255)
                                             aGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACg                   >  1:2210902/1‑70 (MQ=255)
                                               aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATa                                      <  1:3755525/49‑1 (MQ=38)
                                               aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGa                >  1:1786504/1‑71 (MQ=255)
                                               aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGa                >  1:2181206/1‑71 (MQ=255)
                                                  gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt             >  1:3187318/1‑71 (MQ=255)
                                                  gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt             >  1:2336563/1‑71 (MQ=255)
                                                   ccAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGAt                                       <  1:2057913/44‑1 (MQ=38)
                                                   ccAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGAc                           >  1:714107/1‑56 (MQ=38)
                                                    cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTa           <  1:3957165/71‑1 (MQ=255)
                                                      gCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                     <  1:2247518/43‑1 (MQ=37)
                                                      gCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTAcc         >  1:3387390/1‑71 (MQ=255)
                                                        aaCCTGCATGGCAGAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCg        >  1:2608452/1‑71 (MQ=255)
                                                        aaCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCgg       >  1:3422000/1‑71 (MQ=255)
                                                         aCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                              >  1:2865139/1‑47 (MQ=37)
                                                         aCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                  <  1:1956823/59‑1 (MQ=39)
                                                           cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATgg                                  <  1:831419/41‑1 (MQ=25)
                                                           cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCg                       <  1:3997135/52‑1 (MQ=255)
                                                           cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGa                >  1:71690/1‑59 (MQ=39)
                                                           cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGtt    >  1:1331065/1‑71 (MQ=255)
                                                           cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt     <  1:3337696/70‑1 (MQ=255)
                                                            tGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGtt    >  1:3379346/1‑70 (MQ=255)
                                                             gCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGAt               >  1:4152883/1‑58 (MQ=39)
                                                              cATGGCAG‑‑AAATAGCGGCTATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt   <  1:4283820/69‑1 (MQ=255)
                                                              cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCggg      <  1:3585564/66‑1 (MQ=255)
                                                              cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt   <  1:3763950/69‑1 (MQ=255)
                                                              cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt   <  1:1114537/69‑1 (MQ=255)
                                                              cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt   <  1:1877237/69‑1 (MQ=255)
                                                              cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa  <  1:1520484/70‑1 (MQ=255)
                                                                      ||                                                              
GGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTA  >  W3110S.gb/4371203‑4371336

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: