Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,998,196 A→G 13.2% Y454H (TAT→CAT)  yhiP ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,998,1960AG13.2% 130.3 / 13.1 53Y454H (TAT→CAT) yhiPpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base G (4/3);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGCACCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCC  >  W3110S.gb/3998126‑3998265
                                                                      |                                                                     
ggCACCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                       <  1:2492293/71‑1 (MQ=255)
   aCCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGAc                                                                    <  1:2373233/71‑1 (MQ=255)
      aTCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACggg                                                                 >  1:427950/1‑71 (MQ=255)
       tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                       <  1:1117831/64‑1 (MQ=255)
       tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                       <  1:200305/64‑1 (MQ=255)
       tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                                <  1:2384354/71‑1 (MQ=255)
       tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                                <  1:71217/71‑1 (MQ=255)
        cAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                                >  1:2371781/1‑70 (MQ=255)
         aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCaa                                                              <  1:194714/71‑1 (MQ=255)
         aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCaa                                                              <  1:500804/71‑1 (MQ=255)
         aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCaa                                                              <  1:2414662/71‑1 (MQ=255)
         aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGa                                                                     >  1:2724321/1‑64 (MQ=255)
         aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGa                                                                     >  1:436185/1‑64 (MQ=255)
         aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCaa                                                              >  1:1401279/1‑71 (MQ=255)
               aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGtctc                                                        <  1:889782/71‑1 (MQ=255)
               aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCATACCAAACACGTTGGTATAGACGGGCAACGtct                                                         >  1:1211823/1‑70 (MQ=255)
                  actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                       >  1:1906630/1‑53 (MQ=255)
                  actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                     >  1:216931/1‑71 (MQ=255)
                  actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                     >  1:1423069/1‑71 (MQ=255)
                  actactGCAACGCCCAGCGTGACCAGA‑CAATCTTACCAAACACGTTGGTATAGACGGG‑AACGTCTCAAGCg                                                   >  1:2486797/1‑71 (MQ=255)
                   ctactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                                <  1:1012414/59‑1 (MQ=255)
                      ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                  <  1:2655585/70‑1 (MQ=255)
                      ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                  <  1:1407235/70‑1 (MQ=255)
                      ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGa                                                 >  1:2624921/1‑71 (MQ=255)
                      ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGa                                                 >  1:1265776/1‑71 (MQ=255)
                           aCGTCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                            <  1:1869772/71‑1 (MQ=255)
                           aCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                            <  1:1077460/71‑1 (MQ=255)
                           aCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                            <  1:132059/71‑1 (MQ=255)
                            cGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTa                                           <  1:1056523/71‑1 (MQ=255)
                               ccAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                       <  1:1826789/40‑1 (MQ=255)
                                  gCGTGACCAGTCCAATCTTACCAAACACGTTGGtata                                                                       <  1:1664291/37‑1 (MQ=255)
                                  gCGTGACCAGACCAATCTTACCAAAGACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                     <  1:1288685/71‑1 (MQ=255)
                                  gCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                     <  1:774786/71‑1 (MQ=255)
                                  gCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                     <  1:408818/71‑1 (MQ=255)
                                    gTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                     >  1:428921/1‑53 (MQ=255)
                                    gTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttgt                                    >  1:843255/1‑70 (MQ=255)
                                    gTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttgt                                    >  1:2147194/1‑70 (MQ=255)
                                     tGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                     <  1:496949/68‑1 (MQ=255)
                                     tGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTcc                                  <  1:265104/71‑1 (MQ=255)
                                       accagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                                >  1:664181/1‑71 (MQ=255)
                                         cagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCg                                              >  1:340662/1‑55 (MQ=255)
                                              caATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTAGTCCGGCACcgcg                         >  1:1015672/1‑71 (MQ=255)
                                                 tCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                            <  1:46528/49‑1 (MQ=255)
                                                 tCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                                <  1:2375537/61‑1 (MQ=255)
                                                  cTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttgt                                    >  1:2720286/1‑56 (MQ=255)
                                                   ttACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                    >  1:1513900/1‑71 (MQ=255)
                                                   ttACCAAACACGTTGGTATAGACGGGCAACGtctc                                                        >  1:1978626/1‑35 (MQ=255)
                                                    tACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                     <  1:74674/53‑1 (MQ=255)
                                                     aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                    <  1:2649076/69‑1 (MQ=255)
                                                     aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                    <  1:1661337/69‑1 (MQ=255)
                                                     aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                      <  1:464348/67‑1 (MQ=255)
                                                      ccAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGTGGTaaa                    <  1:364844/68‑1 (MQ=255)
                                                           cacGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTaa                                          >  1:1049181/1‑41 (MQ=255)
                                                           cacGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTaa                                          >  1:647033/1‑41 (MQ=255)
                                                           cacGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca            <  1:2656616/71‑1 (MQ=255)
                                                            acGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca            >  1:1703497/1‑70 (MQ=255)
                                                                 ggTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCcaca          >  1:1387314/1‑67 (MQ=255)
                                                                  gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgcc     <  1:346995/71‑1 (MQ=255)
                                                                  gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgc      <  1:1091969/70‑1 (MQ=255)
                                                                     taGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgccgcc  <  1:1661757/71‑1 (MQ=255)
                                                                      |                                                                     
GGCACCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCC  >  W3110S.gb/3998126‑3998265

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