Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,831,664	C→T	E107K (GAA→AAA)	astC ←	succinylornithine transaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,831,664	0	C	T	97.5%	129.9 / ‑4.1	41	E107K (GAA→AAA)	astC	succinylornithine transaminase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/21);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTCCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACC  >  NC_000913/1831422‑1831908                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     cACCGCACAGGTAGAGTCGTCATTCAGCGCGCTGGCAGAGTTAATAGCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCgg                                                                                                                                                                                                                                              <  2:328523/251‑1 (MQ=255)     gCACAGGTAGAGTCGTCACTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCCGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGtt                                                                                                                                                                                                                                          >  2:171117/1‑251 (MQ=255)                     cAATCAGCGCCCTGGCAGAGTGAATATCGTTAGTTGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTTTTAATAGGCTGTCTGCCCCCCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTCCGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGAcgcca                                                                                                                                                                                                                          <  2:219286/251‑3 (MQ=255)                       aTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATc                                                                                                                                                                                                                        <  2:262286/251‑1 (MQ=255)                                     agagTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGCAGCGGTCGTGCGCGAATTTACGCGCCAGTTTTCGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACATAAGAAGACGCGATCGGCAAACGAGGCGt                                                                                                                                                                                                          >  2:346076/1‑251 (MQ=255)                                              atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCACAAAATCCTGTGACTAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGACGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCCGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAGAAAGACGCGCTCGGCAAACGTAGCGTCGATCAAtt                                                                                                                                                                                                 >  2:182993/1‑251 (MQ=255)                                                       tGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGAGACCATGAAACGAATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGGGAGCGAATTTACGCGCCAGATTTAGCGCCGCTACGGTGGCTTACGCACCGGAGGTACAAAAGAAGACGCGATCGGCAAACGGGGCGACGAGAAATGTTTTCGca                                                                                                                                                                                         >  1:145236/1‑249 (MQ=255)                                                                        ccgccgGCATTGGCGCAAAATCCTGTGAATATGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCt                                                                                                                                                                         <  2:361264/249‑1 (MQ=255)                                                                        ccgccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCGTGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCg                                                                                                                                                                       <  1:160877/251‑1 (MQ=255)                                                                        ccgccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCTCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCGTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCg                                                                                                                                                                       <  2:236518/251‑1 (MQ=255)                                                                         cgccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTcct                                                                                                                                                                      >  2:260446/1‑249 (MQ=255)                                                                         cgccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCgg                                                                                                                                                                      >  1:95169/1‑250 (MQ=255)                                                                          gccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGtt                                                                                                                                                                     >  2:149842/1‑251 (MQ=255)                                                                          gccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGtt                                                                                                                                                                     >  1:161655/1‑251 (MQ=255)                                                                           ccgGCAGTGGCGCAAAATCCTGTGAAGAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCGGCCGTAGCGGTCGTGAGCGAAGTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCTATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTg                                                                                                                                                                    <  2:31104/251‑1 (MQ=255)                                                                             ggCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGgt                                                                                                                                                                  <  2:19631/251‑1 (MQ=255)                                                                               cAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGgt                                                                                                                                                                  >  1:121838/1‑249 (MQ=255)                                                                               cAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGgt                                                                                                                                                                  <  2:121838/249‑1 (MQ=255)                                                                                             ccTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGt                                                                                                                                                   <  2:334159/250‑1 (MQ=255)                                                                                                       ggctggatGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCgggattcccgacct                                                                                                                                        >  1:71908/1‑239 (MQ=255)                                                                                                          tggctgCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGGATGCCAGAACTTAc                                                                                                                                     >  1:92579/1‑251 (MQ=255)                                                                                                                      cGCACTGACAGGAAACAGCGTGCGACCATGAAACGCATTTGTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTCCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTATCCGTTGCCGGTATGCCAGAACTGACTCGCCTgttcgt                                                                                                                         <  2:58244/251‑1 (MQ=255)                                                                                                                           tGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAgc                                                                                                                    >  1:362783/1‑251 (MQ=255)                                                                                                                                gTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTc                                                                                                               <  2:92579/251‑1 (MQ=255)                                                                                                                                 aaaaCAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCa                                                                                                              <  2:228348/250‑1 (MQ=255)                                                                                                                                  aaaCAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCAc                                                                                                             >  1:193273/1‑251 (MQ=255)                                                                                                                                         gTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTATTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGGCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAAACGTTGCCGGTATGCCCGAACTTACTCGCCTGTTCGTTTAGCGCTTCACGCAgctc                                                                                                      >  2:201721/1‑248 (MQ=255)                                                                                                                                           gCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCATAACTTTCTCTCCTGTTCGTTCAGCGCTTCACGCAGTTcc                                                                                                     <  1:314653/250‑1 (MQ=255)                                                                                                                                                catgcaaagCATTTTTGAACTCCTCCAGGCTGCTCTTATGGCTGCCGTAGCGGTATTGAGTGATTCTACGCGCCAGTTCTAGCGCAGCTTCGTTGGCTTTCGCACTGGACTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTCTTCGCCAGTCTTAGTCCCGGCTCGTTGGTGTACACGTTGCCTGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTGCCCGCAGTTACGGATgc                                                                                               <  2:387125/243‑1 (MQ=255)                                                                                                                                                          aTTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCTTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCa                                                                                     <  2:42979/251‑1 (MQ=255)                                                                                                                                                                 aaCGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCTATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGtt                                                                              <  1:8037/251‑1 (MQ=255)                                                                                                                                                                     ccACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCAcc                                                                          <  1:240762/251‑1 (MQ=255)                                                                                                                                                                           ggCCGCTCTTATGGCTGCCGTAGCGGTCGTGCGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTAGGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAAt                                                                     <  1:182993/249‑1 (MQ=255)                                                                                                                                                                            gCCGCTCTTATGGCGGCCGTACCGGTCGCGAGCGAATTTACACGCCAGGTTTAGCGCCGCTTCGTTGTCTTTCGCACCGGAGTTACAAAAGGAGACGCGATCGGGAAACGTGGCGTAGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTGGCCGATATGCCAGATCTTGCTCGCCTGTTCGTTCAGCCCGTCACGCAGTTCCGGAGGCGTGTGGCCCACCGCGTTAACCGCactg                                                                    <  2:85898/250‑4 (MQ=255)                                                                                                                                                                                    tATGGCTGCCGTAGCGGTCTTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCcg                                                             <  2:190751/249‑1 (MQ=255)                                                                                                                                                                                          tGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGGTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGGTCACCGCAATGCCACCCGCGAAggcg                                                     >  2:117551/1‑248 (MQ=255)                                                                                                                                                                                            ccGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGat                                                   >  2:187706/1‑251 (MQ=255)                                                                                                                                                                                                aGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGCtata                                                >  2:358620/1‑250 (MQ=255)                                                                                                                                                                                                           gCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCctg                                    >  1:269283/1‑251 (MQ=255)                                                                                                                                                                                                                                   gccgcGTCGTTGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCACACAAGCGCGAACCTTc            <  1:201721/251‑1 (MQ=255)                                                                                                                                                                                                                                             tGGCTTTCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGGCGATATACTCTTTCCCCTGCTGATCCCACAACCGCGAACCTTCGCCACGTAcc  >  2:215849/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTCCGCACCGGAGTTACAAAAGAAGACGCGATCGGCAAACGTGGCGTCGATCAATTTTTTCGCCAGTCGCAGTACCGGCTCGTTGGTGTAACCGTTGCCGGTATGCCAGAACTTACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACC  >  NC_000913/1831422‑1831908

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