Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,130,490 A→G 9.4% F192S (TTT→TCT)  yghQ ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,130,4900AG 9.4% 138.1 / 5.4 53F192S (TTT→TCT) yghQpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base G (2/3);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TAATTCGCGCGCGGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGC  >  W3110S.gb/3130420‑3130560
                                                                      |                                                                      
tAATTCGCGCGCGGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACaa                                                                        <  1:3705761/71‑1 (MQ=255)
 aaTTCGCGCGCGGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACaa                                                                        >  1:3654197/1‑70 (MQ=255)
  aTTCGCGCGCGGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACaaaa                                                                      >  1:2444624/1‑71 (MQ=255)
     cgcgcgcgGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAAccc                                                                   >  1:3758043/1‑71 (MQ=255)
     cgcgcgcgGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAAccc                                                                   >  1:2874031/1‑71 (MQ=255)
     cgcgcgcgGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAAccc                                                                   >  1:913580/1‑71 (MQ=255)
     cgcgcgcgGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAAccc                                                                   >  1:3481945/1‑71 (MQ=255)
     cgcgcgcgGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAAccc                                                                   >  1:750562/1‑71 (MQ=255)
         cgcgGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCaa                                                               <  1:522437/71‑1 (MQ=255)
              cAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCa                                                          <  1:1444364/71‑1 (MQ=255)
              cAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCa                                                          <  1:486474/71‑1 (MQ=255)
                aaCCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCaaa                                                        <  1:1746334/71‑1 (MQ=255)
                aaCCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCaaa                                                        <  1:3672380/71‑1 (MQ=255)
                  ccaccaGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACaaaa                                                                      >  1:2956370/1‑55 (MQ=255)
                  ccaccaGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGt                                                      >  1:34436/1‑71 (MQ=255)
                   caccacTACATGGTGCCGCCAACCAGATTCGACACGTACCACGCAATAACAAAACCCg                                                                  >  1:2623602/1‑58 (MQ=255)
                   caccaGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGt                                                      >  1:883698/1‑70 (MQ=255)
                   caccaGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGt                                                      >  1:1231577/1‑70 (MQ=255)
                   caccaGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGt                                                      >  1:3669816/1‑70 (MQ=255)
                     ccaGTACATGGTGCCGCCAACCAGATTCGACAGGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCaa                                                   <  1:1782193/71‑1 (MQ=255)
                     ccaGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCTCAAAACCAAAGTCaa                                                   <  1:3101061/71‑1 (MQ=255)
                         tACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCaaaa                                                 <  1:1369292/69‑1 (MQ=255)
                          aCATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATAc                                              >  1:1784809/1‑71 (MQ=255)
                           cATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATAc                                              >  1:2059128/1‑70 (MQ=255)
                           cATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATAc                                              >  1:3571883/1‑70 (MQ=255)
                           cATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATAc                                              >  1:3169259/1‑70 (MQ=255)
                               gTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCaaaa                                                 >  1:2361499/1‑63 (MQ=255)
                                  ccgccAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTac                                      >  1:3201242/1‑71 (MQ=255)
                                  ccgccAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTac                                      <  1:3715542/71‑1 (MQ=255)
                                  ccgccAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTac                                      <  1:341042/71‑1 (MQ=255)
                                    gccAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCt                                        <  1:3361885/67‑1 (MQ=255)
                                     ccACCCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGTTacgac                                   >  1:1010381/1‑71 (MQ=255)
                                     ccAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGTTacgac                                   >  1:2570063/1‑71 (MQ=255)
                                          cAGATTCGACACGTACCAGGCAATAACAAAACCCGCaaaa                                                             >  1:2077684/1‑40 (MQ=255)
                                             aTTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCaa                                                   <  1:1758796/47‑1 (MQ=255)
                                             aTTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCaa                                                   <  1:3585140/47‑1 (MQ=255)
                                             aTTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCaa                                                   <  1:135255/47‑1 (MQ=255)
                                              ttCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAg                                           >  1:2293637/1‑54 (MQ=255)
                                               tCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCt                         >  1:1214401/1‑71 (MQ=255)
                                               tCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCt                         >  1:2328547/1‑71 (MQ=255)
                                                    acGTACCAGGCAATAACAGAACCCGCAAAACCAAAGTCAAAAt                                                <  1:1045707/43‑1 (MQ=255)
                                                    acGTACCAGGCAATAACAGAACCCGCAAAACCAAAGTCAAAAt                                                <  1:925462/43‑1 (MQ=255)
                                                    acGTACCAGGCAATAACAGAACCCGCAAAACCAAAGTCAAAAt                                                <  1:2398262/43‑1 (MQ=255)
                                                    acGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTgcgcg                    >  1:2248078/1‑71 (MQ=255)
                                                      gTACCAGGCAATAACAGAACCCGCAAAACCAAAGTCAAAATa                                               >  1:766089/1‑42 (MQ=255)
                                                      gTACCAGGCAATAACAGAACCCGCAAAACCAAAGTCAAAATAc                                              >  1:2957267/1‑43 (MQ=255)
                                                             gCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAgg           <  1:3717920/71‑1 (MQ=255)
                                                              cAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGt          >  1:2366384/1‑71 (MQ=255)
                                                               aaTAACAAAACCCGCAAAACCAAAGTCAAAATACCa                                            <  1:1652142/36‑1 (MQ=255)
                                                                aTAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGtt         >  1:2712992/1‑70 (MQ=255)
                                                                 tAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGaaaa             <  1:1614581/65‑1 (MQ=255)
                                                                 tAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGaaaa             <  1:3827485/65‑1 (MQ=255)
                                                                 tAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGa                <  1:3556115/62‑1 (MQ=255)
                                                                  aaCAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGt          >  1:2944464/1‑67 (MQ=255)
                                                                   aCAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTacga                                    >  1:868041/1‑40 (MQ=255)
                                                                   aCAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGtt         >  1:711048/1‑67 (MQ=255)
                                                                    cAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTtcgtc    >  1:3339688/1‑71 (MQ=255)
                                                                      aaaCCCGCAAAACCAAAGTCAAAATACCAGGCTacga                                    >  1:2966914/1‑37 (MQ=255)
                                                                      aaaCCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGc  >  1:92279/1‑71 (MQ=255)
                                                                      |                                                                      
TAATTCGCGCGCGGCAAACCACCAGTACATGGTGCCGCCAACCAGATTCGACACGTACCAGGCAATAACAAAACCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGC  >  W3110S.gb/3130420‑3130560

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑