Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,185,679	C→T	R38K (AGA→AAA)	U712_10950 ←	Hypothetical Protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,185,679	0	C	T	97.1%	98.0 / ‑4.3	34	R38K (AGA→AAA)	U712_10950	Hypothetical Protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/16);  total (17/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 = 8.06e-01

ATTCAGACGAAAGCCCCCGGCAAATAAGATGCTTGCTGTACGGCAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGTTCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAACAGCCCTCTGACCTGCCGTTTTGTA‑‑TAAAGCGCATCCCCCCTTTTCGTTTTTTCTTTTCTCTCTATAAGACAAGCTATGTTTTTTTGCCTCGGATTATGATTCTC  >  CP006881/2185544‑2185776                                                                                                                                        |                                                                                                    aTTCAGACGAATGCCCCCGGCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGAc                                                                                                 >  4:316446/1‑139 (MQ=255)          aaaGCCCCCGGCAAATAAGATGCTTGCTGTACGGCAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGTTCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAACAGCCCTCTGACCTGCCGtt                                                                                         <  6:101749/139‑1 (MQ=255)           atgCCCCCGGCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCggctt                                                                                       <  8:25456/137‑5 (MQ=255)              cccccGGCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGt                                                                                     <  5:72504/138‑1 (MQ=255)                caccgCAAATAAGATGCTT‑CTGGACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGtaag                                                                                  <  6:609/135‑3 (MQ=255)                 ccGGCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGtaag                                                                                  >  2:134920/1‑136 (MQ=255)                  cGGCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGtaagta                                                                                <  4:303727/139‑5 (MQ=255)                   ggCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGt                                                                                     <  8:168507/133‑1 (MQ=255)                   ggCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGt                                                                                     >  7:168507/1‑133 (MQ=255)                    gCAAATAAGATGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGtaagtaaa                                                                              <  2:35055/139‑7 (MQ=255)                            gATGCTT‑CTGTACGGAATGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATcc                                                                      <  7:175680/139‑1 (MQ=255)                             aTGCTT‑CTGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATccc                                                                     >  8:309738/1‑139 (MQ=255)                                 cttctGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATcccccc                                                                  >  8:202496/4‑139 (MQ=255)                                   tctGTACGGAAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCtt                                                                >  8:366414/2‑139 (MQ=255)                                            aaGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCGttttt                                                       >  4:253231/2‑139 (MQ=255)                                            aaGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCGttttt                                                       >  5:83651/2‑139 (MQ=255)                                            aaGGGAATTATTGAGTCTTATTAGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCGttttt                                                       >  8:388446/2‑139 (MQ=255)                                                               attaGTTTACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGa                                     <  6:83651/137‑1 (MQ=255)                                                                     ttACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTa                              >  6:80400/1‑138 (MQ=255)                                                                     ttACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTa                              <  5:80400/138‑1 (MQ=255)                                                                     ttACCAAAAGATCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTAt                             >  3:157767/1‑139 (MQ=255)                                                                              gatCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGttttttcg                    >  5:337110/3‑137 (MQ=255)                                                                              gatCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGttttttcg                    >  2:392884/3‑137 (MQ=255)                                                                                 cTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATACGTTCTATAAGACAAGCTATGttttttctcct                 >  3:413169/1‑134 (MQ=255)                                                                                    gagaTAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGcc                  <  5:291619/135‑1 (MQ=255)                                                                                    gagaTAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGcc                  <  7:140423/135‑1 (MQ=255)                                                                                    gagaTAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGcc                  <  7:202496/135‑1 (MQ=255)                                                                                    gagaTAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGcc                  <  7:309738/135‑1 (MQ=255)                                                                                         aGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTAt         >  4:318150/1‑139 (MQ=255)                                                                                                 aaaaaaGGAAGATCGCCCTCTGTACAGTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGAtctgc  <  1:392884/138‑5 (MQ=255)                                                                                                                           gTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGcc                  <  2:45223/96‑1 (MQ=255)                                                                                                                           gTAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGcc                  >  1:45223/1‑96 (MQ=255)                                                                                                                            tAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGa                                     >  4:193921/1‑76 (MQ=255)                                                                                                                            tAAAGCAGCTCTTTGACCTGCCGGCTTGTAAGTAAAGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGa                                     <  3:193921/76‑1 (MQ=255)                                                                                                                                        |                                                                                                    ATTCAGACGAAAGCCCCCGGCAAATAAGATGCTTGCTGTACGGCAGGGAATTATTGAGTCTTATTAGTTTACCAAAAGTTCTTGAGATAGAAAATGAAAAAAGGAAGATCGCCCTCTGTACAGTAAAACAGCCCTCTGACCTGCCGTTTTGTA‑‑TAAAGCGCATCCCCCCTTTTCGTTTTTTCTTTTCTCTCTATAAGACAAGCTATGTTTTTTTGCCTCGGATTATGATTCTC  >  CP006881/2185544‑2185776

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