Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 45,813 C→T 100% P3S (CCG→TCG)  yaaU → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb45,8130CT89.2% 96.1 / 4.6 37P3S (CCG→TCG) yaaUpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (18/15);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGC  >  W3110S.gb/45748‑45878
                                                                 |                                                                 
tGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTc                                                               >  1:1371495/1‑70 (MQ=255)
         tttGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACttt                                                     <  1:2429980/71‑1 (MQ=255)
              tACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGa                                                >  1:3404325/1‑71 (MQ=255)
              tACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGa                                                >  1:2966664/1‑71 (MQ=255)
                           aCTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTctcc                                   >  1:1119260/1‑71 (MQ=255)
                           aCTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTctcc                                   >  1:2908063/1‑71 (MQ=255)
                               aCTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTctcctc                                 <  1:548958/69‑1 (MQ=255)
                                   tcATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTAtt                             <  1:788542/69‑1 (MQ=255)
                                   tcATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTAtt                             <  1:3456883/69‑1 (MQ=255)
                                   tcATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTAtt                             <  1:338264/69‑1 (MQ=255)
                                     aTTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:2678875/68‑1 (MQ=255)
                                     aTTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:2649885/68‑1 (MQ=255)
                                      ttGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCa                                           <  1:527915/52‑1 (MQ=255)
                                        gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGAc                                                  >  1:1433577/1‑43 (MQ=255)
                                        gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                      >  1:901414/1‑71 (MQ=255)
                                        gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                      >  1:274505/1‑71 (MQ=255)
                                        gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                      >  1:1986016/1‑71 (MQ=255)
                                        gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                      >  1:1070718/1‑71 (MQ=255)
                                        gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                      >  1:1438882/1‑71 (MQ=255)
                                                cAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTAtt                             <  1:1480250/56‑1 (MQ=255)
                                                 aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc             <  1:853548/71‑1 (MQ=255)
                                                 aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc             <  1:2895423/71‑1 (MQ=255)
                                                    aaTAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCtgt          <  1:496115/71‑1 (MQ=255)
                                                     aTAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCtgt          <  1:745539/70‑1 (MQ=255)
                                                     aTAACAATGCAATAGTCCAGAAAATTTGACGATCTCAAATTATCCTCTa                               >  1:1407605/1‑49 (MQ=25)
                                                      tAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTgg        >  1:810732/1‑71 (MQ=255)
                                                      tAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTgg        >  1:2698076/1‑71 (MQ=255)
                                                      tAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTgg        >  1:1239448/1‑71 (MQ=255)
                                                        aCCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTggg       >  1:102871/1‑70 (MQ=255)
                                                        aCCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTggg       >  1:2014580/1‑70 (MQ=255)
                                                         ccATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:37334/48‑1 (MQ=255)
                                                         ccATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:3643671/48‑1 (MQ=255)
                                                         ccATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:1501016/48‑1 (MQ=255)
                                                         ccATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:1689032/48‑1 (MQ=255)
                                                         ccATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTc                            <  1:2606409/48‑1 (MQ=255)
                                                            tGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGc  >  1:1085144/1‑71 (MQ=255)
                                                            tGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGc  >  1:2163553/1‑71 (MQ=255)
                                                                 |                                                                 
TGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGC  >  W3110S.gb/45748‑45878

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: