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.8%	140.1 / NA	45	A293A (GCC→GCT)	astA	arginine N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/26);  minor base T (1/0);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTG  >  NC_000913/1829652‑1830123                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              aaCACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGAATCGGTACGCACCAGACCCACGCGGAAATGGTGGTAATTTTCATTAGCGACCAGGC‑AGGCt                                                                                                                                                                                                                               >  1:83723/1‑251 (MQ=255) aaCACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCGCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCt                                                                                                                                                                                                                               <  1:326535/251‑1 (MQ=255)           ccGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCCGTCACCGTTACTCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGCACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCACACGCTCGTTTGCCGGGTCGGTACCCACCCGCACCACGCGGACATTGGGATACTTTTCTTTTGCCACCAGGCCAGCCTGGGaatcgc                                                                                                                                                                                                                      >  2:97294/1‑247 (MQ=255)            cGATACCGGATTACGCTTCACACGCGATGCGCACTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGGGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGcc                                                                                                                                                                                                                    >  1:403635/1‑251 (MQ=255)                             tCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGccc                                                                                                                                                                                                     <  2:269760/249‑1 (MQ=255)                                 acGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCACGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGCTAAAATCAAACGCTCGGTTGCCGGCACGGTACGCACCAGCACCACGCGGAAATGGGGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGGAATCGCCCTGCGCAGGCTGCCCTTCTGc                                                                                                                                                                                               >  1:230057/1‑251 (MQ=255)                                     gATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACt                                                                                                                                                                                           >  1:360114/1‑251 (MQ=255)                                                       aTCCAGTCACCGTTAATCCATAAAGTCATGCTGTGTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCTATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTa                                                                                                                                                                         <  1:501177/251‑1 (MQ=255)                                                            gTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGGATCGGGACGCACCAGCACCACGCGGAAATGGTGATAATTTTTATTAGCGACCAGGC‑AGGCTGGGGAATCGCCCTGCGCAGGCTGCCCTTCTGCCCCTTCCCCCAGCCGACTTTTACGGAt                                                                                                                                                                    >  2:148342/1‑251 (MQ=255)                                                                             aaGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCCGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTTCGGATGGCGCGCACGCgggcgg                                                                                                                                                   >  1:36229/1‑247 (MQ=255)                                                                                 cATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTc                                                                                                                                               <  1:241827/251‑1 (MQ=255)                                                                                    gCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATAAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCCCGCGGAAATGGGGATAATTTTAATTAGCGACCAGGC‑AGGCTGGGAAATAGCCCTGCGCAGGGTGCCCTTCTGCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTcac                                                                                                                                             >  2:340037/1‑250 (MQ=255)                                                                                          gtcgcctcTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCGCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTTCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCaa                                                                                                                                       <  2:139253/246‑1 (MQ=255)                                                                                              cctcTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCGCTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCg                                                                                                                                    <  2:34652/249‑1 (MQ=255)                                                                                               ctctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCCCCAGCCGACTTTTACg                                                                                                                                                                       >  1:126298/1‑213 (MQ=255)                                                                                               ctctGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACg                                                                                                                                                                       <  2:126298/213‑1 (MQ=255)                                                                                                              cACCAGCCGATCGCCATCCCCGCCGTGCCATTTGGGGGCATCCAGTGCTGCGGCGGTGAAAATCAAACGCCCGGTTGCCGGATCGGTACGCCCGAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGGCCTGCGCAGGCTGCCCTTCCGCCACTTCCACCAGCCGACGTTTGCGGATGGCGCGCACGCGGTCGATGACACACTCTAGCGTCGGCCCACCCCCAAAGa                                                                                                                  <  2:191818/251‑1 (MQ=255)                                                                                                                  aGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCGGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCGAGCGTCGGCCCACCGTCAAAGATGTc                                                                                                              <  1:415657/251‑1 (MQ=255)                                                                                                                               ccccGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCTGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGTCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACgg                                                                                                 <  1:237529/251‑1 (MQ=255)                                                                                                                                               gAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                                  <  1:9884/250‑1 (MQ=255)                                                                                                                                                aGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGGTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCCCCCTGCACAGGCCGCCCTTCTGACACTTCCACCAACCGGCTTTTACGGATTGCGCGCACCCGGTCTATGTCACACTCAAGCGTCGGCCCACCGTCAAAGGTGTCCATTTAGTTAAGGTAGCGAAAACCttcttt                                                                                >  1:376634/1‑251 (MQ=255)                                                                                                                                                 gggCATCCAGTTGTGCTGAGGTTAAAATCGAACGCTCGGTATCCGGAGCGGTACGCACCAGCACCCAGCGGAAATGGTGATAACTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCGCTGCGCATGATCCCCTTCTGCCACTCCCCCCTGCCGACCTTTGCGGATGGCGCGCACGCGGTCGATGTCACACTCAATCGTCGGCCCGCCAGCAAAGATGTGGATATAGTTGCGGTATGTAAAACCttctttc                                                                               <  2:75917/251‑1 (MQ=255)                                                                                                                                                     aTCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCCGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCCAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTcc                                                                            >  2:252041/1‑250 (MQ=255)                                                                                                                                                            tgtgCGGCGGTTAAAACCAAACGCTCGGTTGCCGGATCGGTACGCCCAAGCCCCACGCGGAAATGGTGATAATTTTCATTAGCGCCCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCACGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACcg                                                                     <  2:378897/250‑1 (MQ=255)                                                                                                                                                                cggcggTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCgg                                                                <  2:378650/251‑1 (MQ=255)                                                                                                                                                                    ggTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCTCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATATAGTTa                                                                                                    <  1:369598/211‑1 (MQ=255)                                                                                                                                                                    ggTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTa                                                                                                    >  2:369598/1‑211 (MQ=255)                                                                                                                                                                     gTTAAAATCAAAAGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCAACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTAAAGCGTCGGCCCACCGTCAAAGCTGTCGCTATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCCCCGCGCGGGCAgg                                                           >  1:421193/1‑251 (MQ=255)                                                                                                                                                                       tagaagcatacttgctgtttccGGGTCGGTACGCCCCAGCACCACGCGGAAATGGTGATAGTTTTAATTAGCGACCAGGC‑AGGCTGGGGAATCGCCCTGCGCAGGCTTCCCTTCTGCCACTTCTACCAGCCGTCTTTTACGGATGGCGCGCACGCGGTCGATGTCACACGCGAGAGTCGGCCCACCGTCAAAGATGGCGATATAGTGACGGTCGCGGAATCCTCCTGTCTTCAGCACCGCGCGGGCAGgc                                                          <  2:455783/230‑1 (MQ=255)                                                                                                                                                                       tAAAATCAAACGCGCGGTGGCCGGATCGGTACGCACCAGCACCCCGCGGAAATGGTGATAATTTTCATTAGCGACCATGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTGCGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTATCGAAAACCTTCTTTCTCCAGCTCCGCGCTGGCAgtc                                                          <  2:316349/250‑3 (MQ=255)                                                                                                                                                                          aaTCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGCTGATAATTTTCATTAGCGACCAGAC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGACCTTCTGCCACTTACACCAGCCGACTTTTACGGATGGCGCGCAAGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAACGATGTCGATATAGTTACGGTATCGAACACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCgg                                                      >  1:249333/1‑251 (MQ=255)                                                                                                                                                                                      ggTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                           >  1:17650/1‑250 (MQ=255)                                                                                                                                                                                           ccGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTgg                                     >  2:454740/1‑250 (MQ=255)                                                                                                                                                                                                   gTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGAc                             <  1:150601/251‑1 (MQ=255)                                                                                                                                                                                                       gcaccagcaacaCGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTGTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGGCAAAGATGGCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCGGCGCCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTcc                         <  1:317996/251‑1 (MQ=255)                                                                                                                                                                                                         accagcaccaCGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTg                       >  1:127387/1‑251 (MQ=255)                                                                                                                                                                                                          ccagcaccaCGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGTGCAGGCGCGGTTTGCGGATGTACCTGACCGAACACGTCCTgg                      >  1:47314/1‑251 (MQ=255)                                                                                                                                                                                                             gcaccaCGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCt                   <  1:252041/251‑1 (MQ=255)                                                                                                                                                                                                                 cacGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAATCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCt               <  1:271436/251‑1 (MQ=255)                                                                                                                                                                                                                 cacGCGGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCGGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCt               <  2:452360/251‑1 (MQ=255)                                                                                                                                                                                                                   cgcgGAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTgg             <  1:454740/251‑1 (MQ=255)                                                                                                                                                                                                                      ggAAATGGTGATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAt          <  1:184129/251‑1 (MQ=255)                                                                                                                                                                                                                        aaaTGGTGATAATTTTCATTAGCGACCAGGC‑GGGCTGGGAAATCGCCCTGCGCAGTCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGTAGACGTCCTGGGCTTCCTGGGATa         >  2:84286/1‑250 (MQ=255)                                                                                                                                                                                                                             gggagaaTTTTCATTAGCGCCCGTGC‑ATGCTGGGAAATCGCCCTGCGCATGCTTCCCTTCTTCCTCTTCCACCAGCCGACTTTTACGTATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTATCGAAACCCTTCTTGCTCCAGCGCCGCGCGGGCAGGCGCGGTGTGCGGATTTACCTGACCGAGGACGTCTTGGGTTTCCTTGGATAAAAAGt   <  2:443508/246‑1 (MQ=255)                                                                                                                                                                                                                               gATAATTTTCATTAGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTg  <  2:17650/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              AACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTTGAGGGCATCCAGTTGTGCGGCGGTTAAAATCAAACGCTCGGTTGCCGGATCGGTACGCACCAGCACCACGCGGAAATGGTGATAATTTTCATTGGCGACCAGGC‑AGGCTGGGAAATCGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTG  >  NC_000913/1829652‑1830123

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