Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,152,494	G→T	A87S (GCG→TCG)	metC →	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,152,494	0	G	T	88.9%	80.0 / ‑4.2	38	A87S (GCG→TCG)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (28/4);  total (30/8)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCG  >  NC_000913/3152273‑3152739                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        cAGGACGCAGCAAACAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCt                                                                                                                                                                                                                          >  1:437957/1‑251 (MQ=255) cAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTTAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGAAGCG‑GTTGCT‑AATTCCATTCTTGCt                                                                                                                                                                                                                          >  1:50451/1‑251 (MQ=255) cAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAACCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCt                                                                                                                                                                                                                          >  1:111644/1‑251 (MQ=255)        cagcaAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTTCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCg                                                                                                                                                                                                                   <  1:272916/251‑1 (MQ=255)           caAAAAATACACTCTCGGCGCGGCAAATAGCGTGGTTCAGCGCGAATCTTCGCTGATCTTTGACAGTGTAGAAGCCAAAGAACACGCGACACGTAATCGCGCCAATGGGGAGTGGGTCTATGGACGGCGCGGCACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGACTTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCCTTCTTGCTTTTATCGAAc                                                                                                                                                                                                                >  2:528455/1‑251 (MQ=255)                 aTACACGCGCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGGGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTTTTGCTTTGATCGAACAGGGCg                                                                                                                                                                                                          <  2:507530/251‑1 (MQ=255)                      ctctCGACGCGGTAAATAGCGTGATTCAGCACGCTTCTTCGCTGGTCTTTGCCAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGAGCGGAACGTTAACCCATTTCTCCTTACAACCAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCCTTCTTGCTTTTATCGAACAGGGCGATCAt                                                                                                                                                                                                     >  2:478885/1‑251 (MQ=255)                              gcgGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgat                                                                                                                                                                                             >  2:420021/1‑251 (MQ=255)                                ggTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatg                                                                                                                                                                                            >  1:516458/1‑250 (MQ=255)                                     aTAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAac                                                                                                                                                                                      >  1:433553/1‑251 (MQ=255)                                     aTAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAac                                                                                                                                                                                      <  1:434439/251‑1 (MQ=255)                                      tAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAac                                                                                                                                                                                      >  1:459061/1‑250 (MQ=255)                                      tAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAacc                                                                                                                                                                                     >  2:272387/1‑250 (MQ=255)                                           tGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACATGGCGATCATGTGTTGATGACCAACACCGcc                                                                                                                                                                                >  1:387013/1‑251 (MQ=255)                                                    gcgcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCAGTCATGTGTTTATGAACAACACCGCCTATGAAc                                                                                                                                                                         >  2:89574/1‑249 (MQ=255)                                                         tcttcGCTGGTCGTTGACAGTATAGAAGCCCAAAAACACGCGACACGTAATCGCGCCCATGGAGAGTTGTTCTATGGAAGGCGAGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGG‑GTGAACTGGCGGGAGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GGTGAT‑AATTCCATTCTTGCTTTTATCGAACAGTGCGATCTTGTGTTGTTGACAAACACCGGCTTTGAACCgcgac                                                                                                                                                                   >  2:387371/1‑246 (MQ=255)                                                                     tttGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGt                                                                                                                                                        >  1:175821/1‑249 (MQ=255)                                                                         aCAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATTACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCaa                                                                                                                                                   >  1:382220/1‑250 (MQ=255)                                                                         aCAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTATATGGAAGGCGCGGAACGTTAAACCATATCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGAAGGCTGCGTGCTAATTCACTGCGGGTCGGCAGCG‑GTTGAT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATTACCACAAACGCATATGAACAGAGTAAGGATTTCTGTAggaa                                                                                                                                                   >  1:382244/1‑247 (MQ=255)                                                                             cactcacgatcaaaaaCACGCACCACTTAGTCGCGCAAAGGGGGATTTTTTCTAAGGCCGGCGCGGACTTTTAACCAATTGCCCTTTACAACAAGCGTGGT‑GTGATCTGGAAGGGGGCGCAGGCTGCGTATTATTTCCCGGGGGGTCGTCGGCG‑GTTGAT‑ATTTCCATTCTGGTTTTTATCAAAGAGGGCTTTCTTGTTTTAATGACAAACGCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATc                                                                                                                                               <  1:523734/239‑1 (MQ=255)                                                                              ggaaaagtaaaaaaaCACGCGAACCGTAAACGCTCTAAGGGAGAGTTTTTCTATGGACGGGGCGTAACTGTAACCCATTTCCTCTTACAACAAGCGATGT‑GTGAACGGGAAGGTGGCGCCGGTTGTGTGTTATTCCCCTG‑GGGGGGGC‑GTGTGGTGTTGAATTCCATCCTTTCTTTTATCAAAAGGGGCTATCTTGTGTTGTTGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                                                                                                             <  1:528455/242‑1 (MQ=255)                                                                                  aaGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                                                                                                         >  2:144356/1‑251 (MQ=255)                                                                                     ccAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTAACAACAAGCGATGT‑GTGAACTGGACGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTAGGCAGCG‑GTTGCT‑CATTCCATTCTTGCTTTTATCGAACAGGGCGCTCATGTATTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTCGCAAACTCCTCAGCaa                                                                                                                                       >  2:436346/1‑250 (MQ=255)                                                                                         aaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTTTAGCAAAATCCTCAGCAAACTgg                                                                                                                                  >  2:234047/1‑251 (MQ=255)                                                                                         aaaaCAAGCGACACGTAATCGCGCCAATGGCGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCAAGTCAGGATTTCTGTAACAAAATc                                                                                                                                               >  2:523734/1‑238 (MQ=255)                                                                                                gcgACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCATGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                                                                                                           >  1:491933/1‑251 (MQ=255)                                                                                                        taagcgcgCCAATGGAGAGCTGGTCTATGGACCGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATTTGTTGATGACCAACACCGCCTATGAACCAAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGCCGTAACGAcatcat                                                                                                                   >  2:47519/5‑251 (MQ=255)                                                                                                                 cAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGCTGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTGCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGATTGATc                                                                                                          >  1:342125/1‑251 (MQ=255)                                                                                                                   aTGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCg                                                                                                        >  2:185493/1‑251 (MQ=255)                                                                                                                                  aTGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGat                                                                                         >  1:543226/1‑251 (MQ=255)                                                                                                                                                      aCCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGCAGGCGCCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCATATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGCAATCATGGTTTGATCCGCTGATTGGTGCCGCTATCGTTAAGCATCTGCAGcc                                                                     >  2:538670/1‑251 (MQ=255)                                                                                                                                                               cccTTACAACAAGCG‑TGTTGGGAGCTGGAAGGGGGCGCAGGCTGCGTGTTTTTCCCCGGCGGGCCGGCTGCG‑GTGGTT‑ATTCCCATTCTTTCTTTTATCAAACGGGGCGATCATGTTGTGAGAACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaa                                                            <  2:25022/250‑1 (MQ=255)                                                                                                                                                                   taaaacaaGCGGTGT‑GTGAGCGGGAAGGGGGCGCGGGCCGCGTGCTTTTTCCCTGCGGGGCGGCGGCG‑TTTGCT‑ATTTCCATTCTTGTTTTTACCAAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaaaa                                                          <  1:185493/246‑1 (MQ=255)                                                                                                                                                                             gATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGACCAGGGCGATCATGTGTTGATGACCAACACCGCCTATGACACGAGTCAGGCTTTCTGTAGCAAAATCCTCAGCAAACTGGTCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGCTATCGTTAAGCATCTGCAGCAAAACCCTAAACTCGTGTTTCTgg                                              >  2:144570/1‑251 (MQ=255)                                                                                                                                                                                        cTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAgg                                    >  2:91145/1‑251 (MQ=255)                                                                                                                                                                                         gggAAGGGGGCGCGGGCTGCGTGCTTGTTCCCTGCGGGCCGGCAGCG‑GTTGCA‑AATTCCATTTTTGCTTTTATCGAACGGGGCGATCTTGTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGc                                   <  1:241016/250‑1 (MQ=255)                                                                                                                                                                                                               gCTATTTCCCTGCGGGGCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGGTCATGTGTTGATGACCAACACCGCCTATGAAACGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTCACGACATCATGGTTTGCTCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCTAGTTTGTGGTATCACCAGGCTCCATACCCCTGGAAGTccccg             >  2:152534/1‑248 (MQ=255)                                                                                                                                                                                                                           cGGGTCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCg  >  1:92974/1‑250 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        CAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCG‑GTTGCT‑AATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCG  >  NC_000913/3152273‑3152739

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