Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,736	C→T	R629R (CGC→CGT)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,736	0	C	T	98.0%	165.9 / NA	49	R629R (CGC→CGT)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/24);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGATGCGCTGGAAGTCGCTTTCGACCGTTATCTGGTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAG  >  NC_000913/3062442‑3063015                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                        ttcagATGCGCTGGAAGTCGCTTTCTACCGTTATCTGGTGCCAAGCCAGTGTGTNACCGTCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCTACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGNTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGc                                                                                                                                                                                                                                                                                    <  2:755813/296‑1 (MQ=255)                        atgtataatagacagTGCCAAGCCAGTGTGTTACCGGCGTGATTGCACAAAGCTNTCATCAGTCTGAGAAATCTGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGAGGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGANGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAg                                                                                                                                                                                                                                                             <  2:769377/286‑1 (MQ=255)                          tacctttATCTGGTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATNATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGNTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTg                                                                                                                                                                                                                                                           <  2:757694/295‑1 (MQ=255)                            ccGTTATCTGGTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTggg                                                                                                                                                                                                                                                         >  1:827267/1‑300 (MQ=255)                              gTTATCTGGTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGgcg                                                                                                                                                                                                                                                       >  1:484182/1‑300 (MQ=255)                                      gTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTTGACGTAGATTTCAGCCCGATGTTCTCTACACCTGAAGAGATCGCCAGTCTGGCAGTAGAAAACGGCGTTCACGTAGTGGGCGCATCCTcc                                                                                                                                                                                                                                               >  2:92580/1‑299 (MQ=255)                                      gTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCa                                                                                                                                                                                                                                               >  1:437783/1‑300 (MQ=255)                                       tGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAgctgtctcttatacac                                                                                                                                                                                                                                              >  1:788473/1‑285 (MQ=255)                                       ggcaaaaccAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACTGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAGGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCAc                                                                                                                                                                                                                                              <  2:274289/293‑1 (MQ=255)                                           aaGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTCAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCAAGTAGTGGGCGCATCCTCACTGGc                                                                                                                                                                                                                                          >  2:496799/1‑300 (MQ=255)                                           aaGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCCCGTAGTGGGCGCATCCTCACTGGc                                                                                                                                                                                                                                          >  1:536038/1‑300 (MQ=255)                                           aaGCCAGTGGGTTACCGGCGTGATGGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGGGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGACCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGc                                                                                                                                                                                                                                          <  2:432547/300‑1 (MQ=255)                                                 gtgtgtTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGNGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCNTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCgg                                                                                                                                                                                                                                    >  2:755767/1‑300 (MQ=255)                                                   gtgtTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAACAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGCTGGGCCAGGATGGACACGATCGCGGCTCGAAAGTGATCGCCAGCGCCTATTCCGATCCCGGTTTCGACGTAGCTTTCCGCCCGATGTTCTCTACACCTGCAGAGATCGCCCGTCTTGCCGTAGAAAACGACGTTCACGTAGTGGGAGCATCCTCACTGGCTGCCgttc                                                                                                                                                                                                                                  >  2:32615/1‑297 (MQ=255)                                                               tGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTATACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGa                                                                                                                                                                                                                      >  1:447670/1‑300 (MQ=255)                                                                        aaaGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAAc                                                                                                                                                                                                             <  1:516716/300‑1 (MQ=255)                                                                                     gTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCAGGAACTGGTCGAAGCGCt                                                                                                                                                                                                >  2:735235/1‑300 (MQ=255)                                                                                                              gATGCCATTGTGGCTCAAACGGAGCAGTTCCTTGCCGACAATGTTTGTCGACCGCTCATTCTGATTACTAAGACGGGCCAGGAGGGGCACGATCGCGGCGCGAAGTTGATCCCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACTTAGTGGGCGCATCCTCACTTGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATc                                                                                                                                                                       <  2:211935/300‑1 (MQ=255)                                                                                                                  ccATGGTTGCGAAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCCCGCACTCGTATCACTAAGATAGGTCAGGAGGGACACGATCGCGGCGCGAAGGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACTTAGATTTAAGCCCGATGTTCTCTACACCTTAAGAGATNGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACGGGTCGAAGCNCTGGAAAAATGGGGACGCGAAGATATTTGCg                                                                                                                                                                   <  2:586400/300‑1 (MQ=255)                                                                                                                                                                  cgcgCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAag                                                                                                                   <  1:703785/300‑1 (MQ=255)                                                                                                                                                                    cgcATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTNTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAagag                                                                                                                 <  1:566645/300‑1 (MQ=255)                                                                                                                                                                    cgcATCCTGATCGCTAAGATGGGCAAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGNCGTTCACGTAGTGGGCGCATCCTCACNGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAagag                                                                                                                 <  2:567534/300‑1 (MQ=255)                                                                                                                                                                                gCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggcg                                                                                                     >  1:297962/1‑300 (MQ=255)                                                                                                                                                                                 cTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggcgg                                                                                                    <  1:689951/300‑1 (MQ=255)                                                                                                                                                                                         ggCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAtttct                                                                                            >  1:52405/1‑298 (MQ=255)                                                                                                                                                                                         ggCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTAt                                                                                            >  1:90543/1‑300 (MQ=255)                                                                                                                                                                                         ggCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCCACGTCGCTCTCAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGCCGCGGGTGGCGTCATTCCGCAGCAGAATTACGCCTTCCTGCACGAGCGCGGCGTGGCGGCGAGTTAt                                                                                            >  2:215029/1‑300 (MQ=255)                                                                                                                                                                                           ccAGGAGGGACACGAGCGCGGGGCGAAAGTGGTCGACAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATgg                                                                                          <  2:50727/300‑1 (MQ=255)                                                                                                                                                                                              ggATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGCTTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCAGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTcc                                                                                       >  2:480864/1‑300 (MQ=255)                                                                                                                                                                                              ggATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTcc                                                                                       >  1:75173/1‑300 (MQ=255)                                                                                                                                                                                                       cgagatcggcgCGAAAGTGATCGCCGGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATATTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGTTCTAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGGCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATGTATGGTCCAGGTACAcc                                                                              <  2:181726/294‑1 (MQ=255)                                                                                                                                                                                                       cGATCGNGGCACGAAATTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTNGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACAcc                                                                              <  2:668347/300‑1 (MQ=255)                                                                                                                                                                                                       cGATCGCTGCGCGAAATTGATCGCCATCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTTAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGCCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACAcc                                                                              <  2:6305/300‑1 (MQ=255)                                                                                                                                                                                                       cGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACAcc                                                                              <  1:498086/300‑1 (MQ=255)                                                                                                                                                                                                       cGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCAGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACAcc                                                                              <  1:806473/300‑1 (MQ=255)                                                                                                                                                                                                         aTCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGGAGGTTTCAGCACGATGTTTTATACACGTGAAGAGGTTGACCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAAATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGAATTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTa                                                                            >  2:83848/1‑300 (MQ=255)                                                                                                                                                                                                               gcgcGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTAAACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCCGGTACACCTATGCTCg                                                                      >  1:283097/1‑300 (MQ=255)                                                                                                                                                                                                                 gcgAAAGTGGTCGCCAGCTCCTATTCCGAGCTCGGTTTCGACGTAGATGTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAGACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAANGCTGATCCCGGAACTGGTCGAAGCGCNGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCCTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGAc                                                                    <  2:567512/300‑1 (MQ=255)                                                                                                                                                                                                                 gcgAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGAc                                                                    <  1:41456/300‑1 (MQ=255)                                                                                                                                                                                                                      aGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTANTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAgtgt                                                               <  1:564907/300‑1 (MQ=255)                                                                                                                                                                                                                           tCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGcgcg                                                          <  1:409679/300‑1 (MQ=255)                                                                                                                                                                                                                                  cgcCTATTCCGATCTCGGTTTCGACGTCGATTTCAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCCGGTACACCTATGCTCGACAGTGTGCGCGACGTACt                                                   >  2:509874/1‑300 (MQ=255)                                                                                                                                                                                                                                  cgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGNTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACt                                                   >  1:553412/1‑300 (MQ=255)                                                                                                                                                                                                                                  cgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACt                                                   >  1:724417/1‑300 (MQ=255)                                                                                                                                                                                                                                      tATTCCGATCTCGGTTTCGACGTAGATTTAAGTCCGATGTTCTCTACACCTGAAGAAATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAGCGCTGATCCCGGAACTAGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAAt                                               <  2:155172/300‑1 (MQ=255)                                                                                                                                                                                                                                                gcGGTTTCGACGTAAATTTAAGCCCGATGGTCTCTACACCTGAATAGATCGCCCGCCTGGCCGTAGAAAACGACGGTCACGGAGTGGGCGCATCCTCACTGGCTGCCGGGCAGAAAACGTTGATCCCGGAACTGGTCGAAGCACTGAAAAAATGGGGACGCGAAGATATCTGCGGGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGcc                                     <  2:311532/299‑1 (MQ=255)                                                                                                                                                                                                                                                                              tctcTACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCa       >  2:99803/1‑300 (MQ=255)                                                                                                                                                                                                                                                                               ctctACACCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGGTTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTCTGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAg      >  2:643992/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                   acacCTGAAGAGATCGCCCGTCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAg  <  1:369214/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                        TTCCGATGCGCTGGAAGTCGCTTTCGACCGTTATCTGGTGCCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGTCGTCGCCCGCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAG  >  NC_000913/3062442‑3063015

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