Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	406,585	A→G	12.3%	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	12.3%	226.1 / 20.4	81	D8D (GAT→GAC)	ybeB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/44);  new base G (2/8);  total (29/52)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TGGATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGT  >  Exported/406515‑406649                                                                       |                                                                 tGGATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACaaaa                                                                  >  1:4635296/1‑71 (MQ=255)   gATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAAt                                                                 <  1:4569195/70‑1 (MQ=255)   gATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAAt                                                                 <  1:1402359/70‑1 (MQ=255)   gATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAAt                                                                 <  1:1777813/70‑1 (MQ=255)    aTTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATcc                                                               <  1:1635496/71‑1 (MQ=255)    aTTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATcc                                                               <  1:1640913/71‑1 (MQ=255)     tttGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCt                                                              >  1:1003388/1‑71 (MQ=255)     tttGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCt                                                              >  1:2179302/1‑71 (MQ=255)        gCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTgg                                                            <  1:3741250/70‑1 (MQ=255)         cccTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAg                                                          <  1:4684378/71‑1 (MQ=255)         cccTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAg                                                          <  1:1362087/71‑1 (MQ=255)         cccTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAg                                                          <  1:2532706/71‑1 (MQ=255)           cTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTg                                                        <  1:3676372/71‑1 (MQ=255)           cTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTg                                                        >  1:3303855/1‑71 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCGTCAATTTTGTCGATAACAAAGTCCTGGAGTGCtt                                                     <  1:4221199/71‑1 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTg                                                             <  1:2981549/63‑1 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTg                                                             <  1:1569433/63‑1 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCtt                                                     <  1:4702282/71‑1 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCtt                                                     <  1:1218317/71‑1 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCtt                                                     <  1:1267448/71‑1 (MQ=255)              aaCGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGc                                                       <  1:828900/69‑1 (MQ=255)                cGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAAt                                                                 >  1:189966/1‑57 (MQ=255)                cGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAAt                                                                 >  1:807541/1‑57 (MQ=255)                cGTCTAAGGCGATGATGTCCTGACCTTTAAGGTCATCAATTTTGTCGATAACAAAAt                                                                 >  1:1794421/1‑57 (MQ=255)                  tCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCt                                                              <  1:2838573/58‑1 (MQ=255)                    tAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGa                                                           >  1:3031482/1‑59 (MQ=255)                      aGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGc                                             <  1:3510138/71‑1 (MQ=255)                      aGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGc                                             <  1:2611574/71‑1 (MQ=255)                      aGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGc                                             <  1:2168353/71‑1 (MQ=255)                      aGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGc                                             <  1:2858627/71‑1 (MQ=255)                       ggCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCa                                            <  1:78348/71‑1 (MQ=255)                        gCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTa                                                   >  1:4714632/1‑63 (MQ=255)                        gCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTa                                                   >  1:4670722/1‑63 (MQ=255)                         cGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAg                                          <  1:4133886/71‑1 (MQ=255)                             gatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAcc                                                 >  1:4190851/1‑60 (MQ=255)                             gatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAcc                                                 >  1:4736644/1‑60 (MQ=255)                             gatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGtttt                                      >  1:2251629/1‑71 (MQ=255)                             gatgTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGtttt                                      >  1:3297455/1‑71 (MQ=255)                                gTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTccc                                   >  1:3942917/1‑71 (MQ=255)                                gTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTccc                                   >  1:1679103/1‑71 (MQ=255)                                  ccTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATcc                                                               >  1:1950462/1‑41 (MQ=255)                                   cTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCtt                                                     <  1:1377126/50‑1 (MQ=255)                                   cTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCt                                <  1:3280898/71‑1 (MQ=255)                                   cTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCt                                <  1:1463963/71‑1 (MQ=255)                                     gcccTTTGAGGTGATCAATTTTTTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCAgg                              >  1:492041/3‑71 (MQ=255)                                     gACCTTTGAGGTCATCAATTTTGTCGATAACAAAATACTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTgg                              >  1:588068/1‑71 (MQ=255)                                      aCCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTggg                             >  1:2499772/1‑71 (MQ=255)                                         tttGAGGTCATCAATTTTGTCGATAACAAAATCCt                                                              <  1:3272249/35‑1 (MQ=255)                                         tttGAGGTCATCAATTTTGTCGATAACAAAATCCt                                                              <  1:498776/35‑1 (MQ=255)                                         tttGAGGTCATCAATTTTGTCGATAACAAAATCCt                                                              <  1:793369/35‑1 (MQ=255)                                          ttGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGc                                             <  1:1643249/51‑1 (MQ=255)                                           tGAGGTCATCAATTTTGTCGATAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTggg                             >  1:1489510/1‑66 (MQ=255)                                            gAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGa                          <  1:3142444/68‑1 (MQ=255)                                             aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCa                      <  1:4283915/71‑1 (MQ=255)                                             aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCa                      <  1:4525109/71‑1 (MQ=255)                                             aGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCa                      <  1:2097213/71‑1 (MQ=255)                                                 catcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGt                            >  1:4582131/1‑61 (MQ=255)                                                 catcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGt                            >  1:2450093/1‑61 (MQ=255)                                                 catcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaaa                    >  1:2179370/1‑69 (MQ=255)                                                  atcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGt                            <  1:3361062/60‑1 (MQ=255)                                                  atcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGt                            <  1:3959961/60‑1 (MQ=255)                                                  atcaATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGt                            <  1:2964374/60‑1 (MQ=255)                                                    caATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGt                            <  1:1817425/58‑1 (MQ=255)                                                    caATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGa                          >  1:4405579/1‑60 (MQ=255)                                                    caATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGa                >  1:4113722/1‑70 (MQ=255)                                                     aaTTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGAt               >  1:3595887/1‑70 (MQ=255)                                                     aaTTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGAt               >  1:3544434/1‑70 (MQ=255)                                                     aaTTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGAt               >  1:3382374/1‑70 (MQ=255)                                                       ttttGTCGATAACAAAATCCTGGAGTGCTTTACCCt                                               <  1:3389442/36‑1 (MQ=255)                                                       ttttGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATaaa            <  1:3929604/71‑1 (MQ=255)                                                            tCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATg        <  1:1257417/70‑1 (MQ=255)                                                             cGATAACAAAATCCTGGAGTGCTTTACCCTGCAAg                                          <  1:1889415/35‑1 (MQ=255)                                                              taTAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTg                               <  1:1184553/44‑1 (MQ=255)                                                              gATAACAAAATCCTGGAGTGCTTTACCCTGCAAGtttt                                      <  1:559033/38‑1 (MQ=255)                                                              gATAACAAAATCCTGGAGTGCTTTACCCTGCAAGtttt                                      <  1:2326202/38‑1 (MQ=255)                                                              gATAACAAAATCCTGGAGTGCTTTACCCTGCAAGtttt                                      <  1:932513/38‑1 (MQ=255)                                                              gATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATg        <  1:916325/68‑1 (MQ=255)                                                              gATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGt      <  1:2174579/70‑1 (MQ=255)                                                                tAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg   <  1:4476063/71‑1 (MQ=255)                                                                tAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg   >  1:4198354/1‑71 (MQ=255)                                                                tAACAAAGTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg   <  1:2874483/71‑1 (MQ=255)                                                                  aCAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGt  <  1:4251557/70‑1 (MQ=255)                                                                     aaaTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGaa                         >  1:3967522/1‑44 (MQ=255)                                                                     aaaTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGaa                         >  1:3767806/1‑44 (MQ=255)                                                                     aaaTCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGaa                         >  1:3001388/1‑44 (MQ=255)                                                                       |                                                                 TGGATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGT  >  Exported/406515‑406649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑