Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,829,887	G→A	A293A (GCC→GCT)	astA ←	arginine N‑succinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,829,887	0	G	A	97.7%	135.3 / ‑6.3	45	A293A (GCC→GCT)	astA	arginine N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (26/16);  minor base C (1/0);  total (28/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 = 9.47e-01

GCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAG  >  NC_000913/1829647‑1830130                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    gccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCCCCAGCACCCCGCGGGA‑ATGGTGATAATTTTTATTAGCGACcaggc                                                                                                                                                                                                                                           >  2:355352/1‑251 (MQ=255)  ccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACcaggca                                                                                                                                                                                                                                          <  2:59925/251‑1 (MQ=255)       acacCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAAAGCTCGGTTGCCGGAACGGTACGCACCCGCACCACGCGGAA‑ATGGTGATAATTTTCCTTAGCGACCAGGCAGGCTg                                                                                                                                                                                                                                     >  2:350678/1‑251 (MQ=255)        cacctagaccGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATACATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGGATCGGTACGCCACAGCACCACGCGGGA‑ATGGGGATAATTTTCATTCGCGACCcggcaggccgg                                                                                                                                                                                                                                    >  2:389128/9‑248 (MQ=255)                              cGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGATAAAATCAAACGCTCGGTTGCCGGTTCGGCACGAACCAGCACCACCCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCCGGGAAATCGCCCTGCGCAGGCTGc                                                                                                                                                                                                              >  1:185604/1‑251 (MQ=255)                                   cacaCGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTc                                                                                                                                                                                                         <  1:77206/251‑1 (MQ=255)                                         cgATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTAAGCACCCGCACCACGCGGAA‑ATGGGGATAATTTTCATTAGCGCCCAGGAAGGCTGGGAAATCGACCTGCGCCAGATGCCCTTCTGCCAc                                                                                                                                                                                                   >  1:331024/1‑251 (MQ=255)                                         cgATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGGATCGGGACGCAACAGCAACACGCGGAA‑ATGGGGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCCAGCTGCCCCTCCGCCAc                                                                                                                                                                                                   >  1:331030/1‑251 (MQ=255)                                                         gTTAGCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGGGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGGGGAA‑ATGGTGATAATTTTCATTGGCGACCAGGAAGGCTGGGAAATCGCCCTGCGCAGGAGGCCCGTCGGCCACTTCCACCAGCCGACtt                                                                                                                                                                                   <  2:152419/251‑1 (MQ=255)                                                          ttATCCAGTCACCGTTAATCCTTAAAGTCATTCTGTTTTCTCCTCTGCGCCCAGGCGAAACAGACGAAAGCGATCCCCGGCGTGGCATTTGCGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTAGGTTGCCGGATCGGTACGCACCAGCCCCACGCGGAACATG‑CGCTAAATTTCATTAGCGACCAGGCAGGCCGGGGAATCAGCCAGCGGAGGCTGCCCTTCTgacccctccaaccacagactct                                                                                                                                                                                  >  2:428207/1‑230 (MQ=255)                                                            aTCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCCAACGCTCGGTTGCCGGATAGGTACGCACCCGCACCACGCGGAA‑ATGGGGATAATTTTCCTTAGCGACCAGGCAGGCTGGGAAATAGCCCTGGGCAGGCGGCCCTTCTGCCACTTCCCCCAGCCCACtttt                                                                                                                                                                                 >  2:197101/1‑250 (MQ=255)                                                                 gTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGGGATAATTTTCATTAGCGACCAGGCAGGCTGGGGAAACGCCCTGGGCCGTCTGCCCTTCTGCCACTTCCACCCGCCGACTTTTCCGGAt                                                                                                                                                                           >  2:94804/1‑251 (MQ=255)                                                                  tCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGGACGCACCAGCACCCCGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATg                                                                                                                                                                          >  1:206047/1‑251 (MQ=255)                                                                       gTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCACACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTTGCGACCAGGCAGGCTGGGAAAACGCCCTGCGCAGGCTGCCCCTCTTCCCCTTTCACCCGCCGACTTTTTCGGATcgcgcg                                                                                                                                                                     >  2:23566/1‑251 (MQ=255)                                                                             ccATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAcgcg                                                                                                                                                               >  1:70423/1‑251 (MQ=255)                                                                                             ttttCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCCCGCGGTCGATGTCACACTc                                                                                                                                                >  1:409627/1‑250 (MQ=255)                                                                                                tctcCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGGGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCCCGCGGTCGATGTCACACTCCAGc                                                                                                                                            >  1:56205/1‑251 (MQ=255)                                                                                                tctcCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGGCTTTTAAGGATGGCGCGCACGCGGTCGATGTCACCCCCAAGc                                                                                                                                            >  1:183555/1‑251 (MQ=255)                                                                                                                                 gATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTa                                                                                                           <  1:197101/251‑1 (MQ=255)                                                                                                                                      ccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATTTCGATATAGTTACGGTa                                                                                                      <  1:133419/251‑1 (MQ=255)                                                                                                                                      ccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTa                                                                                                      <  1:133425/251‑1 (MQ=255)                                                                                                                                            tGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGaa                                                                                                 >  1:164664/1‑250 (MQ=255)                                                                                                                                              gCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAAGCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGGCACGCTCAAGCGTCGGCTCACCGTGAAAGATTAAGATATAGGTACGGTAGCGAAAAc                                                                                              <  2:29342/251‑1 (MQ=255)                                                                                                                                                 tttGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGGAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACCCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACctt                                                                                           >  2:102969/1‑251 (MQ=255)                                                                                                                                                        gCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATACAGTTACGGTAGCGAAAACCTTCTTtctc                                                                                    >  1:334295/1‑251 (MQ=255)                                                                                                                                                                 tgtgCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCCCGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGGAGCGAAAACCTTCTTTCTCCAGCCCcgc                                                                           >  1:286410/1‑251 (MQ=255)                                                                                                                                                                              aaaTCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGcgcg                                                              >  1:262890/1‑251 (MQ=255)                                                                                                                                                                                           ggTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                  >  2:384108/1‑250 (MQ=255)                                                                                                                                                                                                      cGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                                      >  2:136226/1‑251 (MQ=255)                                                                                                                                                                                                      cGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATACAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGAt                                       <  2:334295/250‑1 (MQ=255)                                                                                                                                                                                                          aCGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGt                                  <  2:144151/251‑1 (MQ=255)                                                                                                                                                                                                             caccagcaccaCGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCt                               <  2:409627/251‑1 (MQ=255)                                                                                                                                                                                                              accagcaccaCGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTg                              <  1:165990/251‑1 (MQ=255)                                                                                                                                                                                                                     ccacGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTcc                       >  2:134629/1‑251 (MQ=255)                                                                                                                                                                                                                       acGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTg                     <  1:200545/251‑1 (MQ=255)                                                                                                                                                                                                                       acGCGGAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTg                     <  1:200522/251‑1 (MQ=255)                                                                                                                                                                                                                           ggAA‑ATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAt                 >  1:268721/1‑251 (MQ=255)                                                                                                                                                                                                                             aa‑aTGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGCAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCTGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATa                >  1:91980/1‑250 (MQ=255)                                                                                                                                                                                                                              a‑aTGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaa               >  1:154407/1‑250 (MQ=255)                                                                                                                                                                                                                                  ggTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaaaa             >  1:286128/1‑249 (MQ=255)                                                                                                                                                                                                                                  ggTGATAATTTTCATTAGCGACAAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCGCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAg           <  2:69199/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGTCTGCCCTTCTGCCACTTCCCCCAGCCGACTTGTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAATCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGt          <  2:378506/251‑1 (MQ=255)                                                                                                                                                                                                                                     gATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTgg        <  1:389128/251‑1 (MQ=255)                                                                                                                                                                                                                                     gATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTgg        <  2:397882/251‑1 (MQ=255)                                                                                                                                                                                                                                           tttCATTAGCGACAAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAACGTGGGtatcg  >  1:330170/1‑249 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    GCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAA‑ATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAG  >  NC_000913/1829647‑1830130

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