Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	980,321	G→A	G412D (GGC→GAC)	rpsA →	30S ribosomal protein S1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	980,321	0	G	A	96.2%	195.2 / ‑3.8	53	G412D (GGC→GAC)	rpsA	30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/25);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAGCTCGC  >  REL606/980254‑980396                                                                    |                                                                            cTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGaaa                                                                     <  2:1392909/76‑1 (MQ=255)   ggTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATc                                                                   >  1:1180290/1‑76 (MQ=255)   ggTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATc                                                                   <  2:344313/76‑1 (MQ=255)   ggTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATc                                                                   >  1:889742/1‑76 (MQ=255)   ggTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATc                                                                   >  1:774164/1‑76 (MQ=255)   ggTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATc                                                                   >  2:844235/1‑76 (MQ=255)       cACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTg                                                               <  2:1274016/76‑1 (MQ=255)           tGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAgt                                                           <  1:507039/76‑1 (MQ=255)               tGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAgttgtt                                                       <  2:1432963/76‑1 (MQ=255)                gACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTc                                                      >  2:943702/1‑76 (MQ=255)                gACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTc                                                      >  2:88016/1‑76 (MQ=255)                 aCATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCt                                                     <  2:22634/76‑1 (MQ=255)                 aCATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCt                                                     <  1:1407586/76‑1 (MQ=255)                 aCATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCt                                                     >  1:517327/1‑76 (MQ=255)                     ctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAg                                                 <  2:689711/76‑1 (MQ=255)                      tcCTGGAATGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAgg                                                >  1:657317/1‑76 (MQ=255)                        cTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGtt                                              >  1:74906/1‑76 (MQ=255)                        cTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGtt                                              <  1:724607/76‑1 (MQ=255)                         tGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTg                                             >  2:1558783/1‑76 (MQ=255)                         gggAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTg                                             <  1:924408/75‑1 (MQ=255)                          ggAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGa                                            >  1:410911/1‑76 (MQ=255)                          ggAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGa                                            >  2:349225/1‑76 (MQ=255)                          ggAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGa                                            >  2:820684/1‑76 (MQ=255)                            aaCGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACg                                          <  2:1018475/76‑1 (MQ=255)                               gTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAg                                       <  1:1261915/76‑1 (MQ=255)                               gTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAg                                       <  1:267217/76‑1 (MQ=255)                               gTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAg                                       >  1:1348378/1‑76 (MQ=255)                                ttGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAga                                      >  1:140580/1‑76 (MQ=255)                                  gCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaac                                    <  1:566939/76‑1 (MQ=255)                                  gCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaac                                    <  1:1316104/76‑1 (MQ=255)                                        gaagaaGCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaac                              <  1:930261/76‑1 (MQ=255)                                          agaagCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAgaacgtgaacgt                            <  2:1343371/76‑1 (MQ=255)                                            aagCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTAt                          <  2:867798/76‑1 (MQ=255)                                              gCAGTTCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTAtct                        <  2:505577/76‑1 (MQ=255)                                                   tCGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTg                   <  1:1445906/76‑1 (MQ=255)                                                    cGTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAAAGTATCTTCCttg                  >  2:952744/1‑74 (MQ=255)                                                     gTGAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTggg                 >  2:1218431/1‑76 (MQ=255)                                                       gAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCg               >  2:774310/1‑76 (MQ=255)                                                       gAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCg               <  2:210091/76‑1 (MQ=255)                                                       gAATACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCACCCTGGGCg               >  1:1170524/1‑76 (MQ=255)                                                        aaTACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGt              >  2:770060/1‑76 (MQ=255)                                                         aTACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGtt             <  1:192503/76‑1 (MQ=255)                                                          tACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTa            <  2:263256/76‑1 (MQ=255)                                                          tACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGGTa            >  1:137642/1‑76 (MQ=255)                                                          tACAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGGATCTCCCCGGGCGttt            >  1:893206/1‑75 (MQ=255)                                                           aCAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTaa           <  2:1295023/76‑1 (MQ=255)                                                           aCAAAAAAGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTaa           >  2:1233393/1‑76 (MQ=255)                                                             aaaaaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAAc         >  2:1539089/1‑76 (MQ=255)                                                              aaaaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACa        >  1:1570016/1‑76 (MQ=255)                                                               aaaaGGCGGGTAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAg       <  2:307584/76‑1 (MQ=255)                                                               aaaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAg       >  2:1239946/1‑76 (MQ=255)                                                               aaaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGATAAACAg       <  1:1024976/76‑1 (MQ=255)                                                                aaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAGc      >  2:1173722/1‑76 (MQ=255)                                                                aaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAGc      >  1:65102/1‑76 (MQ=255)                                                                aaaGACGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAGc      <  1:400176/76‑1 (MQ=255)                                                                    gCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAGCTCGc  >  2:1404505/1‑76 (MQ=255)                                                                    |                                                                            CTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAATCGCTGCAGTTGTTCTGCAGGTTGACGCAGAACGTGAACGTATCTCCCTGGGCGTTAAACAGCTCGC  >  REL606/980254‑980396

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑