Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,591,605 G→A 33.1% intergenic (+64/+125) 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,6050GA33.1% 158.6 / 120.8 145intergenic (+64/+125)yfeS/cysMconserved hypothetical protein/cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base G (43/54);  new base A (29/19);  total (72/73)
Fisher's exact test for biased strand distribution p-value = 7.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

AAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAG  >  Exported/1591535‑1591675
                                                                      |                                                                      
aaaTTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                        >  1:2273374/1‑71 (MQ=255)
aaaTTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                        >  1:1867160/1‑71 (MQ=255)
   ttAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa                                                                     <  1:2699824/71‑1 (MQ=255)
       aTACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTTTTTATGCCGGATGCGGCGTAAACGc                                                                 <  1:1482433/71‑1 (MQ=255)
            gCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                    >  1:644159/1‑63 (MQ=255)
             gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATaaa                                                                    <  1:3281216‑M1/15‑1 (MQ=255)
             gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCtt                                                              >  1:257903‑M1/48‑68 (MQ=255)
             gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTa                                                             >  1:284181‑M1/48‑69 (MQ=255)
               tttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                      <  1:5493286/58‑1 (MQ=255)
               ataggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGc                                                                 >  1:1068605‑M1/46‑63 (MQ=255)
                 tGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                    >  1:5612222/1‑58 (MQ=255)
                 aggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAAc                                                                   <  1:1522535‑M1/16‑1 (MQ=255)
                 aggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAAc                                                                   <  1:5605964‑M1/16‑1 (MQ=255)
                 aggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAAc                                                                   <  1:4449442‑M1/16‑1 (MQ=255)
                  ggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTAtccgacc                                                      <  1:5406054‑M1/29‑4 (MQ=255)
                  ggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTAtccgacc                                                      <  1:4881947‑M1/29‑4 (MQ=255)
                  ggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTAtccgacc                                                      <  1:3460856‑M1/29‑4 (MQ=255)
                  gTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                       <  1:2131209/70‑1 (MQ=255)
                    cGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                      <  1:2060160/53‑1 (MQ=255)
                     ggCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                      <  1:3732465/52‑1 (MQ=255)
                     ggCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                   <  1:4785641/71‑1 (MQ=255)
                     ctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTa                                                    <  1:41480‑M1/31‑1 (MQ=255)
                     ctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTa                                                    <  1:2746008‑M1/31‑1 (MQ=255)
                     ctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                   >  1:2185675‑M1/40‑71 (MQ=255)
                      tgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                   <  1:2792159‑M1/32‑1 (MQ=255)
                      tgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                   >  1:4427802‑M1/39‑70 (MQ=255)
                       gataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACa                                                  <  1:2192035‑M1/33‑1 (MQ=255)
                       gataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACa                                                  <  1:1148611‑M1/33‑1 (MQ=255)
                       gataagacgcggcaagcgtctcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCg                                                         >  1:146277‑M1/38‑63 (MQ=255)
                       gataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTa                                                                      >  1:3476384‑M1/38‑50 (MQ=255)
                       gataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACa                                                  <  1:653082‑M1/33‑1 (MQ=255)
                       ctgataagacggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACa                                                  <  1:1112508‑M1/33‑1 (MQ=255)
                        ataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaa                                                 <  1:4676571‑M1/34‑1 (MQ=255)
                        ataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACaa                                                 <  1:4663289‑M1/34‑1 (MQ=255)
                        ataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACa                                                  >  1:3968541‑M1/37‑69 (MQ=255)
                        ataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACa                                                  >  1:3133308‑M1/37‑69 (MQ=255)
                         taagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTAc                                                   >  1:2505612‑M1/36‑67 (MQ=255)
                         taagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTAc                                                   >  1:3689311‑M1/36‑67 (MQ=255)
                         taagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                               <  1:673634‑M1/36‑1 (MQ=255)
                          aagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                               <  1:3316933‑M1/36‑1 (MQ=255)
                          aagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAt                                              <  1:3698561‑M1/37‑1 (MQ=255)
                          aagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAt                                              <  1:1652032‑M1/37‑1 (MQ=255)
                          aagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                              <  1:5355782‑M1/37‑1 (MQ=255)
                          aagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                              <  1:4347264‑M1/37‑1 (MQ=255)
                           agacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATc                                             >  1:3931994‑M1/34‑71 (MQ=255)
                           agacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATc                                             >  1:455065‑M1/34‑71 (MQ=255)
                           agacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCg                                                         >  1:2936739‑M1/34‑59 (MQ=255)
                           aacaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                             >  1:4471248/1‑71 (MQ=255)
                           aacaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                             >  1:1993540/1‑71 (MQ=255)
                            gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAtca                                            >  1:1472609‑M1/33‑70 (MQ=255)
                            gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAtca                                            <  1:2919820‑M1/39‑2 (MQ=255)
                            gacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAtca                                            <  1:2232205‑M1/39‑2 (MQ=255)
                            acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                                             <  1:4865689/54‑1 (MQ=255)
                               acTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                                                      <  1:4041016/42‑1 (MQ=255)
                                cTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                    >  1:553422/1‑43 (MQ=255)
                                cTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACg                                                                  <  1:3071542/45‑1 (MQ=255)
                                  ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACaaaa                                               >  1:670278‑M1/27‑62 (MQ=255)
                                  ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa                                      >  1:4081563‑M1/27‑71 (MQ=255)
                                  gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                      >  1:3292707/1‑71 (MQ=255)
                                   ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                         >  1:235094/1‑67 (MQ=39)
                                   ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAt                                     >  1:4599416/1‑71 (MQ=255)
                                   caagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                     >  1:209120‑M1/26‑70 (MQ=255)
                                    tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAtt                                    >  1:2591310/1‑71 (MQ=40)
                                    tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAtt                                    >  1:2474244/1‑71 (MQ=40)
                                    aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                     >  1:2446982‑M1/25‑69 (MQ=255)
                                    aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                     >  1:2920001‑M1/25‑69 (MQ=255)
                                     agcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa                                     <  1:2228188‑M1/46‑2 (MQ=255)
                                     agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaatc                                   >  1:4238166‑M1/24‑68 (MQ=255)
                                     agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                    <  1:2082957‑M1/47‑3 (MQ=255)
                                      gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa                                  >  1:4811518‑M1/23‑71 (MQ=255)
                                      aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                   <  1:1177903/70‑1 (MQ=39)
                                      aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                  <  1:3609408/71‑1 (MQ=38)
                                      aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                  <  1:2427778/71‑1 (MQ=38)
                                       cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                 >  1:187267‑M1/22‑70 (MQ=255)
                                       cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                 >  1:3585021‑M1/22‑70 (MQ=255)
                                       cgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                 <  1:519117‑M1/50‑2 (MQ=255)
                                        gtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAAtcagg                                <  1:1299606‑M1/51‑3 (MQ=255)
                                         tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                    >  1:3384337‑M1/20‑64 (MQ=255)
                                         tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat                                    >  1:3391076‑M1/20‑64 (MQ=255)
                                         tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                               >  1:5118810‑M1/20‑68 (MQ=255)
                                         tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcagga                               >  1:1080715‑M1/20‑68 (MQ=255)
                                          cgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat                                    >  1:754138‑M1/19‑63 (MQ=255)
                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatata                             <  1:3988206/71‑1 (MQ=39)
                                             gcgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATAtt                           >  1:1995596/1‑71 (MQ=39)
                                                  gcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                 >  1:3826488‑M1/11‑59 (MQ=255)
                                                    gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                   <  1:2608126/56‑1 (MQ=35)
                                                       tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                   >  1:2781268/1‑69 (MQ=255)
                                                       tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatca                 >  1:305977/1‑68 (MQ=11)
                                                       tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatc                  <  1:4377490/70‑3 (MQ=11)
                                                       gatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                 <  1:4317859‑M1/50‑2 (MQ=255)
                                                       gatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag                                 <  1:228643‑M1/50‑2 (MQ=255)
                                                              gATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGc          >  1:2859712/1‑71 (MQ=11)
                                                              gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                        <  1:321210/41‑1 (MQ=32)
                                                              gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg           >  1:161687/1‑70 (MQ=255)
                                                              gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc          >  1:3033291/1‑71 (MQ=255)
                                                              gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc          <  1:487477/71‑1 (MQ=255)
                                                               aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                         >  1:1538705/1‑55 (MQ=255)
                                                                tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCt        <  1:4258487/71‑1 (MQ=255)
                                                                 gcggcgTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTa             >  1:4840718/1‑65 (MQ=11)
                                                                      gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg  >  1:1172367/1‑71 (MQ=255)
                                                                      |                                                                      
AAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAG  >  Exported/1591535‑1591675

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: