Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 880,753 G→A 66.0% intergenic (+72/‑122) yciM → / → pyrF conserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE880,7530GA66.0% 51.7 / 47.5 50intergenic (+72/‑122)yciM/pyrFconserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (4/13);  new base A (10/23);  total (14/36)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCT  >  minE/880693‑880819
                                                            |                                                                  
aaTCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaaa                                                          <  1:2018017/71‑1 (MQ=255)
  tCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGc                                                        <  1:1658539/71‑1 (MQ=255)
        tttAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaaa                                                          <  1:2526127/63‑1 (MQ=255)
        tttAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    >  1:4061281/1‑69 (MQ=255)
            gTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAggg                                               <  1:3638266/70‑1 (MQ=255)
                caacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTg                                                     <  1:1353119/60‑1 (MQ=255)
                caacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGaa                                          <  1:3662649/71‑1 (MQ=255)
                   caTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAAt                                        >  1:1435265/1‑70 (MQ=255)
                    aTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCg                                      <  1:2001735/71‑1 (MQ=255)
                     tACTAATTATTAATGTTCCATTGTGCTTCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                     <  1:628437/71‑1 (MQ=255)
                     tACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                     <  1:1286730/71‑1 (MQ=255)
                       cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCa                                    <  1:1693031/70‑1 (MQ=255)
                         aaTTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    <  1:737123/52‑1 (MQ=255)
                         aaTTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    <  1:2940022/52‑1 (MQ=255)
                         aaTTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    <  1:3780275/52‑1 (MQ=255)
                         aaTTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    <  1:2019831/52‑1 (MQ=255)
                         aaTTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    <  1:601841/52‑1 (MQ=255)
                         aaTTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGc                                                    <  1:2727078/52‑1 (MQ=255)
                          attattAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAggg                                               <  1:1151154/56‑1 (MQ=255)
                             attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga                                           <  1:1145647/57‑1 (MQ=255)
                             attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                     <  1:4103888/63‑1 (MQ=255)
                             attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGtt                              >  1:3765661/1‑70 (MQ=255)
                              ttaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                     <  1:814687/62‑1 (MQ=255)
                                 aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAttt                         <  1:383897/71‑1 (MQ=255)
                                    ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAttttt                       <  1:2371494/70‑1 (MQ=255)
                                    ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAttttt                       <  1:505225/70‑1 (MQ=255)
                                    ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAttttt                       >  1:17839/1‑70 (MQ=255)
                                    ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAttttt                       <  1:2928400/70‑1 (MQ=255)
                                      ccATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTa                     <  1:2505342/70‑1 (MQ=255)
                                       cATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTAtt                   >  1:3221164/1‑71 (MQ=255)
                                         ttGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgag                                            <  1:1447155/44‑1 (MQ=255)
                                         ttGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgag                                            <  1:3073212/44‑1 (MQ=255)
                                          tgtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTAttt                  >  1:291834/1‑69 (MQ=255)
                                          tgtgCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTc                 >  1:873255/1‑70 (MQ=255)
                                          tgtgCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTc                 >  1:3199666/1‑70 (MQ=255)
                                           gtgCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTc                 >  1:1302132/1‑69 (MQ=255)
                                            tgCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc              <  1:1040397/71‑1 (MQ=255)
                                                ccGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACggg           >  1:2849932/1‑70 (MQ=255)
                                                 cGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACggg           >  1:465657/1‑69 (MQ=255)
                                                 cGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACg             <  1:2663887/67‑1 (MQ=255)
                                                 cGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACg             <  1:2221353/67‑1 (MQ=255)
                                                 cGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTa         <  1:3949606/71‑1 (MQ=255)
                                                 cGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTa         <  1:2284736/71‑1 (MQ=255)
                                                   gCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACggg           >  1:3632574/1‑67 (MQ=255)
                                                   gCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAg        >  1:2712670/1‑70 (MQ=255)
                                                     aaCGCCGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAAt     >  1:2196243/1‑71 (MQ=255)
                                                     aaCGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAAt     <  1:3049270/71‑1 (MQ=255)
                                                     aaCGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAAt     <  1:2100411/71‑1 (MQ=255)
                                                        gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCt  <  1:3801442/71‑1 (MQ=255)
                                                        gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCt  <  1:1237018/71‑1 (MQ=255)
                                                            |                                                                  
AATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCT  >  minE/880693‑880819

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: