Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,191,644	C→T	E115K (GAA→AAA)	purB ←	adenylosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,191,644	0	C	T	97.5%	129.7 / ‑3.0	40	coding (343/1371 nt)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/21);  minor base . (1/0);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGGTACTGGCGCTCCATACGGTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACTTGTACA  >  NC_000913/1191400‑1191879                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            cGGTACTGGCGCTCCATACGGTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAAc                                                                                                                                                                                                                                       <  2:226750/251‑1 (MQ=255)             tCCATACGGTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCCGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGCTGAATTTAGAAACCGCGTGCAGTTc                                                                                                                                                                                                                           >  1:191750/1‑251 (MQ=255)              ccATACGGTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTcc                                                                                                                                                                                                                          <  2:270258/251‑1 (MQ=255)                  aCGGTAGGCGACGTTTGCCCTCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGCAGTACAGGCAGAGTGGATGAATTTAGAAACCGCGTGCAGTTCCgggg                                                                                                                                                                                                                      >  2:418516/1‑250 (MQ=255)                  aCGGTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGa                                                                                                                                                                                                                      <  2:549390/251‑1 (MQ=255)                     gTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGAtct                                                                                                                                                                                                                   <  2:464648/251‑1 (MQ=255)                                ttGCCATCTCTTTACCGATGGTTGACGGCGTGGTTGGCTGACCGTGGGTACGAGACAGCATCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCGCGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTAGAGGCAAAGTGGATGAAGTTAGAAACCGCGTTCAGTTCCGGGATCGCCGCCACttt                                                                                                                                                                                                         <  2:228810/250‑1 (MQ=255)                                ttGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGGGCAGTTCCGGGATCTCCGCCACtttt                                                                                                                                                                                                        >  2:179003/1‑251 (MQ=255)                                          tttACAGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGGTATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGGTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAgg                                                                                                                                                                                               >  2:84538/1‑250 (MQ=255)                                               cgttGGTTGACGGCTTGGCTGGCTGACCGTGGGTACGAGACACCAGCGGGATAGCGCGATACAGAACGGCGAGATCTTTAATGCCATCAATCATTTGACGCCAGTATGGCAGGATCACCTCATCACGCGCGGCTGTCAGCATTAATGCGTGGGAGAGTTTTTTTATATCTTCAGAAGTACATGCAAAGTGGATGAATTTAGAAACAGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATAc                                                                                                                                                                                         <  2:397686/248‑1 (MQ=255)                                                   ggTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCa                                                                                                                                                                                      >  1:156785/1‑250 (MQ=255)                                                           gCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAAGCAGTTGCCGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCTAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGGTCCGGGATCTCCGCCGCTTTTTCTTTCAGGAAATACTCAACCGCTTTa                                                                                                                                                                             <  2:572595/251‑1 (MQ=255)                                                                               gTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATT‑AGAAACCGCGTGGAGTTCCGGGATCTCCGCCAATTTTTCTTTCAGGAAATAATCAACCGCTTTAACGTCGTGGTTAGTGGGACGc                                                                                                                                                        >  2:54797/1‑251 (MQ=255)                                                                                           agtgggATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGtttt                                                                                                                                             <  2:109976/248‑1 (MQ=255)                                                                                              gggATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGAt                                                                                                                                          <  1:143163/251‑1 (MQ=255)                                                                                                                   gCGAGATCTTTAATGCCAGCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGGGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCt                                                                                                                     <  2:537139/251‑1 (MQ=255)                                                                                                                      agaTCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGa                                                                                                                   >  1:145486/1‑250 (MQ=255)                                                                                                                           tttAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCAtcttc                                                                                                                           >  1:523832/1‑237 (MQ=255)                                                                                                                           tttAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGGGGGAGAGTTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTCCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGGTTTGATGCGCGCCGCAtcttc                                                                                                                           <  2:523832/237‑1 (MQ=255)                                                                                                                                   cATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCTCTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGAtt                                                                                                     >  1:28115/1‑251 (MQ=255)                                                                                                                                   cATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGAtt                                                                                                     >  1:379733/1‑251 (MQ=255)                                                                                                                                   cATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGAtt                                                                                                     >  1:106380/1‑251 (MQ=255)                                                                                                                                                 cGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAAcc                                                                                       >  1:386018/1‑251 (MQ=255)                                                                                                                                                     aGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGAtt                                                                                   <  2:379733/251‑1 (MQ=255)                                                                                                                                                        atgggagtaTCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAATGCCAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAAGCGATTGCg                                                                                <  1:467422/243‑1 (MQ=255)                                                                                                                                                             aGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGCTTGCATCAAGGTAACCGATTGCGTCGgc                                                                           >  2:171644/1‑251 (MQ=255)                                                                                                                                                                                 ggTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCt                                                       >  2:341147/1‑251 (MQ=255)                                                                                                                                                                                  gTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCtt                                                      <  2:20491/251‑1 (MQ=255)                                                                                                                                                                                         gCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAATCAGGAACTTCCTTGATCGCt                                               <  1:371355/251‑1 (MQ=255)                                                                                                                                                                                         gCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCt                                               >  1:114312/1‑251 (MQ=255)                                                                                                                                                                                           aTTAATGCGTGGGAGAGGTTATTGATATCTCCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTGTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCtgc                                             <  2:45642/251‑1 (MQ=255)                                                                                                                                                                                               aTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCtgcgtgc                                         <  1:84538/251‑1 (MQ=255)                                                                                                                                                                                                                   atatCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGcc                      >  1:54568/1‑250 (MQ=255)                                                                                                                                                                                                                    tatattCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGAAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGcaa                     >  2:347937/5‑248 (MQ=255)                                                                                                                                                                                                                      tCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGcgacg                                                                                                        >  2:27985/1‑165 (MQ=255)                                                                                                                                                                                                                      tCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGcgacg                                                                                                        <  1:27985/165‑1 (MQ=255)                                                                                                                                                                                                                            gAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTc            >  2:14977/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCAGGCAAAGGGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCGCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACtt       <  1:418516/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACTTGTa    <  2:156785/251‑1 (MQ=255)                                                                                                                                                                                                                                      cAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACTTGTaca  <  1:179003/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            CGGTACTGGCGCTCCATACGGTAGGCGACGTTTGCCATCTCTTTACCGATGGTTGACGGCGTGGCTGGCTGACCGTGGGTACGAGACAGCAGCGGGATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACTTGTACA  >  NC_000913/1191400‑1191879

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