Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,943,314	C→A	14.9%	R489L (CGT→CTT)	aspS ←	aspartyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,943,314	0	C	A	14.9%	71.7 / 4.8	27	R489L (CGT→CTT)	aspS	aspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base A (2/2);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

CAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTC  >  CP009273/1943043‑1943547                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          tagcgcagcCAGTGCATGCGGGTTAGCAAAGCTCGGTGCTTCAGGCATCAGACACTCGGCCGGCGTGGTGTTCGGGAAGTCTAGAACGTCACGGATATTGGCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTGCGTTGATACTCAGAATACCAAACACCGTCTGCTGCATATCACCATTTTGGATACGTa                                                                                                                                                                                                                                        <  1:735230/274‑1 (MQ=255)          aggTGCAGTCGGGTTAGCCAAGCTCGGTGCTTCAGTCATCAGCCACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGCTATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAacca                                                                                                                                                                                                                                <  2:231660/272‑1 (MQ=255)                        tAGCAAAGCTCGGTGCTTCAGTCATCAGCCACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGGTATTGTCGGTGCCTGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCNGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTaa                                                                                                                                                                                                                 <  1:622384/275‑1 (MQ=255)                               gCTCGGTGCTTCATTCATCCGACACGCTTCCGCCGTGGTTTTCGGGAAGGCGATAACCTCACTGATAGTGTCTGTGCCGGTCNGCAGCATGNTCAGACGGTCAATACCGAATGCCAGACCTGCGTGCGTCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGCACCACCGCCCACTTCGTAACCATtga                                                                                                                                                                                                          <  1:475860/275‑1 (MQ=255)                                                                      tttCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCNGAGTACCGTATTTCAGAGCATCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTg                                                                                                                                                                   <  1:562281/275‑1 (MQ=255)                                                                                                  ttGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATAAGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCata                                                                                                                                                                                                  >  1:254986/1‑216 (MQ=255)                                                                                                      cGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATAAGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGt                                                                                                                                                                                              <  2:254983/216‑1 (MQ=255)                                                                                                                                                             cggcggAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAg                                                                                                                                                 <  1:709564/206‑1 (MQ=255)                                                                                                                                                                   aGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCa                                                                                                                                                                                                              >  1:282754/1‑139 (MQ=255)                                                                                                                                                                       ccGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATtga                                                                                                                                                                                                          <  2:282749/139‑1 (MQ=255)                                                                                                                                                                              gtagtgaggcgggCATGAAGCCAAATTTCTCGCGCTGTTCCTCTTCGTTGATAACCAGAATACCAAACACCGTCTGCTGCAGATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATa                                                           <  2:152593/263‑1 (MQ=255)                                                                                                                                                                                   gCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCCCCTTCACCGTCGTCTTCAAACATCGGTAAGTCGATAAccc                                                       >  1:152594/1‑274 (MQ=255)                                                                                                                                                                                    cGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTTCATTGCCGTCCGGCCGCCTTCACCGTCGGCTTCACACATCGGGAGTNTGATAACCcag                                                     >  2:75346/1‑274 (MQ=255)                                                                                                                                                                                                            cgcgCCGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCTCCATTATGGATACGTACTGAACCACCGCCCACGCCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGCTGCATCTTTAAGTTCTGCAGCCGTCATATCTTGCGGTGAGGTGAACGTATGGTGCATTGCCAGAAGGCCGCCATCAACGTCGTCTTCAAAAAACGGGAAGTCGATAACCAACAGACGTGCCCATGTGCTTtcgtcg                             >  2:305062/1‑275 (MQ=255)                                                                                                                                                                                                                 tGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGt                           <  1:75347/272‑1 (MQ=255)                                                                                                                                                                                                                   ttCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAg                                                                                                                                                                                           >  2:562768/1‑110 (MQ=255)                                                                                                                                                                                                                       tctTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGtt                                                                                                                                                                                       <  1:562780/110‑1 (MQ=255)                                                                                                                                                                                                                                        aaTACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCa                         >  1:292508/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaTACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGAAGAGGTGAACGNATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGCCGATAACCCACAGCGGGGCCCATTTGGTTTCGTCGGTCAGACCAAGGACTTTACCCACATTc  >  1:668681/1‑274 (MQ=255)                                                                                                                                                                                                                                            cccaccacCGTCTGCTGCATATCCCCATTATGGTTACGTACTGAACCACCCCCCACTTCGTAACCATTGATGACCATATCGTAATCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGAcc                     <  2:292503/246‑1 (MQ=255)                                                                                                                                                                                                                                               aaCACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGaa                                                                                                                        >  1:168951/1‑149 (MQ=255)                                                                                                                                                                                                                                                   ccGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGa                                                                                                                    <  2:168949/149‑1 (MQ=255)                                                                                                                                                                                                                                                      tCTGCTGCATATCACCATTATGGATAAGTACTGAACCACCGCCCACTTCGTAACCATtgatg                                                                                                                                                                                                        >  1:249043/1‑62 (MQ=255)                                                                                                                                                                                                                                                        tgctgcATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGNAACCATTGATGAACATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATa                                                           <  1:635000/201‑1 (MQ=255)                                                                                                                                                                                                                                                          ctgcATATCACCATTATGGATAAGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCa                                                                                                                                                                                                    <  2:249040/62‑1 (MQ=255)                                                                                                                                                                                                                                                                       ttATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCata                                                                                                                                         >  2:294589/1‑108 (MQ=255)                                                                                                                                                                                                                                                                           ggATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCtt                                                                                                                                     <  1:294594/108‑1 (MQ=255)                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          CAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTC  >  CP009273/1943043‑1943547

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