Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,591,647 T→C 7.3% intergenic (+106/+83) yfeS → / ← cysM conserved hypothetical protein/cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,591,6470TC 7.3% 177.0 / 4.6 68intergenic (+106/+83)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base T (36/27);  new base C (3/2);  total (39/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.91e-01

GTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCA  >  Exported/1591578‑1591709
                                                                     |                                                              
gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatata                                                               <  1:3988206/71‑1 (MQ=39)
 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatatat                                                              <  1:2293409/71‑1 (MQ=39)
  gcgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATAtt                                                             >  1:1995596/1‑71 (MQ=39)
     cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                           <  1:4304830/70‑1 (MQ=40)
        tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGa                                                        >  1:922555/1‑70 (MQ=255)
         gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaat                                                      >  1:2029427/1‑70 (MQ=11)
            tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                     >  1:2781268/1‑69 (MQ=255)
            tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt                                                   >  1:2820729/1‑71 (MQ=255)
            tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatca                                                   >  1:305977/1‑68 (MQ=11)
            tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatc                                                    <  1:4377490/70‑3 (MQ=11)
                   gATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGc                                            >  1:2859712/1‑71 (MQ=11)
                   gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg                                             >  1:161687/1‑70 (MQ=255)
                   gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc                                            >  1:3033291/1‑71 (MQ=255)
                   gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc                                            <  1:487477/71‑1 (MQ=255)
                    aTGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaa                                                       <  1:4827526/59‑1 (MQ=11)
                    aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                           >  1:1538705/1‑55 (MQ=255)
                     tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCt                                          <  1:4258487/71‑1 (MQ=255)
                      gcggcgTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTa                                               >  1:4840718/1‑65 (MQ=11)
                       cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                     >  1:5689638/1‑58 (MQ=255)
                          cttAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa                                     <  1:1974287/69‑1 (MQ=255)
                           gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                    <  1:291041/71‑1 (MQ=255)
                           gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                    >  1:1172367/1‑71 (MQ=255)
                            tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat                                                     <  1:5390949/53‑2 (MQ=11)
                            tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat                                                     <  1:3676233/53‑2 (MQ=11)
                            tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAg                                              <  1:1589502/60‑1 (MQ=11)
                            tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCt                                          <  1:3975788/64‑1 (MQ=11)
                            tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAACCCTTGt                                                <  1:3206/58‑1 (MQ=255)
                            tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                    <  1:579092/70‑1 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTg                                                            >  1:2682872/1‑44 (MQ=38)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                  >  1:3680747/1‑70 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                  >  1:3221703/1‑70 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                         >  1:2702139/1‑47 (MQ=11)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                         >  1:2013041/1‑47 (MQ=11)
                               aCGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTa                                <  1:2269605/71‑1 (MQ=11)
                                 gCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAgc                              >  1:2944531/1‑71 (MQ=11)
                                 gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaatt                                                     <  1:4338094/48‑3 (MQ=11)
                                   cTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcg                             >  1:234498/1‑70 (MQ=255)
                                   attATCCGGCCTACAAAATCGTGCTAATTCAATATATTGGAGAAACCTTg                                                 <  1:4279856/49‑1 (MQ=255)
                                     tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                               >  1:3572312/1‑50 (MQ=255)
                                     tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                          >  1:1372159/1‑71 (MQ=255)
                                     tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                          >  1:1676644/1‑71 (MQ=255)
                                     tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                          >  1:1725128/1‑71 (MQ=255)
                                      ttCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                       <  1:329634/40‑1 (MQ=11)
                                      aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCtt                                                  <  1:136632/46‑1 (MQ=255)
                                      aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                          >  1:1384668/1‑70 (MQ=255)
                                      aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                          >  1:596644/1‑70 (MQ=255)
                                       tCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCATCa                        <  1:4161649/71‑1 (MQ=11)
                                              tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTg                                         <  1:559899/47‑1 (MQ=255)
                                              tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcgc                            <  1:5312983/60‑1 (MQ=255)
                                                cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa                                     >  1:5266805/1‑49 (MQ=255)
                                                cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCa                           >  1:2144260/1‑59 (MQ=255)
                                                  aaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                 >  1:4591587/1‑67 (MQ=255)
                                                   aaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                  >  1:5019737/1‑49 (MQ=255)
                                                   aaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg            >  1:5153723/1‑71 (MQ=255)
                                                    aTCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                               >  1:2087793/1‑35 (MQ=255)
                                                        tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt       <  1:61967/71‑1 (MQ=255)
                                                        tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt       <  1:5566025/71‑1 (MQ=255)
                                                          cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca     <  1:1215325/71‑1 (MQ=255)
                                                          cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca     <  1:2307079/71‑1 (MQ=255)
                                                          cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc      <  1:1925143/70‑1 (MQ=255)
                                                           tAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                >  1:3941392/1‑59 (MQ=255)
                                                           tAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                >  1:3445288/1‑59 (MQ=255)
                                                             aTTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                               >  1:5526055/1‑42 (MQ=255)
                                                             aTTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcat    <  1:5241376/69‑1 (MQ=255)
                                                             aTTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcat    <  1:5247975/69‑1 (MQ=255)
                                                               tCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt               <  1:664481/56‑1 (MQ=255)
                                                                cAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa                                >  1:3956066/1‑38 (MQ=255)
                                                                cAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt               >  1:324939/1‑55 (MQ=255)
                                                                cAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                >  1:3762801/1‑54 (MQ=255)
                                                                cAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcat    >  1:324839/1‑66 (MQ=255)
                                                                    atatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt       <  1:1526085/59‑1 (MQ=255)
                                                                    atatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCACATCAGGCAGTTTTGCGTTTGt       <  1:672142/59‑1 (MQ=255)
                                                                     tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcatca  >  1:4745054/1‑63 (MQ=255)
                                                                     tatTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcatca  >  1:2110136/1‑63 (MQ=255)
                                                                     |                                                              
GTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCA  >  Exported/1591578‑1591709

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: