Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,497,165 G→A 10.8% A65T (GCG→ACG)  fsaB → fructose‑6‑phosphate aldolase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,497,1650GA10.8% 136.5 / 5.4 46A65T (GCG→ACG) fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base G (23/18);  new base A (3/2);  total (26/20)
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

AGTGCTGCCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTGTGGTGAAAATCCCGGT  >  W3110S.gb/3497095‑3497235
                                                                      |                                                                      
aGTGCTGCCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcg                                                                        >  1:1089710/1‑71 (MQ=255)
aGTGCTGCCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcg                                                                        >  1:3381667/1‑71 (MQ=255)
   gctgcCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgc                                                                       >  1:3801753/1‑69 (MQ=255)
    ctgcCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgcg                                                                      >  1:3236935/1‑69 (MQ=255)
    ctgcCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACgcgc                                                                     <  1:880556/70‑1 (MQ=255)
       ccGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAggg                                                                 >  1:1426321/1‑71 (MQ=255)
       ccGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACACGCAgg                                                                  <  1:2321030/70‑1 (MQ=255)
           gTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgcg                                                                      >  1:318640/1‑62 (MQ=255)
              tGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGAtggtg                                                          >  1:1895175/1‑71 (MQ=255)
                 aaaaaGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCATGGGAtggtgg                                                         <  1:70825/69‑1 (MQ=255)
                 aaaaaGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGAtggtgg                                                         <  1:3810924/69‑1 (MQ=255)
                 aaaaaGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGgaa                                                       <  1:2317414/71‑1 (MQ=255)
                   aaaGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAgg                                                                  <  1:1147766/58‑1 (MQ=255)
                   aaaGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGAtg                                                             <  1:1652981/63‑1 (MQ=255)
                   aaaGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGgaaga                                                     <  1:3269085/71‑1 (MQ=255)
                     aGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGgaagaa                                                    >  1:1831150/1‑70 (MQ=255)
                      gCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgcg                                                                      <  1:1841436/51‑1 (MQ=255)
                      gCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagc                                                  <  1:3144252/71‑1 (MQ=255)
                         aTTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgc                                                                       <  1:3110117/47‑1 (MQ=255)
                         aTTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgc                                                                       <  1:1012020/47‑1 (MQ=255)
                         aTTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGAcgc                                                                       <  1:3606348/47‑1 (MQ=255)
                          ttGGTTATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagc                                                  <  1:122841/67‑1 (MQ=255)
                          ttGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagcgaag                                              >  1:4075628/1‑71 (MQ=255)
                          ttGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagcgaag                                              >  1:2821485/1‑71 (MQ=255)
                             gTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagcgaa                                               <  1:1590649/67‑1 (MQ=255)
                               gatgaGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCt                                         >  1:2089674/1‑71 (MQ=255)
                                atgaGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTg                                        >  1:755376/1‑71 (MQ=255)
                                atgaGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTg                                        >  1:571212/1‑71 (MQ=255)
                                 tgaGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTgc                                       >  1:1402213/1‑71 (MQ=255)
                                   tgggCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGTGAtggtgg                                                         <  1:1431051/50‑1 (MQ=255)
                                   aGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGAtggtg                                                          >  1:425707/1‑50 (MQ=255)
                                   aGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGgaag                                                      >  1:429334/1‑54 (MQ=255)
                                   aGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGCAGAAGCGAAGCGCCTgcgc                                     <  1:3502646/71‑1 (MQ=255)
                                    gggCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGcgcg                                    <  1:2463085/71‑1 (MQ=255)
                                        aTTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGc                                <  1:2327922/71‑1 (MQ=255)
                                         ttCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGc                                >  1:2922103/1‑70 (MQ=255)
                                           cTGTTTGCTCAGACCATGAGCCGCGACACGCAGGGGATGGTGGAAGAAGCGAAgcgc                                           >  1:2576805/1‑57 (MQ=255)
                                            tGTTTTCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTAtt                            >  1:3138620/1‑71 (MQ=255)
                                            tGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTAtt                            >  1:1186928/1‑71 (MQ=255)
                                            tGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTAt                             >  1:3153153/1‑70 (MQ=255)
                                                  cTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGga                                                        <  1:2392634/37‑1 (MQ=255)
                                                    cAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTgcgc                                     <  1:3168681/54‑1 (MQ=255)
                                                     aGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTgcgc                                     >  1:3164945/1‑53 (MQ=255)
                                                     aGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTAtt                   >  1:2995690/1‑71 (MQ=255)
                                                     aGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTAtt                   <  1:3943004/71‑1 (MQ=255)
                                                      gACCATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagcgaag                                              >  1:4172383/1‑43 (MQ=255)
                                                        ccATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATtgtg                >  1:3455454/1‑71 (MQ=255)
                                                        ccATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATtgt                 <  1:943613/70‑1 (MQ=255)
                                                         cATGAGCCGCGACGCGCAGGGGATGGTGGAAgaagcgaagccc                                           >  1:4033642/1‑41 (MQ=255)
                                                          aTGAGCCGCGACGCGCAGGGGATGGTGGAAgaagcg                                                 <  1:1395996/36‑1 (MQ=255)
                                                           tGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTgtggt              >  1:20430/1‑70 (MQ=255)
                                                           tGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTgtggt              >  1:1141318/1‑70 (MQ=255)
                                                           tGAGCCGCGACACGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGt                      >  1:2866291/1‑62 (MQ=255)
                                                           tGAGCCGCGACACGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTgtggt              >  1:1031164/1‑70 (MQ=255)
                                                                cgcgACACGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATtgtg                <  1:1436635/63‑1 (MQ=255)
                                                                 gcgACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTGTGGTGa            >  1:425226/1‑66 (MQ=255)
                                                                      gcgcAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTGTGGTGAAAATCCCGGt  <  1:2104332/71‑1 (MQ=255)
                                                                      |                                                                      
AGTGCTGCCGCGTCTGCAAAAAGCGATTGGTGATGAGGGCATTCTGTTTGCTCAGACCATGAGCCGCGACGCGCAGGGGATGGTGGAAGAAGCGAAGCGCCTGCGCGACGCTATTCCGGGTATTGTGGTGAAAATCCCGGT  >  W3110S.gb/3497095‑3497235

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: