Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,591,660 T→A 11.9% intergenic (+119/+70) 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,6600TA11.9% 157.3 / 12.3 59intergenic (+119/+70)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base A (3/4);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAT  >  Exported/1591590‑1591723
                                                                      |                                                               
tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt                                                                 >  1:2820729/1‑71 (MQ=255)
       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)
         tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCt                                                        <  1:4258487/71‑1 (MQ=255)
          gcggcgTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTa                                                             >  1:4840718/1‑65 (MQ=11)
              cttAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa                                                   <  1:1974287/69‑1 (MQ=255)
               gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                                  <  1:291041/71‑1 (MQ=255)
               gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                                  >  1:1172367/1‑71 (MQ=255)
                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)
                  aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                                >  1:3680747/1‑70 (MQ=255)
                  aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                                >  1:3221703/1‑70 (MQ=255)
                   aCGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTa                                              <  1:2269605/71‑1 (MQ=11)
                     gCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAgc                                            >  1:2944531/1‑71 (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:1725128/1‑71 (MQ=255)
                         tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                                        >  1:1676644/1‑71 (MQ=255)
                         tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                                        >  1:1372159/1‑71 (MQ=255)
                          aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCtt                                                                <  1:136632/46‑1 (MQ=255)
                          aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                                        >  1:596644/1‑70 (MQ=255)
                          aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt                                        >  1:1384668/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:5566025/71‑1 (MQ=255)
                                            tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt                     <  1:61967/71‑1 (MQ=255)
                                              cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca                   <  1:2307079/71‑1 (MQ=255)
                                              cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca                   <  1:1215325/71‑1 (MQ=255)
                                              cTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc                    <  1:1925143/70‑1 (MQ=255)
                                               tAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                              >  1:3445288/1‑59 (MQ=255)
                                               tAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                              >  1:3941392/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)
                                                           ttGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCa                                      >  1:813325/1‑39 (MQ=255)
                                                            tGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt                     <  1:3012577/55‑1 (MQ=255)
                                                            tGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt                     <  1:2280961/55‑1 (MQ=255)
                                                               aGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAt  <  1:5002906/71‑1 (MQ=255)
                                                               aGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAt  <  1:87188/71‑1 (MQ=255)
                                                                gAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTcc          >  1:2853965/1‑62 (MQ=255)
                                                                      ttGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCg         <  1:178184/57‑1 (MQ=35)
                                                                      ttGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCg         <  1:4474811/57‑1 (MQ=35)
                                                                      |                                                               
TATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAT  >  Exported/1591590‑1591723

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: