Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,682,955 G→A 8.6% R673R (CGC→CGT yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,682,9550GA 8.6% 179.2 / 5.9 58R673R (CGC→CGTyfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base A (2/3);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

AGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCTTGCGCATCGCTATAAAC  >  Exported/1682890‑1683023
                                                                 |                                                                    
aGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGtt                                                                  <  1:1387218/70‑1 (MQ=255)
aGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTc                                                                 <  1:2369328/71‑1 (MQ=255)
aGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTc                                                                 <  1:1734531/71‑1 (MQ=255)
 gATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGtt                                                                  >  1:1648911/1‑69 (MQ=255)
 gATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTcg                                                                >  1:1327482/1‑71 (MQ=255)
 gATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTc                                                                 >  1:3499686/1‑70 (MQ=255)
  aTCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTcgc                                                               >  1:1388647/1‑71 (MQ=255)
  aTCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTcgc                                                               >  1:2531200/1‑71 (MQ=255)
    cTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTcgc                                                               >  1:1593027/1‑69 (MQ=255)
       tttGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAAt                                                           >  1:3487150/1‑70 (MQ=255)
        ttGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATc                                                          <  1:3828676/70‑1 (MQ=255)
               aaCTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCAc                                                  <  1:2872646/71‑1 (MQ=255)
                 cTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCa                                                   <  1:2292902/68‑1 (MQ=255)
                  tgtgACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGaa                                                            >  1:1949799/1‑58 (MQ=255)
                  tgtgACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATc                                                          >  1:323800/1‑60 (MQ=255)
                  tgtgACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATc                                                          >  1:2297311/1‑60 (MQ=255)
                  tgtgACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATc                                                          >  1:1770714/1‑60 (MQ=255)
                  tgtgACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAgg                                                       <  1:1725243/63‑1 (MQ=255)
                    tgACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTg                                               <  1:1502088/69‑1 (MQ=255)
                        tGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGgcgc                                                    <  1:2139421/60‑1 (MQ=255)
                        tGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGgcgc                                                    <  1:3318919/60‑1 (MQ=255)
                        tGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGgcgc                                                    <  1:282752/60‑1 (MQ=255)
                        tGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGgcgc                                                    <  1:2166296/60‑1 (MQ=255)
                        tGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGcc                                         >  1:336674/1‑71 (MQ=255)
                         gCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGccc                                        >  1:2771720/1‑71 (MQ=255)
                          ccAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGccc                                        <  1:230868/70‑1 (MQ=255)
                          ccAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGccc                                        <  1:394350/70‑1 (MQ=255)
                           cAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGACGTTCGCGAAt                                                           <  1:2724926/50‑1 (MQ=255)
                           cAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGACGTTCGCGAAt                                                           <  1:15625/50‑1 (MQ=255)
                           cAGCCTTCATCTTCTGACCAGATCCAGTAGTAATCGCGACGTTCGCGAAt                                                           <  1:674248/50‑1 (MQ=255)
                             gCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGaa                                                            <  1:2841289/47‑1 (MQ=255)
                             gCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAg                                                        <  1:3561060/51‑1 (MQ=255)
                             gCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGACGTTCGCGAATCa                                                         >  1:2568145/1‑50 (MQ=255)
                             gCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGACGTTCGCGAATCa                                                         >  1:3416583/1‑50 (MQ=255)
                              ccTTCATCTTCTGACCAGTTCCAGTAGTAATCGCTGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGacac                                   <  1:3782329/71‑1 (MQ=255)
                               cTTCATCTTCTGACCAGTTCCAGTAGTAATCGCTGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGacac                                   <  1:3777918/70‑1 (MQ=255)
                                ttCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTcgcg                                                              <  1:879844/42‑1 (MQ=255)
                                ttCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTcgcg                                                              <  1:304555/42‑1 (MQ=255)
                                ttCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTc                                                                 <  1:623332/39‑1 (MQ=255)
                                  cATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATc                                                          >  1:1049008/1‑44 (MQ=255)
                                  cATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCt                               >  1:288989/1‑71 (MQ=255)
                                  cATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCt                               >  1:2540663/1‑71 (MQ=255)
                                   aTCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCtt                              >  1:1599442/1‑71 (MQ=255)
                                     cttctGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCtttt                            >  1:8107/1‑71 (MQ=255)
                                      ttctGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCtttt                            >  1:1217471/1‑70 (MQ=255)
                                        ctGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCtttttt                          <  1:272888/70‑1 (MQ=255)
                                         tGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCa                        >  1:1013766/1‑71 (MQ=255)
                                         tGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCa                        >  1:2131126/1‑71 (MQ=255)
                                         tGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCa                        >  1:2156869/1‑71 (MQ=255)
                                             cAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGc                     <  1:3118069/70‑1 (MQ=255)
                                              aGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCt                   <  1:3568376/71‑1 (MQ=255)
                                              aGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCt                   <  1:2730330/71‑1 (MQ=255)
                                                ttCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTg                                               >  1:2482827/1‑41 (MQ=255)
                                                ttCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTg                                               >  1:266128/1‑41 (MQ=255)
                                                ttCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACg                      <  1:485366/66‑1 (MQ=255)
                                                  ccAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCtttttt                          >  1:34819/1‑60 (MQ=255)
                                                           aTCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCTTGCGCATCGCta       >  1:2095056/1‑70 (MQ=255)
                                                           aTCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCTTGCGCATCGCta       >  1:577868/1‑70 (MQ=255)
                                                              gcggcgTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCTTGCGCATCGCTATaa    <  1:3739288/70‑1 (MQ=255)
                                                                ggcgTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCTTGCGCATCGCTATAAAc  <  1:3260098/70‑1 (MQ=255)
                                                                 gcgTTCGCGAATCAGGAGCACCTGCAAGCCCGACACGGCttt                             >  1:392478/1‑42 (MQ=255)
                                                                 |                                                                    
AGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCTTTTTTCACGCCTTGCGCATCGCTATAAAC  >  Exported/1682890‑1683023

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: