Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,540	A→G	6.7%	noncoding (304/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,763,540	0	A	G	6.7%	198.6 / 6.1	75	noncoding (304/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (25/45);  new base G (2/3);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT  >  Exported/1763501‑1763608                                        |                                                                     aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                       >  1:2080024/1‑71 (MQ=14) aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                       >  1:449594/1‑71 (MQ=14) aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                       >  1:423641/1‑71 (MQ=14) aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                       >  1:3045702/1‑71 (MQ=14) aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                       >  1:2701320/1‑71 (MQ=14)  ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                      >  1:31987/1‑71 (MQ=32)  ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                      <  1:2274331/71‑1 (MQ=17)   gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTgg                                     <  1:490577/71‑1 (MQ=11)   gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTgg                                     <  1:2481990/71‑1 (MQ=11)    aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                    <  1:1269786/71‑1 (MQ=34)     gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTg                                   <  1:566110/71‑1 (MQ=11)     gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                    >  1:3175186/1‑70 (MQ=25)      tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                  <  1:2400803/71‑1 (MQ=21)      tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                  <  1:187055/71‑1 (MQ=21)        tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                     >  1:2321549/1‑66 (MQ=25)        tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                     >  1:407457/1‑66 (MQ=25)        tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                >  1:468084/1‑71 (MQ=21)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                              >  1:2667952/1‑71 (MQ=34)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                              >  1:1050636/1‑71 (MQ=34)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                               <  1:3385288/70‑1 (MQ=25)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                               <  1:2383514/70‑1 (MQ=25)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                <  1:896841/69‑1 (MQ=25)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                <  1:2823864/69‑1 (MQ=25)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                <  1:2316001/69‑1 (MQ=25)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                              >  1:2768645/1‑71 (MQ=21)          gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                              >  1:268553/1‑71 (MQ=21)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:1759088/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:158616/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:1755242/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:593331/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:1094390/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:2343415/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:2628856/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:1250629/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:3369470/71‑2 (MQ=11)          gACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGcca                              <  1:344131/71‑2 (MQ=11)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                   >  1:1817984/1‑65 (MQ=35)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                   >  1:2592988/1‑65 (MQ=35)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                   >  1:3395954/1‑65 (MQ=35)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTg                                   >  1:1622022/1‑65 (MQ=35)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                             <  1:2657554/71‑1 (MQ=21)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                             >  1:2791192/1‑71 (MQ=21)           aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                             <  1:2984759/71‑1 (MQ=21)            ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                            <  1:2269088/71‑1 (MQ=34)            ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                            <  1:1534099/71‑1 (MQ=34)             cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                            <  1:1459344/70‑1 (MQ=35)             cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                            >  1:2770499/1‑70 (MQ=35)             cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                           >  1:1172289/1‑71 (MQ=34)                 tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcc                               >  1:3054889/1‑63 (MQ=11)                 tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                        >  1:3082721/1‑70 (MQ=25)                   tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                     >  1:2117497/1‑71 (MQ=34)                     aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                            <  1:2906679/62‑1 (MQ=25)                     aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                           >  1:2429005/1‑63 (MQ=25)                         ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                <  1:1784568/70‑1 (MQ=255)                         ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                <  1:2341536/70‑1 (MQ=255)                         ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                <  1:2534359/70‑1 (MQ=255)                         ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                 >  1:143677/1‑69 (MQ=40)                             tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa           <  1:2716738/71‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa           <  1:123096/71‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa           <  1:203697/71‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa           <  1:3391515/71‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACAAGCTa           <  1:1654433/71‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                      <  1:1004419/44‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                      <  1:1141152/44‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                      <  1:2172058/44‑1 (MQ=255)                             tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                      <  1:2456802/44‑1 (MQ=255)                               tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa          <  1:2333982/70‑1 (MQ=255)                               tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt         <  1:2620077/71‑1 (MQ=255)                                 gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                          <  1:3071650/52‑1 (MQ=255)                                 gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                          <  1:1472666/52‑1 (MQ=255)                                 gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc        <  1:736120/70‑1 (MQ=255)                                 gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc        <  1:2296588/70‑1 (MQ=255)                                   tGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt         >  1:42192/1‑67 (MQ=11)                                     atCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATccga     <  1:1856091/68‑3 (MQ=11)                                      tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt  <  1:2828359/71‑1 (MQ=255)                                        |                                                                     AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCT  >  Exported/1763501‑1763608

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