Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	760,166	C→T	18.2%	Y487Y (TAC→TAT)	sucA →	2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	760,166	0	C	T	18.2%	126.5 / 20.3	55	Y487Y (TAC→TAT)	sucA	2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base T (6/4);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

CAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATG  >  NC_000913/759869‑760457                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                    caacaaCCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTAc                                                                                                                                                                                                                                                                                                     >  1:255509/1‑298 (MQ=255)                               tCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCATGCCCCGATTTGCCACGTTAACGCTGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCa                                                                                                                                                                                                                                                                                    <  2:208004/285‑1 (MQ=255)                                                 tACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGAAGGAAAAAGTGGCGACGCTGGAAGATGCCACCg                                                                                                                                                                                                                                                    >  1:179786/1‑299 (MQ=255)                                                    gCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCgaga                                                                                                                                                                                                                                                 >  1:333802/1‑299 (MQ=255)                                                     ccGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGAt                                                                                                                                                                                                                                                >  1:274702/1‑299 (MQ=255)                                                      cttactgtactgATATCGGTAAGATGGTTCAGGCCCCGTTTTTCCACGTTAACGCGTACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATGTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGaa                                                                                                                                                                                                                                                              <  2:480172/282‑1 (MQ=255)                                                       tccgatctctgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATgg                                                                                                                                                                                                                                              <  1:386808/291‑1 (MQ=255)                                                           tgtactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGCGCCGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTa                                                                                                                                                                                                                                           >  1:69944/1‑298 (MQ=255)                                                               ctgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTaa                                                                                                                                                                                                                                          >  1:516116/1‑295 (MQ=255)                                                                tgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGATGGAGCAGGAACAAGTGGCGACGCTGGAAGAGGcaccaga                                                                                                                                                                                                                                                   >  2:339068/1‑279 (MQ=255)                                                                           aGATGGTTCGGGCCCCGCGTGGACGCGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTTGCGCTCGATTTCCGTAACACCTTTAAACGAGATGTCTCCATCGACCTGTTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAAc                                                                                                                                                                                                                                         <  2:349077/284‑1 (MQ=255)                                                                                  tCAGGCCCCGATTTTCCACGTTAACGCGGACGCTCCGGAAGCCGTTGGCTTTGTGACCCGTCTGGCGCTCGATGTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCAGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGAt                                                                                                                                                                                                                    <  1:289800/298‑1 (MQ=255)                                                                                     ggCCCCGATTTTCAACGTGAACGCGGACGATCCGGAAGCCGTTGCCTTGGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCGGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACcgcg                                                                                                                                                                                                                               <  2:86660/284‑1 (MQ=255)                                                                                         ccGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGCAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGCAAAGCTGGAGCAGGAAAAAGTGGCGACGCTGCAAGATGACAACGAGAAGGTTAACCTGTACCGCGATg                                                                                                                                                                                                                            >  2:342785/1‑283 (MQ=255)                                                                                                tCCACGTTAACTCGGCCGATCCGGAAGCCGTTTCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGAt                                                                                                                                                                                                                    <  2:293117/284‑1 (MQ=255)                                                                                                   aCGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTg                                                                                                                                                                                                  >  1:518327/1‑299 (MQ=255)                                                                                                     gTTAACGCGGATGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGCGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTgg                                                                                                                                                                                                 >  1:62115/1‑298 (MQ=255)                                                                                                     gTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGc                                                                                                                                                                                                >  1:234000/1‑299 (MQ=255)                                                                                                                  aTCCTGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACTGGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTgg                                                                                                                                                                                                 <  2:225055/285‑1 (MQ=255)                                                                                                                    ccGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGa                                                                                                                                                                                  <  1:468669/298‑1 (MQ=255)                                                                                                                       gAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCGGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTAACGAGATGCGCTGGATGCTGGCGATTGCGTAGTGGCagg                                                                                                                                                                                             >  2:23426/1‑283 (MQ=255)                                                                                                                             gTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGc                                                                                                                                                                         >  1:46012/1‑298 (MQ=255)                                                                                                                                     tgtgACCCGTCTGGCGCTCGATTTCCGTACCACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCa                                                                                                                                                                <  1:63843/299‑1 (MQ=255)                                                                                                                                         aCCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTg                                                                                                                                                            >  1:185114/1‑299 (MQ=255)                                                                                                                                         aCCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTg                                                                                                                                                            >  1:212705/1‑299 (MQ=255)                                                                                                                                             gTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGGGTGGCGTCCGATGAACATGCAc                                                                                                                                                                       >  2:48555/1‑284 (MQ=255)                                                                                                                                                    gcTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTat                                                                                                                                                  >  1:40864/1‑297 (MQ=255)                                                                                                                                                       cGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTATCTCa                                                                                                                                              >  1:592135/1‑299 (MQ=255)                                                                                                                                                          tttCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGAGGACCATGCACCCTTTCCCCTGGt                                                                                                                                                          >  2:118328/1‑284 (MQ=255)                                                                                                                                                          tttCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGGCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAAc                                                                                                                                            <  1:561542/298‑1 (MQ=255)                                                                                                                                                          tttACGTAACACCTTTAAACGTGAGGTCTTCAGGGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAGGCCGACGGGCCGAGCGGAACCGAGCCGCTGATGTATCAAAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTc                                                                                                                                                         <  2:29699/285‑1 (MQ=255)                                                                                                                                                                        gtAAACGGGATGTCTTCATCCACCTGGTGGGCTCCCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAAc                                                                                                                                            <  2:158299/283‑1 (MQ=255)                                                                                                                                                                           aaCGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAagag                                                                                                                          <  1:180519/299‑1 (MQ=255)                                                                                                                                                                              gTGATGTCTTCATCGACCGGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGGGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAACCATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAAt                                                                                                                                     <  2:452473/285‑1 (MQ=255)                                                                                                                                                                                   ttcttctTCGAACTGGTGTGGTACCGCCGTCACGTCCACAACGAAGCTGACGAGCCGAGCTCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGCAAAAGTGTCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATggg                                                                                                                                  <  2:60649/282‑1 (MQ=255)                                                                                                                                                                                    tCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCTGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAAc                                                                                                                 >  1:394928/1‑299 (MQ=255)                                                                                                                                                                                       tCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTATCTCAACCACGAATGGGACGAag                                                                                                                            <  2:162871/285‑1 (MQ=255)                                                                                                                                                                                         aTCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGt                                                                                                            >  1:81941/1‑299 (MQ=255)                                                                                                                                                                                                  gtgtGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGAAGCTGGAAGATGCCACCGAGATGGTTAACCTGGACCGCGATGCGCTGGATGCTGGCGATTGCGCAGTGGCAGAGTGGAGTGCGATGAACAAGCACTCTGTCACCTGGACGCCGTACCACAAACACGAATGGGAAGAagagatcggaaga                                                                                                                 >  2:379801/1‑276 (MQ=255)                                                                                                                                                                                                        tACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACTTCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACaaa                                                                                                              >  2:13829/1‑282 (MQ=255)                                                                                                                                                                                                        tACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGtt                                                                                                           <  2:450504/285‑1 (MQ=255)                                                                                                                                                                                                        tACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGACGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTATCTCAACCACGAATCGGACGAAGAGTACCCGAACAAAGTTGAGATGAAgcg                                                                                                >  1:49332/1‑296 (MQ=255)                                                                                                                                                                                                          ccgccgTCACGGCCAAAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCTGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTga                                                                                                         <  2:345971/285‑1 (MQ=255)                                                                                                                                                                                                           cgccgTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCAACGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCCCGAATGGGACGAAGAGTACCCGAACAAAGTTgag                                                                                                        >  2:328712/1‑285 (MQ=255)                                                                                                                                                                                                                        acaACGACGCGGACGTGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCGATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGGTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTATCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGc                                                                                           <  2:515324/285‑1 (MQ=255)                                                                                                                                                                                                                           aCGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGc                                                                          <  1:257920/299‑1 (MQ=255)                                                                                                                                                                                                                             gAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACTTCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCAt                                                                        <  1:17025/299‑1 (MQ=255)                                                                                                                                                                                                                                gCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTATCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTg                                                                                   <  2:34869/285‑1 (MQ=255)                                                                                                                                                                                                                                          caaacgcAACACAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCTAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCATACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGcc                                                            <  1:600515/294‑1 (MQ=255)                                                                                                                                                                                                                                                         cGCTGATGTATCAGATAATCAAAAAACATCCGACACCGCGCAAAATCGACGCTGACTAGCTGGAGCAGGAAAAAGTGGCGATGCTGGAAGATGCCACCGAGTGGGTTAACCTGTACCGCGATGCGCTGGATGCTGTCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTGGAAAt                                             <  1:297742/298‑1 (MQ=255)                                                                                                                                                                                                                                                         cGCTGATGTATCAGAAAATCAAAAAACATCCGACGCCGCGCAAAATCTATGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGAAACTCAACCACGAATGGGACGAAGAGTACGCGAACAAAGTTGATATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTAGaact                                             >  1:526777/1‑296 (MQ=255)                                                                                                                                                                                                                                                                     aGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGACACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGtt                                >  1:460397/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                aaaCATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGCAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGCTGAACATGCACTCTTTCACCTGGTCGCCGTACCCCAACCACGAATGGGTCGAAGAGTACCCGCACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAACCGCCTCAGCACGGTGCCAGAAGCTGTTGACATgacgtacgcg                                   >  2:69093/1‑276 (MQ=255)                                                                                                                                                                                                                                                                                aaaCATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTcgc                                   <  2:187197/285‑1 (MQ=255)                                                                                                                                                                                                                                                                                        gACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATc             >  1:227082/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                                   aaaTCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATg  <  1:27679/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                                          cGCTGACAAGGTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCa         <  2:39551/285‑1 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                    CAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGCGGACGATCCGGAAGCCGTTGCCTTTGTGACCCGTCTGGCGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCTGGTGTGCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCAGAGTGGCGTCCGATGAACATGCACTCTTTCACCTGGTCGCCGTACCTCAACCACGAATGGGACGAAGAGTACCCGAACAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATG  >  NC_000913/759869‑760457

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