Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,147,792	T→A	S58T (TCG→ACG)	plsX →	putative phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,147,792	0	T	A	95.6%	124.2 / NA	46	S58T (TCG→ACG)	plsX	putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (25/18);  minor base C (1/1);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAG  >  NC_000913/1147546‑1148037                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGCCCAAACCGGGCGGCGACGATACCTTGACCCGTCTACCCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTAcgcg                                                                                                                                                                                                                                                   <  1:194906/251‑1 (MQ=255) cATCTGCGTGATGAAGCTTAGGGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGTCGACGATACCTTTACACGTCTATCCCTGGCGTAAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTAcgcg                                                                                                                                                                                                                                                   <  2:32994/251‑1 (MQ=255)              aaGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGattatt                                                                                                                                                                                                                                      >  1:153836/1‑251 (MQ=255)                                 cccAGGCAACTGGGGAAAGACCAAAACGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCCTTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCATGCGCAGTCAGTTATCg                                                                                                                                                                                                                   >  2:116455/1‑251 (MQ=255)                                  ccAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGACGTTAGATGTCATGGGAGGGGTTTTTGGCCATTCCGTGACAGTGACTGCAGCATTGCAGGCACTGAATTATAATTCGCAACTCACTCTTCTTTTATTCGGCAATTCCGACGCCATCACGCAGTTACTTGCTAAAGCTGACTTTGATCAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGc                                                                                                                                                                                                                  >  2:165410/1‑251 (MQ=255)                                     ggCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACCCGTCTAACCCTGGCGTTAGATGTCAGGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAg                                                                                                                                                                                                               <  1:177075/251‑1 (MQ=255)                                           gggggAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACATTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGc                                                                                                                                                                                                         <  2:73448/250‑1 (MQ=255)                                                aaaGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGc                                                                                                                                                                                                    >  1:124326/1‑251 (MQ=255)                                                     ccAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCg                                                                                                                                                                                               >  2:237438/1‑251 (MQ=255)                                                      cAAACCGGGCGGCGGAGATACCTTGACACGTATAACCCTGGCGTTAGCTGTCCTGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAAAAACGTACGCGTCTGCAGATTATTCCTGAGAAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGc                                                                                                                                                                                              >  2:496036/1‑251 (MQ=255)                                                               gggcgACGATACCTTGACACGTCTAACCCTGGCTTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATcc                                                                                                                                                                                     <  2:421193/250‑1 (MQ=255)                                                                       aTACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGCAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTATTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCGTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCATTCGCAAGCTATCCGCGCCAGt                                                                                                                                                                             >  1:438426/1‑251 (MQ=255)                                                                                       aaCCCTGGCGTTAGATGTCATGGGAGGGGCTTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCATACGTGGGAGTTCAATgc                                                                                                                                                             >  1:14401/1‑251 (MQ=255)                                                                                             ggCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCATCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGt                                                                                                                                                         <  2:150151/249‑1 (MQ=255)                                                                                                    gATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTgg                                                                                                                                                 >  1:416278/1‑250 (MQ=255)                                                                                                            gggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCTGCATTGCAGGCACTGAATCCTAATCCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATGCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTTGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTg                                                                                                                                         <  1:30904/250‑1 (MQ=255)                                                                                                            gggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTg                                                                                                                                         <  1:317867/250‑1 (MQ=255)                                                                                                             ggAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGTATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGGGCTgggga                                                                                                                                        >  2:497571/1‑247 (MQ=255)                                                                                                                             cTTCCGTGACAGTGCCCGAAGCATTGCAGAAACTGACTTATAATTCGCAACTCACTAATCTATTCGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTAACTTTGACCAACAACCGCGTCAGCCGCTTATCACTGCGCCGACCCTTATCGTCAAGGATGCCAGTCCTTCACAAACTATCCGCGCCCGTCGAGGGAGTTCAAAGCGCGTGTCCCTGGAGCTGGTGaaacgaagcatcactac                                                                                                                        >  2:471841/1‑236 (MQ=255)                                                                                                                                  gTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCtg                                                                                                                  >  1:202245/1‑251 (MQ=255)                                                                                                                                                       cAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCgg                                                                                             >  1:122350/1‑251 (MQ=255)                                                                                                                                                           cACTGAATTCTAATTCGCATCTCACTCTTCTTTTAGTCGGCAATTCCTACGCCATCACGCCATTACTTGCTAAAGCTGACTTTTAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCGTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGTGCTGGTGAAATAAGGTCTAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGgcg                                                                                         <  1:283147/251‑1 (MQ=255)                                                                                                                                                             cTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCACTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACAGGGGCGCt                                                                                       >  1:477122/1‑251 (MQ=255)                                                                                                                                                                aatccTAATGCGCAACCCACTCTTCTTTTATGCGGCAATTCCGACGCCATCACGCCATTACTGGCTAAAGCTGACTTTGAACAACGTACGCGTCGGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGGCGTGGGAGTTCAATGCGCGGGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGAt                                                                                    <  2:497874/247‑1 (MQ=255)                                                                                                                                                                  ttCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATgg                                                                                  >  2:122827/1‑251 (MQ=255)                                                                                                                                                                      aagtcccacctcactgTTCTTTTAGTCGGCATTTCCGGCGCCACCACGCCATTATGTGCTAAAGCTGACTTTGAACAACGTACGCTTCTGCAGATTAGTCCTGCGCAGTCAGTTTTCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGTAGTTCAATGCGCGTGGCCCTTGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGGCAGTGCCGGTAATACCGGGGCGCTGATgggtct                                                                              <  2:499085/242‑4 (MQ=255)                                                                                                                                                                      aaTTCGCAACTCCCTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTGTGAACAACGTACTCGTCTGCAGATTATTCCTGCGCAGTCATTTATCGCCAGTTATGCCCAGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAAGACCGGGGCGCTGATGGGGCt                                                                              <  2:22253/251‑1 (MQ=255)                                                                                                                                                                              actcactcTTCTTTTACCCGGCAACTCCGCCGCCCTACCCCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGCCTGCAGATTATTCCTGCGCAGGCAGTTACCGCCAGTGATGCCCGGCCTGCGCAATCCAACCGCTCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGGGCCGGTAATACCGGGGCGCTGATTGGGCTGGCAAAat                                                                      <  2:202245/251‑1 (MQ=255)                                                                                                                                                                                tcactcTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAatt                                                                     >  1:435345/1‑250 (MQ=255)                                                                                                                                                                                           ttAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAg                                                                                                                                                                                                               >  1:31542/1‑101 (MQ=255)                                                                                                                                                                                           ttAGTCGGCAATCCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAg                                                                                                                                                                                                               <  2:31542/101‑1 (MQ=255)                                                                                                                                                                                            tcgcaggCAATTCCGACCCCATCCCGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATGATTCTGGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAGGCTATCCGGTCCAGTCGTGGGAGTTCAATGCGCGTGGCCTTTGGGCTGGGGAAAGAAGGGTGACCGCAAGCCTTTGTCAGTGCCGGTGATACCGGGGGGGGGAGGGGTCTGTCAAAATTATTACGCAAGcc                                                         <  2:207348/245‑1 (MQ=255)                                                                                                                                                                                                    aaTTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATCCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAggg                                                >  2:404012/1‑251 (MQ=255)                                                                                                                                                                                                          gACGCCATGACACCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGGGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGTACGGAGGTCGAGCGCAAGCCTGTGTCCGTGCCGGTAAGACCGGGGCGATGATGGGGCTGGAAAAATTATTACTCAAGCCCCTGGAGGGGAtg                                            >  2:299019/1‑248 (MQ=255)                                                                                                                                                                                                           ccGCCATCCCGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATGCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCATGCTATCCGCGCCAGTCGTGGGAGTTCAATGCTCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGTGCTGGCAAAATTATTCCTCAAGCCCCTGGAGGGGATTGAg                                         <  1:445437/250‑1 (MQ=255)                                                                                                                                                                                                                  cACGCCATTACTTGCTAAAGCTGACTTTGACCAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGACTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCTAGCCACTGGAGGGGGTTGAGCGTCCgg                                  >  2:349808/1‑251 (MQ=255)                                                                                                                                                                                                                  cACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCAGGCCTTCGCAAGCTATCCGCGCAAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGAGCTGATGGGGCTGGCAAAATTATTACTCAAGACCCTGGAGGGGCTTGGGCGGCCgg                                  >  2:349809/1‑251 (MQ=255)                                                                                                                                                                                                                    cGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTACTCCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGGTTGAGCGTCCGgcg                                >  2:11739/1‑251 (MQ=255)                                                                                                                                                                                                                    cGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCAAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAACGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCCAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCgg                                  >  2:11753/1‑249 (MQ=255)                                                                                                                                                                                                                       cATTACTTGCTAAAACTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTg                             >  1:372273/1‑251 (MQ=255)                                                                                                                                                                                                                           cctggCTAAAGCTGACTTTGAACAACGTCCGCGTCGGCAGATTATTCCTGCGCAGTCAGTTACCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAAGGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTg                          <  2:210265/246‑1 (MQ=255)                                                                                                                                                                                                                                     gCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAcc               >  2:402413/1‑251 (MQ=255)                                                                                                                                                                                                                                     gCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAcc               >  1:145174/1‑251 (MQ=255)                                                                                                                                                                                                                                              gAACAACTTACGCGTCTGCAGATTATTCCTGGCCAGTCAGTTAACTCCATTGATGCTCGTCCGTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTTGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGGGTCAGTGCCGGTAAGACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGc       <  1:321410/250‑1 (MQ=255)                                                                                                                                                                                                                                                  aacGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCaaa    >  1:379054/1‑249 (MQ=255)                                                                                                                                                                                                                                                  aacGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAg  <  1:437510/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAG  >  NC_000913/1147546‑1148037

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