Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,359,867	G→A	29.5%	noncoding (1228/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,359,867	0	G	A	29.5%	91.1 / 48.7	61	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (33/10);  new base A (9/9);  total (42/19)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  Exported/2359818‑2359936                                                  |                                                                      gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                  >  1:24964/1‑71 (MQ=34) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                  >  1:1412868/1‑71 (MQ=255) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                  <  1:66456/71‑1 (MQ=255)  cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                 >  1:469020/1‑71 (MQ=255)    ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                               >  1:1976749/1‑71 (MQ=255)     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                               <  1:91623/70‑1 (MQ=35)     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                <  1:1054478/69‑1 (MQ=255)      aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             >  1:1897381/1‑71 (MQ=34)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             <  1:1976499/68‑1 (MQ=35)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                          <  1:458686/71‑1 (MQ=255)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                         <  1:450087/71‑1 (MQ=255)           tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             >  1:502575/1‑66 (MQ=35)           tACATATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                     >  1:148121/1‑58 (MQ=9)            aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                        >  1:742612/1‑70 (MQ=255)              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:878157/1‑68 (MQ=35)              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:963883/1‑68 (MQ=35)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                            >  1:995360/1‑62 (MQ=35)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                     <  1:95/69‑1 (MQ=255)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                    >  1:1343740/1‑70 (MQ=255)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                   >  1:617124/1‑71 (MQ=255)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                    <  1:383399/64‑1 (MQ=11)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                             <  1:1631892/71‑1 (MQ=255)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                >  1:1598157/1‑67 (MQ=35)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                            <  1:117699/71‑1 (MQ=34)                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                           >  1:137019/1‑70 (MQ=35)                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                     <  1:1477282/60‑1 (MQ=11)                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                           >  1:300083/1‑70 (MQ=255)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                     >  1:1403569/1‑42 (MQ=32)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                 >  1:93021/1‑62 (MQ=35)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:844030/1‑70 (MQ=255)                            gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                               <  1:779038/47‑1 (MQ=38)                             cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                      <  1:1940499/71‑1 (MQ=34)                              aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                   >  1:1672312/1‑57 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                    >  1:1596312/1‑71 (MQ=35)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:1823441/1‑65 (MQ=255)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                    >  1:494932/1‑37 (MQ=37)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                    >  1:161980/1‑37 (MQ=37)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                    >  1:444005/1‑37 (MQ=37)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                     >  1:807939/1‑36 (MQ=33)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                  >  1:553395/1‑71 (MQ=39)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                  >  1:1729478/1‑38 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                  >  1:1679362/1‑70 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                  >  1:1799976/1‑70 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                  <  1:1574592/70‑1 (MQ=40)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                 <  1:1512889/68‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt              >  1:366999/1‑71 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt              <  1:1340718/71‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg               >  1:353240/1‑70 (MQ=255)                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt            >  1:343933/1‑71 (MQ=255)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc         >  1:1302178/1‑71 (MQ=255)                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        >  1:392370/1‑71 (MQ=255)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:524157/1‑53 (MQ=35)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                    >  1:520875/1‑58 (MQ=38)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg       >  1:1461840/1‑71 (MQ=255)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg       >  1:1198048/1‑71 (MQ=255)                                            aaGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         >  1:514411/2‑37 (MQ=255)                                              gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                      <  1:2050315/38‑1 (MQ=40)                                              gCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                 <  1:80719/43‑1 (MQ=14)                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:879528/1‑35 (MQ=40)                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg    >  1:1261977/1‑71 (MQ=255)                                                 tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCa  <  1:1081925/71‑1 (MQ=255)                                                  |                                                                      GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  Exported/2359818‑2359936

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