Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,605,956 C→T 33.3% V82I (GTT→ATT)  dnaT ← DNA biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,605,9560CT33.3% 41.0 / 29.9 36V82I (GTT→ATT) dnaTDNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (7/5);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAG  >  W3110S.gb/4605888‑4606022
                                                                    |                                                                  
cAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACgg                                                                  <  1:313464/71‑1 (MQ=255)
 aGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACggg                                                                 <  1:2575079/71‑1 (MQ=255)
                   tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCt                                                >  1:1254369/1‑70 (MQ=255)
                   tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCt                                                >  1:1313087/1‑70 (MQ=255)
                     gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGa                                              >  1:113613/1‑70 (MQ=255)
                     gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGa                                              >  1:909557/1‑70 (MQ=255)
                     gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGa                                              >  1:1892644/1‑70 (MQ=255)
                       gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATag                                           <  1:3534128/71‑1 (MQ=255)
                           ccAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATaga                                          <  1:465188/68‑1 (MQ=255)
                           ccAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATaga                                          <  1:758827/68‑1 (MQ=255)
                               tctggATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGa                                              <  1:2509624/57‑1 (MQ=255)
                               tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCg                                                      <  1:1689556/52‑1 (MQ=255)
                               tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGa                                              <  1:2386709/60‑1 (MQ=255)
                                ccGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCg                                  >  1:989718/1‑71 (MQ=255)
                                 cGGATACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGt                                 >  1:642801/1‑71 (MQ=255)
                                 cGGATACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGt                                 >  1:2395640/1‑71 (MQ=255)
                                 cGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGt                                 >  1:2926866/1‑71 (MQ=255)
                                 cGGATACATGGCGAATTTCCCCATGGGTAAAGCAACGGGAg                                                               >  1:3150103/1‑41 (MQ=255)
                                     tACATGGCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTc                                >  1:917518/1‑68 (MQ=255)
                                        aTGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGc                          <  1:205114/71‑1 (MQ=255)
                                           gCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCg                                                      >  1:1794261/1‑40 (MQ=255)
                                           gCGAATTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGa                        >  1:525423/1‑70 (MQ=255)
                                           gCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATAGTCCAGAGCGa                        <  1:3645904/70‑1 (MQ=255)
                                              aaTTTCCCCATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATc                    <  1:1499776/71‑1 (MQ=255)
                                                tttCCCCATGGGTACAGCAATGGGAGCGGCTTGCTGTTCCTGATa                                            <  1:3899095/45‑1 (MQ=38)
                                                 ttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTCATCGTCCAGAGCGACATCGCt                 >  1:47225/1‑71 (MQ=255)
                                                    cccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATag                                           <  1:3445893/42‑1 (MQ=255)
                                                    cccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATag                                           <  1:1721727/42‑1 (MQ=255)
                                                    cccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTcc              <  1:2598812/71‑1 (MQ=255)
                                                    cccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTcc              <  1:1799673/71‑1 (MQ=255)
                                                      cATGGGTACAGCAATGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCg            <  1:1646677/71‑1 (MQ=255)
                                                        tGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGa          >  1:173444/1‑71 (MQ=255)
                                                            tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTg                                      <  1:386743/39‑1 (MQ=255)
                                                            tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTg                                      <  1:1921937/39‑1 (MQ=255)
                                                             aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGa                        >  1:2240400/1‑52 (MQ=255)
                                                                gCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAg  >  1:1092704/1‑71 (MQ=255)
                                                                    |                                                                  
CAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAG  >  W3110S.gb/4605888‑4606022

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