Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,989,605 G→A 16.0% intergenic (+40/+39) yecI → / ← yecJ predicted ferritin‑like protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,989,6050GA16.0% 59.7 / 5.1 25intergenic (+40/+39)yecI/yecJpredicted ferritin‑like protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/13);  new base A (2/2);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTG  >  W3110S.gb/1989539‑1989674
                                                                  |                                                                     
aaTGTCGTGTCACACCAGCTTCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTgag                                                                     <  1:1255453/69‑1 (MQ=255)
aaTGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTgag                                                                     <  1:406778/69‑1 (MQ=255)
aaTGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTaaga                                                                    <  1:3047300/70‑1 (MQ=255)
    tCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTAAGAAAcc                                                                <  1:696914/70‑1 (MQ=255)
      gtgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCt                                                             <  1:3589646/71‑1 (MQ=255)
               cAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCt                                                    <  1:2958615/71‑1 (MQ=255)
                        tGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGtt                                                 >  1:3393137/1‑65 (MQ=255)
                         gATCATCATCGGCGCTAATGCATTGTGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAagc                                           <  1:1013049/70‑1 (MQ=255)
                         gATCATCATCGGCGCTAATGCATTGTGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAagc                                           <  1:1008179/70‑1 (MQ=255)
                         gATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAag                                            >  1:4058984/1‑69 (MQ=255)
                               catcGGCGCTAATGCATTGCGCCGATGAAGGTTTTAAGAAACCGCTGCCt                                                         >  1:1370577/1‑50 (MQ=255)
                               catcGGCGCTAATGCATTGCGCCGATGAAGGTTTTAAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtt                                    >  1:2477422/1‑71 (MQ=255)
                                 tcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttt                                  <  1:395522/71‑1 (MQ=255)
                                    gcgcTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttttt                               <  1:1573748/71‑1 (MQ=255)
                                          aTGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATgg                          <  1:1917662/70‑1 (MQ=255)
                                           tGCATTGCTCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGtt                                                 >  1:4051256/1‑46 (MQ=255)
                                             cATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAag                                            >  1:2474934/1‑49 (MQ=255)
                                             cATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAtt                      <  1:2937838/71‑1 (MQ=255)
                                             cATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAtt                      <  1:1314361/71‑1 (MQ=255)
                                             cATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAt                       >  1:501582/1‑70 (MQ=255)
                                             cATTGCGCCGATGAAGGTTTTAAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAtt                      <  1:1239570/71‑1 (MQ=255)
                                               ttGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATggg                         <  1:1570006/66‑1 (MQ=255)
                                               ttGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATggg                         <  1:2511714/66‑1 (MQ=255)
                                                 gcgcCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAt                       >  1:2281078/1‑66 (MQ=255)
                                                 gcgcCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAt                       >  1:2276417/1‑66 (MQ=255)
                                                   gcCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATgg                          <  1:2545760/61‑1 (MQ=255)
                                                          aaGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttttt                               <  1:2008781/49‑1 (MQ=255)
                                                                ttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACtt   >  1:3728501/1‑71 (MQ=255)
                                                                  gagaAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTg  <  1:3543444/70‑1 (MQ=255)
                                                                  gagaAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTg  <  1:2121473/70‑1 (MQ=255)
                                                                  |                                                                     
AATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTG  >  W3110S.gb/1989539‑1989674

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: