Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,274,249 A→G 100% T64A (ACT→GCT)  rtn → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,274,2490AG85.7% 89.7 / 7.3 35T64A (ACT→GCT) rtnconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base G (14/16);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCC  >  W3110S.gb/2274184‑2274302
                                                                 |                                                     
aaTACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCgctact                                                  <  1:260167/71‑1 (MQ=255)
  tACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGa                                                >  1:2041847/1‑71 (MQ=255)
           cTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCa                                       <  1:142273/71‑1 (MQ=255)
           cTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCa                                       <  1:2259543/71‑1 (MQ=255)
           cTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGATc                                         >  1:228043/1‑69 (MQ=255)
              aTTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGc                                     >  1:2481417/1‑70 (MQ=255)
               ttACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCg                                   >  1:2584946/1‑71 (MQ=255)
                    gACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGAc                              <  1:2705625/71‑1 (MQ=255)
                    gACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGAc                              <  1:2241724/71‑1 (MQ=255)
                        tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTa                          <  1:2687938/71‑1 (MQ=255)
                        tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTa                          <  1:2786241/71‑1 (MQ=255)
                        tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTa                          <  1:461205/71‑1 (MQ=255)
                        tACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTa                          <  1:1202206/71‑1 (MQ=255)
                            aaaGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATa                      <  1:2317546/71‑1 (MQ=255)
                            aaaGTATCTCGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATa                      >  1:801658/1‑71 (MQ=255)
                                  tctcGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGa                               <  1:2450105/56‑1 (MQ=255)
                                  tctcGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGa                               <  1:1544726/56‑1 (MQ=255)
                                  tctcGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGAt                       <  1:2504505/64‑1 (MQ=255)
                                    tcGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCa               >  1:867761/1‑70 (MQ=255)
                                    tcGATACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAg              >  1:65633/1‑71 (MQ=255)
                                    tcGATACCTATTTTGCCGACCTGAAATCCACTACTGAc                                               >  1:1754069/1‑38 (MQ=255)
                                    tcGATACCTATTTTGCCGACCTGAAATCCACTACTGAc                                               >  1:1689348/1‑38 (MQ=255)
                                    tcGATACCTATTTTGCCGACCTGAAATCCACTACTGAc                                               >  1:2699789/1‑38 (MQ=255)
                                       aTACCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTa                          <  1:2516789/56‑1 (MQ=255)
                                         aCCTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGa                        >  1:1247740/1‑56 (MQ=255)
                                          ccTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt        <  1:2413481/71‑1 (MQ=255)
                                          ccTATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt        <  1:22095/71‑1 (MQ=255)
                                            tATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCa            <  1:657502/65‑1 (MQ=255)
                                            tATTTTGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTaa      >  1:39929/1‑71 (MQ=255)
                                              ttttGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCt        <  1:1572438/67‑1 (MQ=255)
                                              ttttGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcc    >  1:843284/1‑71 (MQ=255)
                                               tttGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGc                                  >  1:2676656/1‑40 (MQ=255)
                                                ttGCCGACCTGAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcccc  >  1:1578765/1‑71 (MQ=255)
                                                          gAAATCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGa                        >  1:156973/1‑39 (MQ=255)
                                                           aaaTCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGa                        >  1:2456472/1‑38 (MQ=255)
                                                           aaaTCCGCTACTGACCGGCTCCAGCCGCTGACCTTAGa                        >  1:1591396/1‑38 (MQ=255)
                                                                 |                                                     
AATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCC  >  W3110S.gb/2274184‑2274302

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: