Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,682,573 G→A L211L (CTG→TTG)  ilvD ← dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,682,5730GA80.0% 67.2 / 14.7 35L211L (CTG→TTG) ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/7);  new base A (9/19);  total (9/26)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

ACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGACA  >  W3110S.gb/3682520‑3682639
                                                     |                                                                  
aCGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAg                                                    <  1:778757/70‑1 (MQ=255)
aCGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAg                                                    <  1:3929602/70‑1 (MQ=255)
aCGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAg                                                    <  1:3440788/70‑1 (MQ=255)
aCGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAg                                                    <  1:2935939/70‑1 (MQ=255)
  ggTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTc                                                 <  1:1610422/71‑1 (MQ=255)
  ggTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGGCAGGCAGTTc                                                 <  1:3715311/71‑1 (MQ=255)
       gCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGa                                            <  1:3309132/71‑1 (MQ=255)
            ggTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAg                                       >  1:719868/1‑71 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:1104287/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:768388/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:4065545/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:3612676/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:1191269/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:2359584/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:1511882/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAg                                           >  1:2148487/1‑60 (MQ=255)
                   agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:1084102/71‑1 (MQ=255)
                       gcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATc                            <  1:2865737/71‑1 (MQ=255)
                       gcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATc                            <  1:3076347/71‑1 (MQ=255)
                       gcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATc                            <  1:1419059/71‑1 (MQ=255)
                       gcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATc                            <  1:2161121/71‑1 (MQ=255)
                       gcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATc                            <  1:1655147/71‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:244949/48‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:805525/48‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:362914/48‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:3516703/48‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:3360935/48‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:2666817/48‑1 (MQ=255)
                                          cGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTaaa                                <  1:2398848/48‑1 (MQ=255)
                                               ggCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGa    <  1:3301040/71‑1 (MQ=255)
                                                 cccAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGaca  <  1:1703341/71‑1 (MQ=255)
                                                 cccAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGaca  <  1:1348639/71‑1 (MQ=255)
                                                 cccAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGaca  <  1:2929345/71‑1 (MQ=255)
                                                 cccAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGaca  <  1:904796/71‑1 (MQ=255)
                                                 cccAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGGCGGaca  <  1:9223/71‑1 (MQ=255)
                                                     |                                                                  
ACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCGGACA  >  W3110S.gb/3682520‑3682639

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: