Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	406,585	A→G	46.7%	D8D (GAT→GAC)	ybeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	406,585	0	A	G	46.7%	1.7 / 61.1	45	D8D (GAT→GAC)	ybeB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/8);  new base G (9/12);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

GATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTC  >  Exported/406517‑406650                                                                     |                                                                  gATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATc                                                                 <  1:1884476/71‑1 (MQ=255)      gCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTg                                                              >  1:1800802/1‑69 (MQ=255)         cTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTg                                                         >  1:2232842/1‑71 (MQ=255)           gAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCt                                                       >  1:313695/1‑71 (MQ=255)            aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCtt                                                      >  1:2966350/1‑71 (MQ=255)            aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTg                                                         <  1:1704100/68‑1 (MQ=255)              cGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTa                                                    <  1:268840/71‑1 (MQ=255)                 cTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTAccc                                                 >  1:1392160/1‑71 (MQ=255)                 cTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAccc                                                 >  1:3244107/1‑71 (MQ=255)                     ggCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCa                                             <  1:2360939/71‑1 (MQ=255)                      gCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGc                                              >  1:2864941/1‑69 (MQ=255)                           gatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGttt                                        >  1:1148960/1‑70 (MQ=255)                           gatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGttt                                        >  1:1492080/1‑70 (MQ=255)                            atgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGtttt                                       <  1:2911575/70‑1 (MQ=255)                              gTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTccc                                    >  1:1579257/1‑71 (MQ=255)                                ccTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTg                                               >  1:2135309/1‑58 (MQ=255)                                 cTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTccccc                                  <  1:2282403/70‑1 (MQ=255)                                 cTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGc                                                        <  1:1630615/48‑1 (MQ=255)                                  tGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTg                                >  1:3529373/1‑71 (MQ=255)                                      cTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTg                            <  1:2494046/71‑1 (MQ=255)                                        ttGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGaa                          >  1:3484746/1‑71 (MQ=255)                                         tGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAAt                         >  1:3050504/1‑71 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCa                       <  1:3079829/71‑1 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTa                                                    >  1:2853465/1‑42 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAc                                                   >  1:2037456/1‑43 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAc                                                   >  1:2723711/1‑43 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAc                                                   >  1:2941475/1‑43 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAc                                                   >  1:694412/1‑43 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAc                                                   >  1:3200224/1‑43 (MQ=255)                                           aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAc                                                   >  1:3201290/1‑43 (MQ=255)                                            ggTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaa                      <  1:3627398/71‑1 (MQ=255)                                             gTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaaa                     <  1:3630924/71‑1 (MQ=255)                                              tcatcaATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCa                       <  1:190850/68‑1 (MQ=255)                                               catcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaaa                     >  1:2769293/1‑69 (MQ=255)                                                   aaTTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGtttt                                       <  1:3194572/47‑1 (MQ=255)                                                     ttttGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATaaa             <  1:852277/71‑1 (MQ=255)                                                     ttttGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATaa              >  1:1413536/1‑70 (MQ=255)                                                       ttGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGaa                          <  1:461119/56‑1 (MQ=255)                                                        tGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATaaa             <  1:1790260/68‑1 (MQ=255)                                                         gTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAAt          <  1:3378555/70‑1 (MQ=255)                                                          tCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATgg        >  1:3460963/1‑71 (MQ=255)                                                            gATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTccc                                    <  1:2242440/41‑1 (MQ=255)                                                            gATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGt       >  1:2711044/1‑70 (MQ=255)                                                             aTAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCt     <  1:104887/71‑1 (MQ=255)                                                                aCAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTc  <  1:770589/71‑1 (MQ=255)                                                                     |                                                                  GATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTC  >  Exported/406517‑406650

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: ‑