Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,352,815	G→A	A92V (GCA→GTA)	rpsI ←	ribosomal protein S9

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,352,815	0	G	A	94.4%	33.9 / ‑3.2	19	A92V (GCA→GTA)	rpsI	ribosomal protein S9
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (10/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

CCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTC  >  CP000730/2352684‑2352945                                                                                                                                    |                                                                                                                                   ccGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCtaata                                                                                                                           <  1:2723/141‑1 (MQ=255)    aCAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTGCTAATACTg                                                                                                                        >  2:345305/1‑141 (MQ=255)     gagttAGCGTTTTGAGAATGGGGGGGAACGACGAGCTGCGTTAAGAGCTGGTTTTTTGCGTTTGTTCAGACGTGGGTGACGAGGAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGc                                                                                                                       <  2:425559/138‑1 (MQ=255)     cAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTCCTTTCATACGTGGGTCACGAGTAAGTACTCCAGCGCGGTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAGTAGTGc                                                                                                                       >  2:71176/1‑141 (MQ=38)               tttGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTCCGTTCTTTCCTACCTGGGTCCCCACTAAGTCATCCAGCGCGTTTTCCAGAACCTCTGTATTCAGGCTCTACTTCTAATCATGCACGAGCGAtt                                                                                                             >  2:177/1‑141 (MQ=37)                          aGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTTTTTTTTTTTTTACGTAGGACACGCGTAAGTAATCCAACGCTGGTTAAAGAATCTCTATCTTCTGGATCTCCTCCCAATAATGCACGAGCGATTCCGTGACggct                                                                                                  >  2:235766/1‑139 (MQ=255)                          aGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAAGATCTACTTCTAATAATGCACGAGCGATTCCATGagcgat                                                                                                  >  2:476664/1‑136 (MQ=255)                                                     tGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccacc                                                                       <  2:408338/141‑1 (MQ=255)                                                       gatctttACGTGCCTGCATACGCGGGTCACGAGCAAGTAATCCAGCGCGTTTTAAAGAACAGCTGGATTCAGGATCTACTGCTAATAAGGCACGAGCCAGCCCGTGACGGATAGCTTGAGCTTGACCAGCGGAacaacccc                                                                     >  2:231912/5‑139 (MQ=255)                                                        ttttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccaccccc                                                                    >  2:461833/1‑141 (MQ=255)                                                         tttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAAccaccacca                                                                   >  2:417347/1‑141 (MQ=255)                                                                            gTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAAcat                                                <  2:213841/141‑1 (MQ=255)                                                                                       gtaagtaaTCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTCCCt                                     >  2:37299/1‑141 (MQ=255)                                                                                        taagtaaTCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCtt                                    <  2:423654/141‑1 (MQ=255)                                                                                                    cgcgTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGCAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTa                        >  2:469581/1‑141 (MQ=255)                                                                                                       gTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCGGTTACAt                     >  2:153260/1‑141 (MQ=255)                                                                                                         tgtAAAGATCCTGTGTATGCAGGATCTACTTTGTAGTATGCACGAGCGATTCCGTGACGGATAGTTTGAGCTTGTCCCGTGAGACCACCACCATGAACGTTAGCTAAAACATCATAGTTAGCTTTAGTTTCAGTTACATCa                   <  2:331967/139‑1 (MQ=255)                                                                                                           tcactacccTCTGTATTCAGGATCTACTTGTAATAATGCACGATCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTATCTTTAGTTTCAGTTACATCaaa                 <  2:275104/134‑1 (MQ=255)                                                                                                                          ttCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTc  <  2:455090/141‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   CCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTC  >  CP000730/2352684‑2352945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑