Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,765,463	Δ1 bp	7.1%	coding (454/1242 nt)	U712_19410 →	Tyrosyl‑tRNA synthetase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,765,463	0	G	.	7.1%	184.0 / 8.3	56	coding (454/1242 nt)	U712_19410	Tyrosyl‑tRNA synthetase 2
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base . (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CTGACTGATGAAGAAGTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAGGGCTATGATTCTGTCGTTCTCGA  >  CP006881/3765331‑3765599                                                                                                                                     |                                                                                                                                          ctgactgaTGAAGAAGTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAAt                                                                                                                                        >  1:44536/1‑136 (MQ=255)    actgaTGAAGAAGTTCAGCACAATGCCAAAACCTACTTTGAGCAATTAGGAAAAGTGCATGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCa                                                                                                                                   <  1:120879/138‑1 (MQ=255)        atgaAGAAGTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGa                                                                                                                               >  8:93093/1‑138 (MQ=255)        atgaAGAAGTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGa                                                                                                                               >  4:253657/1‑138 (MQ=255)           aagaagTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGaaaa                                                                                                                            <  2:44536/138‑1 (MQ=255)                gTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATATAGAAGATGTCATTGAATTAGCAGGGAAAATAACgg                                                                                                                      >  4:108338/1‑139 (MQ=255)                      cacaATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACTGTAGacc                                                                                                                >  2:145841/1‑139 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGc                                                                                                              <  4:257454/136‑1 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGc                                                                                                              >  3:257454/1‑136 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCt                                                                                                            >  8:227875/1‑138 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCt                                                                                                            >  7:3260/1‑138 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCt                                                                                                            >  6:49069/1‑138 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCt                                                                                                            >  3:142881/1‑138 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCt                                                                                                            >  3:228176/1‑138 (MQ=255)                           tGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCt                                                                                                            >  1:45264/1‑138 (MQ=255)                                aaaCCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATg                                                                                                        >  1:73136/1‑137 (MQ=255)                                     tACTTTGAGCAATTCGGAAAAGTGCTTGCTCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATAAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAgcgc                                                                                                  >  5:85/1‑138 (MQ=255)                                     tACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAgcgc                                                                                                  >  7:118068/1‑138 (MQ=255)                                       cTTTGAGCAATTCTGAAAAGTGCTAGATCCATAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGTAAAATAACGTTAGCCCGCCTGATGGAGCGcg‑a                                                                                               <  1:145841/138‑1 (MQ=255)                                       cTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGcg‑ac                                                                                              >  7:114212/1‑139 (MQ=255)                                       cTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGcg‑ac                                                                                              >  6:91390/1‑139 (MQ=255)                                       cTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGcg‑a                                                                                               <  8:10243/138‑1 (MQ=255)                                           gAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATT‑AATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACtt                                                                                         >  7:142507/1‑139 (MQ=255)                                           gAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATT‑AATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACtt                                                                                         >  7:138650/1‑139 (MQ=255)                                            aGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCa                                                                                                                                   <  6:3177/98‑1 (MQ=255)                                            aGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCa                                                                                                                                   >  5:3177/1‑98 (MQ=255)                                            aGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACtt                                                                                         >  4:190943/1‑139 (MQ=255)                                              cAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTg                                                                                       <  4:228176/139‑1 (MQ=255)                                                aatctgAAAAGAGCTTGATCCAGAAAAAGTAGAGATTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTgaa                                                                                     <  6:185855/134‑1 (MQ=255)                                                 ttCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTa                                                                                                                    >  1:30227/1‑108 (MQ=255)                                                 ttCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTa                                                                                                                    <  2:30227/108‑1 (MQ=255)                                                           tGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGCCTTTGAAGAACGCATCGc                                                                          >  8:240144/1‑139 (MQ=255)                                                           tGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGc                                                                          >  6:15647/1‑139 (MQ=255)                                                             cTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGcc                                                                         <  3:190943/138‑1 (MQ=255)                                                             cTTGATACAGAAAAAGGCGAGCGTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGcc                                                                         <  6:234619/138‑1 (MQ=255)                                                                 aTCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCa                                                                    <  5:49069/139‑1 (MQ=255)                                                                     aGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCaaaaa                                                                >  4:200463/1‑139 (MQ=255)                                                                       aaaaaGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAAcc                                                              <  8:118068/139‑1 (MQ=255)                                                                        aaaaGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCa                                                             <  3:200463/139‑1 (MQ=255)                                                                          aaGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAAt                                                           <  4:142881/139‑1 (MQ=255)                                                                                                aaaaTGGCTGAAAACCTTGAATCTAGAAGATGTCATTGAATTTGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTa                                      <  5:149213/138‑1 (MQ=255)                                                                                                aaaaTGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTa                                      <  8:114212/138‑1 (MQ=255)                                                                                                aaaaTGGCTGAAAACATTGAATCTAGAAGATGTCATT‑AATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTAcc                                    <  8:142507/139‑1 (MQ=255)                                                                                                aaaaTGGCTGAAAACATTGAATCTAGAAGATGTCATT‑AATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTAcc                                    <  8:138650/139‑1 (MQ=255)                                                                                                  aaTTGCCGAAAACCATGAATCTATAACAAGTCAATTAAATAGCACGGAAAATAACGGGACCCCGCCAGCTGGAGCGCGCACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTAcc                                    <  2:73136/139‑1 (MQ=255)                                                                                                  aaTGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTAcc                                    >  8:261819/1‑138 (MQ=255)                                                                                                  aaTGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTAcc                                    >  7:238527/1‑138 (MQ=255)                                                                                                    tGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCAt                                 >  4:215471/1‑139 (MQ=255)                                                                                                    tGGCTGAAAACAATGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAAAAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCAt                                 <  3:231957/139‑1 (MQ=255)                                                                                                          aaaaCATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATg                            <  3:253657/138‑1 (MQ=255)                                                                                                            aaCATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAg                         <  7:240144/139‑1 (MQ=255)                                                                                                                  gAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAGGGCTa                    <  5:91390/138‑1 (MQ=255)                                                                                                                          aagaTGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTTCATGAATTCTTTTACCCATTGATGCAGGGCTATGATTCTg            <  3:108338/138‑1 (MQ=255)                                                                                                                               gTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAGGGCTATGATTCTGTCGTtc      >  4:114205/1‑139 (MQ=255)                                                                                                                                   ttGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAGGGCTATGATTCTGTCGTTCTCGa  <  5:15647/139‑1 (MQ=255)                                                                                                                                    tGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAGGGCTATGATTCTGTCGTTCTCGa  >  2:277752/1‑138 (MQ=255)                                                                                                                                     |                                                                                                                                          CTGACTGATGAAGAAGTTCAGCACAATGCCAAAACCTACTTTGAGCAATTCGGAAAAGTGCTTGATCCAGAAAAAGTCGAGCTTCACTATAACTCAAAATGGCTGAAAACATTGAATCTAGAAGATGTCATTGAATTAGCAGGGAAAATAACGGTAGCCCGCCTGATGGAGCGCG‑ACGACTTTGAAGAACGCATCGCCATGCAAAAACCAATCTCACTGCATGAATTCTTTTACCCATTGATGCAGGGCTATGATTCTGTCGTTCTCGA  >  CP006881/3765331‑3765599

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