Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,484,463 T→C 17.6% I38V (ATC→GTC)  azoR ← NADH‑azoreductase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,484,4630TC17.6% 155.9 / 28.5 68I38V (ATC→GTC) azoRNADH‑azoreductase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base C (7/5);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CGCAGAGCGCCAACCAGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACTGAGAG  >  W3110S.gb/1484396‑1484530
                                                                   |                                                                   
cGCAGAGCGCCAACCAGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGACTtc                                                                  <  1:497391/71‑1 (MQ=255)
   agagCGCCAACCAGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTtcatc                                                               <  1:1571332/71‑1 (MQ=255)
      gcgcCAACCAGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGc                                                            >  1:2130016/1‑71 (MQ=255)
       cgcCAACCAGTTCGCCATCCAGTATCGGAATCGGATTTGCAGCCAGGGCGCGAACGGTGATTtca                                                                 <  1:1209987/65‑1 (MQ=255)
           aaCCAGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGt                                                       >  1:3813012/1‑71 (MQ=255)
            aCCAGTACGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTg                                                      >  1:659276/1‑71 (MQ=255)
               aGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTCCATCAGCGGAg                                                        >  1:3670702/1‑66 (MQ=255)
                   cGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGTAGTGCTTTTcg                                               >  1:1206264/1‑71 (MQ=255)
                      cATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGATTGCTTTTCgcgc                                            <  1:3076178/71‑1 (MQ=255)
                      cATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAGTGCTTTTCgcgc                                            <  1:4537910/71‑1 (MQ=255)
                       aTCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCgcgc                                            >  1:2531844/1‑70 (MQ=255)
                       aTCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCgcgc                                            >  1:4666710/1‑70 (MQ=255)
                       aTCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGcc                                           <  1:4553904/71‑1 (MQ=255)
                        tCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTtc                                                                  <  1:1298128/47‑1 (MQ=255)
                        tCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTtc                                                                  <  1:5050564/47‑1 (MQ=255)
                        tCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTtc                                                                  <  1:1081587/47‑1 (MQ=255)
                        tCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTtc                                                                  <  1:4500614/47‑1 (MQ=255)
                        tCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCg                                                           <  1:4152135/54‑1 (MQ=255)
                          cAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCAt                                         <  1:3623675/70‑1 (MQ=255)
                          cAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCAt                                         <  1:4920643/70‑1 (MQ=255)
                           aGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCAt                                         <  1:3639173/69‑1 (MQ=255)
                            gTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTg                                       >  1:1814586/1‑70 (MQ=255)
                               ccGGAATCGGATTTGCATCCAGGTCGCGAACGGTGATTTCATCAGCTGAGTGCTTTTCGCGCCATTGTTca                                   <  1:2080547/71‑1 (MQ=255)
                                 ggAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAg                                                        >  1:937290/1‑48 (MQ=255)
                                   aaTCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTcaacaa                               <  1:2374575/71‑1 (MQ=255)
                                    ttCGGATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAGTGCTTTTCGCGCCATTGt                                      <  1:3524716/62‑1 (MQ=255)
                                    aTCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTcaacaa                               >  1:1669568/1‑70 (MQ=255)
                                    aTCGGATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAGTGCTTTTCGCGCCATTGt                                      <  1:5034998/63‑1 (MQ=255)
                                    aTCGGATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAGTGCTTTTCGCGCCATTGt                                      <  1:4233563/63‑1 (MQ=255)
                                     tCGGATTTGCAGCCAGGTCGCGAACGGTGATTtca                                                                 <  1:4430830/35‑1 (MQ=255)
                                     tCGGATTTGCAGCCAGGTCGCGAACGGTGATTtca                                                                 <  1:3986063/35‑1 (MQ=255)
                                     tCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTc                                                <  1:4454317/52‑1 (MQ=255)
                                     tCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTc                                                <  1:4313505/52‑1 (MQ=255)
                                     tCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACaaa                              >  1:892011/1‑70 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGATTtcatc                                                               >  1:2250317/1‑35 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTcaacaa                               >  1:3339518/1‑67 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAaata                           <  1:4306187/71‑1 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAg                                                        >  1:4426774/1‑42 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAg                                                        >  1:514466/1‑42 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAg                                                        >  1:3335535/1‑42 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAg                                                        >  1:3951733/1‑42 (MQ=255)
                                       ggATTTGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAg                                                        >  1:2468196/1‑42 (MQ=255)
                                         aTTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAaataat                         >  1:5083749/1‑71 (MQ=255)
                                         aTTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAaataa                          <  1:2199057/70‑1 (MQ=255)
                                           ttGCAGCCAGGTCGCGAACGGTGACTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCg                       >  1:1200951/1‑71 (MQ=255)
                                              cAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGAc                    <  1:2809369/71‑1 (MQ=255)
                                              cAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGAc                    <  1:454979/71‑1 (MQ=255)
                                                  cAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACt                <  1:3835836/71‑1 (MQ=255)
                                                    ggTCGCGAACGGTGACTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGa              >  1:2157265/1‑71 (MQ=255)
                                                     gTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTcaacaa                               <  1:3270872/53‑1 (MQ=255)
                                                     gTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTc                                    <  1:4413076/48‑1 (MQ=255)
                                                     gTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTc                                    <  1:3645720/48‑1 (MQ=255)
                                                       cgcgAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACaaaa                             >  1:2174963/1‑53 (MQ=255)
                                                       cgcgAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGAt             >  1:1835376/1‑69 (MQ=255)
                                                        gcgAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGAt             >  1:4135979/1‑68 (MQ=255)
                                                        gcgAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGAt             >  1:4622991/1‑68 (MQ=255)
                                                        gcgAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTa           <  1:2504296/70‑1 (MQ=255)
                                                         cgAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGa         <  1:3567450/71‑1 (MQ=255)
                                                          gAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTcaacaa                               <  1:2456360/48‑1 (MQ=255)
                                                             cGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAaata                           <  1:5138676/49‑1 (MQ=255)
                                                             cGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACTga     <  1:4066656/71‑1 (MQ=255)
                                                              ggTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACt       >  1:1922376/1‑68 (MQ=255)
                                                               gTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACaa                  >  1:765457/1‑56 (MQ=255)
                                                               gTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGAt             >  1:4907910/1‑61 (MQ=255)
                                                                tGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTa           <  1:4657932/62‑1 (MQ=255)
                                                                 gATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACTGagag  >  1:3467090/1‑70 (MQ=255)
                                                                 gATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACTGagag  >  1:3461508/1‑70 (MQ=255)
                                                                 gATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACTGagag  >  1:3346875/1‑70 (MQ=255)
                                                                   |                                                                   
CGCAGAGCGCCAACCAGTTCGCCATCCAGTACCGGAATCGGATTTGCAGCCAGGTCGCGAACGGTGATTTCATCAGCGGAGTGCTTTTCGCGCCATTGTTCAACAAAATAATCGGACAACTGATTAGACTGAGAG  >  W3110S.gb/1484396‑1484530

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: