Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,002,152 T→C 17.1% Y124C (TAC→TGC)  fliY ← cystine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,002,1520TC17.1% 91.9 / 14.1 41Y124C (TAC→TGC) fliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (12/22);  new base C (3/4);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAA  >  W3110S.gb/2002091‑2002221
                                                             |                                                                     
ggCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGtt                                                               >  1:2690791/1‑70 (MQ=255)
       ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATc                                                       >  1:2510370/1‑71 (MQ=255)
       ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGCACGGGGTTGAGAAATc                                                       >  1:1942367/1‑71 (MQ=255)
              tGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTAtttt                                                >  1:360594/1‑71 (MQ=255)
               gCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTa                                                                      >  1:12417/1‑48 (MQ=255)
               gCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTa                                                                      >  1:2533469/1‑48 (MQ=255)
                  ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGt                                                     <  1:4659794/62‑1 (MQ=255)
                  ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCt                                            >  1:620884/1‑71 (MQ=255)
                  ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCt                                            >  1:614780/1‑71 (MQ=255)
                   tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTTTACGGGGTTgag                                                            <  1:2474374/54‑1 (MQ=255)
                    cGTTACCTTTTTTCACCAGCGCCTGAATACCAGGAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTg                                          >  1:4610798/1‑71 (MQ=255)
                    cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCt                                            <  1:3210749/69‑1 (MQ=255)
                     gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTg                                          <  1:1067354/70‑1 (MQ=255)
                        aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGCACgggg                                                                 >  1:59208/1‑44 (MQ=255)
                         ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCt                                      <  1:2401977/70‑1 (MQ=255)
                         ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCt                                      <  1:2441586/70‑1 (MQ=255)
                         ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCt                                      <  1:3033182/70‑1 (MQ=255)
                          cTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACgg                                                                   <  1:1897341/40‑1 (MQ=255)
                           ttttttCACCAGCGCCTGAATACCAGAAATGGTGCAc                                                                     <  1:4206555/37‑1 (MQ=255)
                           ttttttCACCAGCGCCTGAATACCAGAAATGGTGCAc                                                                     <  1:1494050/37‑1 (MQ=255)
                            tttttCACCAGCGCCTGAATACCAGAAATGGTGCACGGGGTTGAGaa                                                          <  1:4196608/47‑1 (MQ=255)
                             ttttCACCAGCGCCTGAATACCAGAAATGGTGCACgg                                                                   <  1:1350578/37‑1 (MQ=255)
                                tCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                               <  1:4151604/70‑1 (MQ=255)
                                tCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                               <  1:189386/70‑1 (MQ=255)
                                tCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                               <  1:612258/70‑1 (MQ=255)
                                      gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                        <  1:450242/71‑1 (MQ=255)
                                      gcgcCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                        <  1:4364269/71‑1 (MQ=255)
                                           tGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                           <  1:3354632/63‑1 (MQ=255)
                                             aaTACCAGAAATGGTGTACGGGGTTGAGAAATCGTAtttttt                                              >  1:1875984/1‑42 (MQ=255)
                                               tACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTa               <  1:4873449/71‑1 (MQ=255)
                                                aCCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa              >  1:4487938/1‑71 (MQ=255)
                                                aCCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa              >  1:3686484/1‑71 (MQ=255)
                                                 ccAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa              <  1:805451/70‑1 (MQ=255)
                                                 ccAGAAATGGTGCACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt             >  1:2366987/1‑71 (MQ=255)
                                                  cAGAAATGGTGTACGGGGTTGAGAAATCGTAtttttt                                              >  1:639839/1‑37 (MQ=255)
                                                    gAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                               <  1:1902746/50‑1 (MQ=255)
                                                    gAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac          <  1:2671161/71‑1 (MQ=255)
                                                    gAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac          <  1:2788466/71‑1 (MQ=255)
                                                     aaaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                               <  1:1533078/49‑1 (MQ=255)
                                                     aaaTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGa                               <  1:2596850/49‑1 (MQ=255)
                                                        tGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac          <  1:2134688/67‑1 (MQ=255)
                                                          gtgtACTGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc    <  1:3310097/71‑1 (MQ=255)
                                                          gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                   >  1:3271657/1‑40 (MQ=255)
                                                          gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                   >  1:2341220/1‑40 (MQ=255)
                                                          gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                   >  1:3767032/1‑40 (MQ=255)
                                                          gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaa                              >  1:1329143/1‑45 (MQ=255)
                                                          gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc    <  1:1438624/71‑1 (MQ=255)
                                                          gtgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc    <  1:3372983/71‑1 (MQ=255)
                                                           tgtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                     <  1:5046609/37‑1 (MQ=255)
                                                            gtACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa  >  1:5005747/1‑71 (MQ=255)
                                                             tACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa  <  1:2911981/70‑1 (MQ=255)
                                                             tACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa  <  1:4699225/70‑1 (MQ=255)
                                                             tACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa  <  1:4205545/70‑1 (MQ=255)
                                                             |                                                                     
GGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAA  >  W3110S.gb/2002091‑2002221

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: