Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,352,815	G→A	85.1%	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	85.1%	68.1 / 5.5	33	A92V (GCA→GTA)	rpsI	ribosomal protein S9
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (14/14);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAA  >  CP000730/2352688‑2352948                                                                                                                                |                                                                                                                                      cAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATga                                                                                                                          >  2:187318/1‑140 (MQ=255)        cGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCg                                                                                                                   >  1:326345/1‑141 (MQ=255)        cGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCg                                                                                                                   >  1:172950/1‑141 (MQ=255)                     gAGGTGAGCGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGa                                                                                                      <  1:301034/141‑1 (MQ=255)                     gAGGTGACCGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGa                                                                                                      <  1:301038/141‑1 (MQ=255)                           aaCCACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGCGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATNCCGTGACGGATagctt                                                                                                <  1:210168/141‑1 (MQ=255)                               acgaGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATagcttgagc                                                                                            <  2:256956/141‑1 (MQ=255)                               acgaGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATagcttgagc                                                                                            <  2:204342/141‑1 (MQ=255)                               acgaGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATagcttgagc                                                                                            <  2:305300/141‑1 (MQ=255)                                   gctgctTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGt                                                                                        <  1:16968/141‑1 (MQ=255)                                    ctgctTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTc                                                                                       <  2:87181/141‑1 (MQ=255)                                             gACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGCGCTTGTCCAGTGAAac                                                                              >  1:84379/1‑141 (MQ=255)                                              aCCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAacc                                                                             >  1:169625/1‑141 (MQ=255)                                                 tGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGACaccacc                                                                          >  2:244554/1‑141 (MQ=255)                                                 tGGTTTGTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccacc                                                                          <  1:157149/141‑1 (MQ=255)                                                   gTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccaccac                                                                        >  1:234005/1‑141 (MQ=255)                                                                tttCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAAc                                                           <  2:48541/141‑1 (MQ=255)                                                                             tCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAg                                              <  2:212864/141‑1 (MQ=255)                                                                                        taaTCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGt                                   >  1:263146/1‑141 (MQ=255)                                                                                          aTCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTACAACATCATAGTTACCTTTAGttt                                 >  2:191642/1‑141 (MQ=255)                                                                                               gcgcGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGtt                            >  1:288379/1‑141 (MQ=255)                                                                                                cgcgTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTa                           >  2:20193/1‑141 (MQ=255)                                                                                                 gcgTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTAc                          >  2:156674/1‑141 (MQ=255)                                                                                                   gTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCTCCACCATGACCGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACAt                        <  2:108167/141‑1 (MQ=255)                                                                                                   gTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACAt                        >  2:305202/1‑141 (MQ=255)                                                                                                   gTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACAt                        >  1:27741/1‑141 (MQ=255)                                                                                                     tttAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCCCCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCa                      <  1:256107/141‑1 (MQ=255)                                                                                                          aGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAAtgg                 <  1:134302/141‑1 (MQ=255)                                                                                                          aGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAAtgg                 <  1:305039/141‑1 (MQ=255)                                                                                                              ccTCTGTATGCAGGAGCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTCCAGTTACaaaaaagaacaac             >  2:12057/1‑129 (MQ=255)                                                                                                                tctGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATggttggtt           <  1:30490/141‑1 (MQ=255)                                                                                                                tctGTATTCAGGATCTACTTCTAATAATGAACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTTAAACCACCACCATGAACGTTAACTAAAACATAATAGTTACCTTTAGTTTCAGTTACATCAAATggttggtt           <  2:69214/141‑1 (MQ=255)                                                                                                                         aGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTaa  >  1:273765/1‑141 (MQ=255)                                                                                                                                |                                                                                                                                      CAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAA  >  CP000730/2352688‑2352948

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