Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,590,693 A→G 16.0% L341L (CTT→CTC ydeT ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,590,6930AG16.0% 56.4 / 5.5 25L341L (CTT→CTCydeThypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base G (2/2);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCGTGTGTGACAA  >  W3110S.gb/1590627‑1590763
                                                                  |                                                                      
gTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGGAGTc                                                                    >  1:477091/1‑71 (MQ=255)
          aTACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATa                                                          <  1:1072644/71‑1 (MQ=255)
              aGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAAcc                                                      >  1:136766/1‑71 (MQ=255)
               gTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCt                                                     <  1:1549907/71‑1 (MQ=255)
                   gAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGAc                                                  >  1:1499381/1‑70 (MQ=255)
                        tttATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAAc                                                       >  1:1436505/1‑60 (MQ=255)
                         ttATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTc                                           >  1:949778/1‑71 (MQ=255)
                           aTCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTc                                           >  1:1332898/1‑69 (MQ=255)
                              tttGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGa                                      >  1:203598/1‑71 (MQ=255)
                                tGCCCCATGAAACCTGTAATTGCCCTGACTGTGGGAGTCCAGTCAGATAAACCTGACCATTTTCCGCGAc                                     >  1:1210881/1‑70 (MQ=255)
                                tGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCcgcg                                       <  1:199992/68‑1 (MQ=255)
                                tGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCcgc                                        <  1:242396/67‑1 (MQ=255)
                                    ccATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAAtgc                                <  1:174344/71‑1 (MQ=255)
                                        gAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGAc                                                  <  1:107305/49‑1 (MQ=255)
                                            ccTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCAttt                        >  1:953248/1‑71 (MQ=255)
                                                   ttGCCCTGACTGTGGGAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACaa                                   <  1:827175/53‑1 (MQ=255)
                                                   ttGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAAtg                                 <  1:1248142/55‑1 (MQ=255)
                                                    tGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGa                                                   >  1:639459/1‑36 (MQ=255)
                                                     gCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAAtgct                               >  1:1215027/1‑55 (MQ=255)
                                                        cTGACTGTGGGAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACaa                                   <  1:138584/48‑1 (MQ=255)
                                                        cTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCg            <  1:1113316/71‑1 (MQ=255)
                                                        cTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCg            <  1:398063/71‑1 (MQ=255)
                                                        cTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCg            >  1:124526/1‑71 (MQ=255)
                                                         tGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTTCCATTTTTATTCTCTCCgt           >  1:1230949/1‑71 (MQ=255)
                                                                 gAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCGTGTGTGAc    <  1:867034/70‑1 (MQ=255)
                                                                  aaGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCGTGTGTGACaa  >  1:1025636/1‑71 (MQ=255)
                                                                  |                                                                      
GTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAACCTGACCATTTTCCGCGACAATGCTGCCATTTTTATTCTCTCCGTGTGTGACAA  >  W3110S.gb/1590627‑1590763

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: