Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,002,183 G→A 15.2% R114C (CGC→TGC)  fliY ← cystine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,002,1830GA15.2% 79.1 / 7.4 33R114C (CGC→TGC) fliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base A (3/2);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

ACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAT  >  W3110S.gb/2002115‑2002249
                                                                    |                                                                  
aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCt                                                                  >  1:1469476/1‑71 (MQ=255)
             aGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGt                                                     >  1:2498389/1‑71 (MQ=255)
             aGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGt                                                     <  1:1005167/71‑1 (MQ=255)
              gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAGTCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                                                    <  1:963999/71‑1 (MQ=255)
              gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                                                    <  1:2420726/71‑1 (MQ=255)
              gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                                                    <  1:16508/71‑1 (MQ=255)
                gcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                                  >  1:2087256/1‑71 (MQ=255)
                       tACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGAtt                                            >  1:336634/1‑70 (MQ=255)
                         ccAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                         >  1:1526362/1‑71 (MQ=255)
                         ccAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                         >  1:1017448/1‑71 (MQ=255)
                               aTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATgg                                                      >  1:1013471/1‑52 (MQ=255)
                               aTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                                >  1:621013/1‑58 (MQ=255)
                               aTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                               >  1:1758264/1‑59 (MQ=255)
                                     tACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                              <  1:1459398/70‑1 (MQ=255)
                                       cGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                              >  1:1661577/1‑52 (MQ=255)
                                       cGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                              >  1:2742675/1‑52 (MQ=255)
                                       cGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAc                           >  1:2523495/1‑71 (MQ=255)
                                       cGGGGTTGAGAAATCGTATTTTTTCTTGCCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACg                          >  1:941743/1‑71 (MQ=255)
                                            ttGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa                      <  1:905968/71‑1 (MQ=255)
                                            ttGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa                      <  1:1739272/71‑1 (MQ=255)
                                            ttGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGAtt                                            <  1:2147069/49‑1 (MQ=255)
                                             tGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa                      <  1:2586069/70‑1 (MQ=255)
                                              gagaAATCGTATTTTTTCTTGCACTCATCAGAAAt                                                        >  1:1058327/1‑35 (MQ=255)
                                               agaAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa                                          <  1:1538317/48‑1 (MQ=255)
                                               agaAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa                                          <  1:1323004/48‑1 (MQ=255)
                                                 aaaTCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc                 >  1:286807/1‑71 (MQ=255)
                                                  aaTCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc                <  1:721321/71‑1 (MQ=255)
                                                     cGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAg              >  1:1208123/1‑70 (MQ=255)
                                                      gTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGc             <  1:2397162/70‑1 (MQ=255)
                                                      gTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGc             >  1:851107/1‑70 (MQ=255)
                                                      gTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGc             >  1:2482515/1‑70 (MQ=255)
                                                         ttttttCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                              <  1:1964040/50‑1 (MQ=255)
                                                           ttttCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                            >  1:1243142/1‑50 (MQ=255)
                                                                  gcgcTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAt  <  1:1774559/69‑1 (MQ=255)
                                                                    |                                                                  
ACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCAT  >  W3110S.gb/2002115‑2002249

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: