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	98.0%	167.2 / NA	50	A293A (GCC→GCT)	astA	arginine N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/34);  minor base T (1/0);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATC  >  NC_000913/1829646‑1830133                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       tgccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGGGATAATTTTCATTAGCGACcagg                                                                                                                                                                                                                                               >  1:454046/1‑251 (MQ=255)  gccATAACCCCGCGCCCGGGACCGGATTACGCTGCACCCGCGATGCGCCCTGGCCCGTGATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCGCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGTCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTTATAATTTTCATTAGCGTCcaggc                                                                                                                                                                                                                                              <  2:98179/251‑1 (MQ=255)  gccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACcaggc                                                                                                                                                                                                                                              <  1:195156/251‑1 (MQ=255)    cATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGCTGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCCAACGCCCGGTTGCCGGATCGGGACGCACCAGCCCCACGCGGAAATGGTGATAATTTTCATTAGCgccccggccg                                                                                                                                                                                                                                            >  1:155375/1‑242 (MQ=255)               gCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCCCGGTTGCCGGAACGGGACGCACCCGCACCACGCGGAAATGGGGATAATTTTCATTAGCGCCcaggcacggtgggcaat                                                                                                                                                                                                                                  >  2:343851/1‑239 (MQ=255)                      aCCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGGGTCGGTACGCACCAGCAACACGCGGAAATGGGGATAATTTTCATTTGCGACCCGGCAGGCTGGGAAATCGccccgc                                                                                                                                                                                                                          >  2:62079/1‑248 (MQ=255)                         ggATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGGTAAAATCAAACCCCCGGTTGCCGGATCGGTACGCACCAGCACCACCCGGCAATGGTGATAATTTTCTTTAGCGACcaggcacgcctggaaacacgcccgcgcc                                                                                                                                                                                                                       >  2:386361/1‑229 (MQ=255)                                  ttCACACGCGATGCGCCCTGGCCAGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGAACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCGGCACCACGCGGCAATGGGGATAATTTTCATTAGCGACCAGGCAGGCTGGGAGATCGCCCTGAGCAGGGTGCCCt                                                                                                                                                                                                              >  1:398563/1‑251 (MQ=255)                                       acGCGATACGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCG‑CGGGTAAAAACAAACGCCCGGGTGCCGGATCGGTACGCCCCAGCACCACGCGGAAAAGGTGATAATTTTCATTATCGACCAGGCCGGCTGGGGAAACACCCTGCGCAGGCTGCCCtttttac                                                                                                                                                                                                        >  2:453966/1‑246 (MQ=255)                                           gATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACt                                                                                                                                                                                                     <  1:339615/251‑1 (MQ=255)                                                  ccttgCCCGTTATCCAGTCACCGTTAATCCAGAAAGTCATGCTGTTTCCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCTTTGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCCGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTccacca                                                                                                                                                                                              <  2:80628/247‑1 (MQ=255)                                                      gCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCGGCCACTTCCACCAGCCg                                                                                                                                                                                          <  2:448286/251‑1 (MQ=255)                                                        ccgttgtacAGTCTCCGTTAAGACAGAAAGTCTGGCTGTGTTCTCCGCTGCTCACAGGCGCACCAGCCGAACGGGAGCCCCTGCGTGGCATGTGAGGGCATCCAGTTGTGCGGCGGTTAAAAGCAAACGCTCGGTTGCCGTATCGGGACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGTCTGGGGAATGGCCCTGATCAGCCTGCCCTTCTGCCACTGCCACCATCCGAc                                                                                                                                                                                        <  2:381268/243‑1 (MQ=255)                                                                                  aaaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGCTAAAATCAAACGCTCGGTTGCCGGATCGGTACCCACCCGCACCCAGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAAACGACTTTTACGGATGGGGCGCACGCGGTCg                                                                                                                                                              >  1:225909/1‑251 (MQ=255)                                                                                                         gcgcACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGc                                                                                                                                         <  1:11304/249‑1 (MQ=255)                                                                                                            cacaGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCAcc                                                                                                                                    <  1:307744/251‑1 (MQ=255)                                                                                                                   gcACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAcg                                                                                                                                                                    <  1:295292/212‑1 (MQ=255)                                                                                                                   gcACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAg                                                                                                                             <  1:44598/251‑1 (MQ=255)                                                                                                                    cACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCGGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGa                                                                                                                            <  1:218168/251‑1 (MQ=255)                                                                                                                     aCCAGACGAACGCGATCCC‑GGCGGGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCGCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATg                                                                                                                          <  2:102449/251‑1 (MQ=255)                                                                                                                      ccAGACGAACGCGAACCCCGGCGTGGCAGTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCCCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCCGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGa                                                                                                                            <  2:121508/249‑1 (MQ=255)                                                                                                                        aGACGAACGCTATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCGCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTTCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCTGCCCACCGTCAAAGATGTc                                                                                                                        <  2:172959/251‑1 (MQ=255)                                                                                                                                               gCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCCGCACCACGCGGGAATGGGGATAATTTTCATTAGCGACCCGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGcc                                                                                                                                                                                                        >  1:184988/1‑148 (MQ=255)                                                                                                                                               gCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGcc                                                                                                                                                                                                        <  2:184988/148‑1 (MQ=255)                                                                                                                                                 aTTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACct                                                                                               <  1:215594/251‑1 (MQ=255)                                                                                                                                                         gCATCCAGTTGTGCGGCGGTTAAAATCAAAAGCTCGGTTGACGGATCGGGACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTtctc                                                                                       >  1:428153/1‑251 (MQ=255)                                                                                                                                                          cATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTcc                                                                                      <  1:11246/251‑1 (MQ=255)                                                                                                                                                          cATCCAGTTGTGCGGCGGTTAAAATAAAACGCTCGGTTGCCGGATGGGTACGCACCAGCACCACGCGGAAATGGGGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGACCTTCTGCCACTTCCACCAACCGCCTTTTACGGATGGCGCGCACGAGGTCGATGTCACACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATATAGTTACGATAGCGAAAACCTTCTTTTTcc                                                                                      >  1:24521/1‑251 (MQ=255)                                                                                                                                                            cccAGTTGTGCGGCGGTTAAAATCAAACGCGCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGGCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAg                                                                                    <  2:454046/250‑1 (MQ=255)                                                                                                                                                             ccAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTGCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAg                                                                                    <  2:177953/250‑1 (MQ=255)                                                                                                                                                                         cgggtAAAATCAAACGCTCGGTGGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCGGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTGGATGTCACACTCAATCGTCGGCCCACCGTCAAAGATGTCGATAGAGTTACGGTAGCGAAAACCTGCTTTCTCCAGCTCCGCGCGGGCa                                                                       <  2:225909/247‑1 (MQ=255)                                                                                                                                                                          ggTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAg                                                                                                                                                                                             <  1:169447/132‑1 (MQ=255)                                                                                                                                                                          ggTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAg                                                                                                                                                                                             >  2:169447/1‑132 (MQ=255)                                                                                                                                                                                     aaCGCTCTGGTGCCGGATCGGTTCGCACCAGCCCCACGCGGAAATGTTGATAATTTTCATGAGCGCCCAGGCAGGCTGGGAAATCGCCCTGCGGAGGCTGCCCTTCGGACCCGTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACGCTCGAGCGACGGCCCACCGTCAAAGATGTCGATATATGTGCGGTTGCGAGATGCTTCTTTCTCCAGCACCGCGCGTGCAGGCGCGGTTTGc                                                           <  1:65195/251‑1 (MQ=255)                                                                                                                                                                                                  cGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAc                                              <  1:62079/251‑1 (MQ=255)                                                                                                                                                                                                    gATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGt                                                                                                                                  >  1:385658/1‑165 (MQ=255)                                                                                                                                                                                                    gATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGt                                                                                                                                  <  2:385658/165‑1 (MQ=255)                                                                                                                                                                                                     aTCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCg                                            >  1:104185/1‑250 (MQ=255)                                                                                                                                                                                                           gcGCACCAGCAACACGCGGAAAGTGTGATAATTTTCATGAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGGCGGCCCACCGTCAAAGATGTCGAGATAGTTACGGTAGCGAAAACCTTCTTTCGCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACg                                      <  2:6435/249‑1 (MQ=255)                                                                                                                                                                                                               cccagcaccaCGCGGAAATGTTGATAATTTCCACTAGCGACCAGGCAGGCTGGGAAATCGCCCGGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCTTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTGCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGTATGTACCTGACCGATGGCGTCCTg                                 <  2:104790/250‑1 (MQ=255)                                                                                                                                                                                                                ccagcaccaCGCGGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCGGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTgg                                <  2:155375/251‑1 (MQ=255)                                                                                                                                                                                                                         cgcgGAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTgg                       <  2:348002/251‑1 (MQ=255)                                                                                                                                                                                                                             gAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACgg                                                                                                           <  2:248592/163‑1 (MQ=255)                                                                                                                                                                                                                             gAAATGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACgg                                                                                                           >  1:248592/1‑163 (MQ=255)                                                                                                                                                                                                                               aaTGGTGATAATTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaa                  >  2:245729/1‑250 (MQ=255)                                                                                                                                                                                                                                  ggTGTTAATTTTCATTATCGTCCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTGCTGCCACTTCCACCAGCCGACTTTGACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCTGCACCGCGCGGGCATGCGCGGTTTGCGGATGTCCCTGACCGATGACGTCCTGGTCTTCCTGTGAGaaaaa               <  2:373776/250‑1 (MQ=255)                                                                                                                                                                                                                                     gATAAGTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCTCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTgg           <  1:403104/251‑1 (MQ=255)                                                                                                                                                                                                                                      aTAATTTTCATTAGCGACCCGGCCGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg          <  2:89292/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaTTTTCATTAGCGACCAGGCAGGCTGGGAAATCGCCCGGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGTCTTCCTGGGATAAAAAGTGGGta        <  2:450315/251‑1 (MQ=255)                                                                                                                                                                                                                                              cATTAGCGACCAGGCAGGCTGGGAACTCGCCCTTCTCCGTCGTCCCTTTTGCCCCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGCCACACTCAAGCTTCGGCCCACCGTCAAAGATGTCGATATAGTTTCGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATc  <  2:241610/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       TGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGCAGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATC  >  NC_000913/1829646‑1830133

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