Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 50,503 T→A 10.5% I240I (ATA→ATT apaH ← diadenosine tetraphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb50,5030TA10.5% 93.6 / 4.7 38I240I (ATA→ATTapaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base T (12/22);  new base A (2/2);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TTTATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGA  >  W3110S.gb/50440‑50573
                                                               |                                                                      
tttATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCt                                                                  <  1:820487/70‑1 (MQ=255)
  tATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGtata                                                                       >  1:1409253/1‑63 (MQ=255)
          cAGCGGAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCg                                                           <  1:3949269/67‑1 (MQ=255)
           aGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCt                                                      >  1:3326912/1‑71 (MQ=255)
             cgcAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCttt                                                    <  1:502629/71‑1 (MQ=255)
              gcaggcagGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTg                                                   <  1:2383703/71‑1 (MQ=255)
                  gcaggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCgg                                                             <  1:2676621/57‑1 (MQ=255)
                   caggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCt                                               <  1:3966419/70‑1 (MQ=255)
                   caggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTc                                              <  1:3087548/71‑1 (MQ=255)
                    aggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTc                                              >  1:1802319/1‑70 (MQ=255)
                      gTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGc                                                            <  1:3683052/54‑1 (MQ=255)
                      gTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGc                                                            <  1:79799/54‑1 (MQ=255)
                      gTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGc                                                            <  1:682337/54‑1 (MQ=255)
                       ttAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGc                                          <  1:946669/71‑1 (MQ=255)
                        tAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGc                                          >  1:3221188/1‑70 (MQ=255)
                         aaTGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCt                                                                  <  1:2416696/45‑1 (MQ=255)
                         aaTGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCg                                         <  1:198773/70‑1 (MQ=255)
                           tGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGt                                                          <  1:1356511/51‑1 (MQ=255)
                               ccaccCCAGCAGCAGACGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGa                                        >  1:1662666/1‑65 (MQ=255)
                                       gcagcagcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCa                           >  1:633231/1‑70 (MQ=255)
                                           cagcCGTTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                      <  1:3365430/71‑1 (MQ=255)
                                           cagcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                      <  1:340133/71‑1 (MQ=255)
                                            agcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                      >  1:375547/1‑70 (MQ=255)
                                            agcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                      <  1:4374201/70‑1 (MQ=255)
                                            agcCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                      >  1:2410122/1‑70 (MQ=255)
                                               cGGTATCCAGCGCGTAAATACCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATg                  >  1:657565/1‑71 (MQ=255)
                                                 gTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGc                 >  1:1456756/1‑70 (MQ=255)
                                                 gTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCt                >  1:1677056/1‑71 (MQ=255)
                                                 gTATCCAGCGCGTATATACCTTCAGGCGTACCTTTGCCCTCCAGCGa                                        >  1:405343/1‑47 (MQ=255)
                                                  tATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg               <  1:717163/71‑1 (MQ=255)
                                                     ccAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAt            >  1:3211087/1‑71 (MQ=255)
                                                      cAGCGCGTAAATACCTTCCGGCGTACCTTTGCCTTcc                                             >  1:740868/1‑37 (MQ=37)
                                                       aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTg                               <  1:449598/50‑1 (MQ=255)
                                                       aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCaa                          <  1:3650893/55‑1 (MQ=255)
                                                       aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCaa                          <  1:1482152/55‑1 (MQ=255)
                                                       aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAttc          >  1:4336703/1‑71 (MQ=255)
                                                         cgcgTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTAttctt        <  1:4341370/71‑1 (MQ=255)
                                                            gTAAATACCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGc                      <  1:19926/54‑1 (MQ=39)
                                                            gTAAATACCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGc                      <  1:762887/54‑1 (MQ=39)
                                                             tatataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGc                      <  1:527912/53‑1 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCGTCCCGCGATGCCCAGTGACCAAAGGCGa                    >  1:408695/1‑53 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGtcc                                   <  1:3203038/38‑1 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCa                           <  1:467169/46‑1 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTa             <  1:4418000/60‑1 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTa             <  1:714358/60‑1 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa  >  1:4415008/1‑71 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa  >  1:2265887/1‑71 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa  >  1:309818/1‑71 (MQ=255)
                                                               tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGa  >  1:2920827/1‑71 (MQ=255)
                                                               |                                                                      
TTTATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGA  >  W3110S.gb/50440‑50573

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: