Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,954,628	C→T	G794S (GGC→AGC)	torZ ←	trimethylamine N‑oxide reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,954,628	0	C	T	100.0%	117.3 / NA	39	G794S (GGC→AGC)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/22);  total (17/22)

AAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCA  >  NC_000913/1954381‑1954872                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     aaaCGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACCCCTTATCCCGGCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTTCATTGCTCGt                                                                                                                                                                                                                                                   >  1:485313/1‑251 (MQ=255)     gcggggTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTatat                                                                                                                                                                                                                                               <  2:217262/247‑1 (MQ=255)              tgattgcgggggaGCCGTAAAGGTCTGGTATTACGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAACCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑GATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGat                                                                                                                                                                                                                                       <  1:395132/238‑1 (MQ=38)                      cgcAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCaa                                                                                                                                                                                                                              <  2:334994/251‑1 (MQ=255)                      cganaGCCGTAAAGGTCTGTTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCaa                                                                                                                                                                                                                              <  2:600339/247‑1 (MQ=255)                       gcAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAAc                                                                                                                                                                                                                             >  1:6470/1‑251 (MQ=255)                           gCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTg                                                                                                                                                                                                                         <  1:626092/251‑1 (MQ=255)                                                  ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGcc                                                                                                                                                                                                  >  1:393283/1‑251 (MQ=255)                                                      aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACTTCTTACGCCTGAACAGCTGGCTGATCAAACGCCCTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCt                                                                                                                                                                                              <  1:306378/251‑1 (MQ=255)                                                              tGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTg                                                                                                                                                                                      >  1:3704/1‑251 (MQ=255)                                                                  cGCCGTATCCGGCCTACAAAAACTTGCTAATTCAAGATTTTGCAGGGCCTATGTAGGCCTGATAATCATAGCGCATCAGGCAGCTTTAGGATTG‑CATAACCTCAGCGCCCGCTTCC‑GGCCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGAGCATTCCTCGTCTATTTTTCGATCTACACCAACGCTGAGTTACCGCCACAGGCATTTGCCAGCTGCGATCTgggg                                                                                                                                                                                   >  2:530644/1‑250 (MQ=255)                                                                        aTCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATcc                                                                                                                                                                             <  2:533618/250‑1 (MQ=255)                                                                                  aaaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAacac                                                                                                                                                                  <  1:244775/250‑1 (MQ=255)                                                                                          tGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCGGCGAGCTGGGGATATCCGCCGTTAACACGTTCGc                                                                                                                                                            <  1:569631/249‑1 (MQ=255)                                                                                                      attttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAGCCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCACCGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTa                                                                                                                                               <  1:213452/247‑1 (MQ=255)                                                                                                                   tatctacgCCTAATAAGCATTGCGCATCAGCCAGCTTTACGTTTG‑CATAACCCCAGCGCCCGTTTCC‑GGGTGCTATTCCGGTTTTACCCCTAATCAGCTGGCTGACTACACGCCATTAATTTCGGTGCTTCGCTCGTATATTTTTCGATCTACACCAACGCTGAGTTACCGTCACAGGCATTTGCCAGCTGCTAGCTGGGGATACCCGCCGTTAACTCGTTCGCACTGCCGTTTTTACACAAGCCAttt                                                                                                                                   <  1:329255/242‑1 (MQ=255)                                                                                                                    aTGTAGGCATGATAAGAATAGCACATCAGGCAGCTTAACGTCTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGACCAGCTGGCTGATCAAACACCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTAGCACTGCCGTTCTTACCCAAGCCAATTTcc                                                                                                                                >  2:450693/1‑251 (MQ=255)                                                                                                                        aGGCGTGATAAGCCTGGCGTATCAGGCATCTGTACTTTTG‑CATACCCCCAGTGCCCGTTTCC‑GGGCGCTTTTCACGTCGTCCCCCTAAACACCAGGCAGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATTTTTTCCGATATCCACCAACGCTGAGTGTCCGGCTCAGCCATTCGCGTGTTGCGATCTGGGAATATCCGCTGTTAACGCGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGAt                                                                                                                            <  2:79378/251‑1 (MQ=255)                                                                                                                        aGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTcc                                                                                                                                >  2:164851/1‑247 (MQ=255)                                                                                                                                  aGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGc                                                                                                                  >  2:579321/1‑251 (MQ=255)                                                                                                                                  aGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGc                                                                                                                  >  1:465099/1‑251 (MQ=255)                                                                                                                                  aGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGAGCCAGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATCTACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGGTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGc                                                                                                                  >  2:121390/1‑251 (MQ=255)                                                                                                                                   gCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCa                                                                                                                 <  1:644180/251‑1 (MQ=255)                                                                                                                                     agtatcgctgctgcgtgatttcGTTTGTCATACCCTCGGCGCC‑GTTTCCCGGGCGCTATTCCCTCCTTCCGCCTGAACAGCGGGCTGACCAAACCCCGTTAACTTCGGTGTTTTCCTCGTATATTTTTCGTTATACACCCACCATGATCTCCCGGCACAGGCATTTGCAAGCGTCGAGCTGGGCATATCCGCCGTTAACACGTTCGCTCTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCAc                                                                                                                <  2:319491/230‑1 (MQ=255)                                                                                                                                       aGCGCATCAGGCAGCTTTCCGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCtt                                                                                                             <  1:611603/251‑1 (MQ=255)                                                                                                                                        gcgcATCAGGTATCTTTACGTTTG‑CATACCCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTc                                                                                                            <  1:445013/251‑1 (MQ=255)                                                                                                                                              cAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTCCGCCTGAACAGCTGGCTGACCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCAc                                                                                                      <  1:167729/251‑1 (MQ=255)                                                                                                                                                                      ccTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATcc                                                                               <  1:111345/251‑1 (MQ=255)                                                                                                                                                                                    tCC‑GGGCGCTATTCACGTCTTAGGCCTGAACAGCTGGCTGATCAAACGCCGTTACCTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCGGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTcaccac                                                                  <  2:6470/250‑1 (MQ=255)                                                                                                                                                                                     cc‑GGGCGCTATTCACTTCTTACGCTTGAACAGCTGGTTAATCAAACGCCTTTAACTTCGGTGCATTGCTAGTATATTTGTCGAAATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATAGCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCGCGCATACCACGCCTTTTTGGATCCCGTCAGTCACCACcg                                                                <  1:561854/251‑1 (MQ=255)                                                                                                                                                                                                   aCGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCGGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAAt                                                   >  2:323958/1‑251 (MQ=255)                                                                                                                                                                                                     gTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACAGCGCCTGTCAAAATCt                                                 >  1:79912/1‑251 (MQ=255)                                                                                                                                                                                                                        ggctgATCAAACGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATCTACACCAACGCTGAGATACCGGACCAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGac                              >  1:605746/1‑251 (MQ=255)                                                                                                                                                                                                                             ataaaacccgtttaacgtgggTGAATTGCTCGTATATTTTTCGATATACACCACCGCTGAGTGACCGGCACAGGCTTTGGCCAGCTGCGAGCTGGGGATTTCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGCCCACGTTTGTTCCAGACACGCa                         <  2:559024/232‑1 (MQ=255)                                                                                                                                                                                                                                aaaCGCCGTTAACTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCAcc                       >  1:144848/1‑250 (MQ=255)                                                                                                                                                                                                                                           aCTTCGGTGCATTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTc           >  1:299541/1‑251 (MQ=255)                                                                                                                                                                                                                                               cGGTGCATTGCTCGTATATTTTTCGATATCCCCCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCGGGGGATTTCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCa        <  2:465099/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gTGCATTGCTCGTATATTTTTCGATATACACCCACGCTGAGTTAACGGCACAGGCATCTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAAt      >  1:631637/1‑250 (MQ=255)                                                                                                                                                                                                                                                     aTTGCTCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCa  >  1:569772/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     AAACGCGCGGTGCTGATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTG‑CATAACCTCAGCGCCCGTTTCC‑GGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTTTACACAAGCCATTTTCCAGATCTGGCCATGCACCTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCA  >  NC_000913/1954381‑1954872

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: ‑