Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,916	C→T	A356V (GCC→GTC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,916	0	C	T	95.3%	189.1 / ‑5.8	64	A356V (GCC→GTC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (34/27);  total (36/28)
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

CGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGG  >  NC_000913/3061670‑3062161                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGACCGACGACCTGTTGCGCGCGGCACGTTATTTCTGAAGCGAAGCGGTCAGTGGATTTGGCCCACAGGACCCGCAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCCCCACCATTGAAGCGCTGGCTGCGCCGCTGGGCGGTACTCAGTCACTGCCTACCAACGCCttt                                                                                                                                                                                                                                                   >  2:38550/1‑251 (MQ=255)  gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCTCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGGTTTGGCGCACAGGACCCGGAATCACTGGCGCCGCGTACCCACTGCCAGACCTCCGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTAATCAGTCACTGCATACCAACGTCtttt                                                                                                                                                                                                                                                  >  1:299728/1‑250 (MQ=255)              catcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCCGTCACTGCATACCAACGTCTTTGACGAAGCGCTTg                                                                                                                                                                                                                                      >  1:208573/1‑251 (MQ=255)               atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTgg                                                                                                                                                                                                                                     <  2:84876/251‑1 (MQ=255)                  ggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGTCGCTGGGCGGTACTCATTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGttt                                                                                                                                                                                                                                  <  2:253086/251‑1 (MQ=255)                         aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTAcc                                                                                                                                                                                                                           >  1:126299/1‑251 (MQ=255)                          tCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATCCCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                                                                                                          >  1:279284/1‑251 (MQ=255)                          tCTGTCTATGAACATTGACATGTTGCGAGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTAGCGCACAGGACCCGAAATCACTGGCGCTGCGTAACCCCTGCAAGACCTAAGGCTGGAGCCTCACTGAACAGGATCCGTATAACAACGTTATCCGCAACACCATTGAAGCGCTGGCTGCGACGCCGGGCGGTACTCAGTCACTGCCTACCAAAGTCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                                                                                                          >  2:23399/1‑251 (MQ=255)                               ttATGTACGTCGCTATGTTTCGTGCGGCACGATCTTTATGGAGCGAAACGGTCCGTGGATTTGCCGCACAGGCCCCGAAATCACTGGCGCTGCGTACCCACTGCACGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGCCGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtc                                                                                                                                                                                                                     >  2:391950/1‑251 (MQ=255)                               ttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCTTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACGCAGTCACTGCATACCAACGTCTTTGACGAATCGCTTGGTTTGCCTACCGAttt                                                                                                                                                                                                                      <  2:78794/250‑1 (MQ=255)                                  ggAACGTCGCCATGTTGCGTGCGGCACGTTAGTTATGGAGCGAAGCGGGCAGGGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATGGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                                                                                                                                                                                                  <  2:86211/250‑1 (MQ=255)                                    aaCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGc                                                                                                                                                                                                                <  1:224688/251‑1 (MQ=255)                                      cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGAGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCAACATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAAAGTCTTTGACGAAGCGCTTGGGTTGCCTACCGATTTCTCCGCAc                                                                                                                                                                                                              >  1:393345/1‑251 (MQ=255)                                         cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGCGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGc                                                                                                                                                                                                            >  2:120154/1‑250 (MQ=255)                                         cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAg                                                                                                                                                                                                                                            >  2:317706/1‑218 (MQ=255)                                         cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAg                                                                                                                                                                                                                                            <  1:317706/218‑1 (MQ=255)                                         cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTAGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGc                                                                                                                                                                                                            >  1:142824/1‑250 (MQ=255)                                          gCCAGGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATGTGGCGCACATGACCCGAAAGCACGGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCGCCACCATTGAAGCGCTGGCTGCGACGCTGTGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACTGATTTCTCAGCACGCAt                                                                                                                                                                                                          <  1:356051/251‑1 (MQ=255)                                                ttGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTTGATTTGGCCCACAGGACCCGAAATCACTGGCGCTGCGTACCCCCTGCCAGACCTCATGCTGGAGCCTGACTTAACAGGATCCGTATAACAACGTTATCCGCCCCACCATTGGAGCGCTGGCTGCGTCGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGccc                                                                                                                                                                                                     <  1:129753/250‑1 (MQ=255)                                                     tgcgGAACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACGGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGAAGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAac                                                                                                                                                                                                >  2:244883/1‑250 (MQ=255)                                                         gCACGGTATTTCTGGAGCGTACAGGTCAGTGAATTAGGCGCACAGGACACGACCTCACTGGAGCTGCGCACCCACTGCCAGACCTCAGGCTGGAGCCTGCCTGAAAAGGTTCCGTATAACAACGCTATCCGCCCCACCATTGAAGAGCTGGCTGCGACGATGGGCAGTACTCAGTCACTGCATACCAACGCCTTTGACGACGCGGTCGGTTTGCCTACCGATTTCTCAGCGCGCCTTGCCCGACACAccc                                                                                                                                                                                            >  2:232000/1‑250 (MQ=255)                                                            cGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACAccc                                                                                                                                                                                            >  1:304880/1‑247 (MQ=255)                                                             gTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAGTGAAGCGCTGGCTGCGACGCTGGTCGGTACTCAGTCACTGCATACCAACGTCGTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCGCGCACCACCCAGatc                                                                                                                                                                                       <  2:167754/251‑1 (MQ=255)                                                               taattatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACGGGCGCTGCTTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACATGATCCGTAGAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCACACCAACGTCTTTGACTAAGCGCTTGGTTTGCCTACCGATTCCTCAGCACGCATTGCCCGCAACACCCAGatca                                                                                                                                                                                      <  2:403256/247‑1 (MQ=255)                                                                    gtGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATcc                                                                                                                                                                                <  1:146889/250‑1 (MQ=255)                                                                                        gATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAGCCGTATAACAACGTTCTCCGCACCACCATTGAAGCGCTGGATGCGCCGCTGGGCGGTCCTCAGTCACTGCATACCTACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGACGAATCAGAACTCTGc                                                                                                                                                            >  1:397102/1‑251 (MQ=255)                                                                                           ttGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAAGTCTGCCg                                                                                                                                                          <  1:106686/250‑1 (MQ=255)                                                                                                      gACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATcc                                                                                                                                              <  2:279661/251‑1 (MQ=255)                                                                                                         ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACt                                                                                                                                           <  1:361851/251‑1 (MQ=255)                                                                                                              aTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACAGAACAGGATTCGTATAATAACGTTATCCGCACCGCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTGTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCAACATCCAGGAAGAATCAGAACTCTGCCGCTCCGTCGATCCACTGGCCg                                                                                                                                      <  1:232000/251‑1 (MQ=255)                                                                                                                     gcgcTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAAGCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTa                                                                                                                               <  2:140525/251‑1 (MQ=255)                                                                                                                       gcTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTAtt                                                                                                                             <  1:165126/251‑1 (MQ=255)                                                                                                                               cccACTGCCAGACCGAAGGCTGGAGCCTGACTGAACAGGAGCCGTATAACAACGTTATCCGCACCACCATGGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTTACGAAGCGATGGGTTGGCCTACCGATTTCGCAGCACGCATGGCCCGCAACACCCAGATCATCATCCAGGAAGAATCATAACTCTGCCGCACCGTCGATCCACTGGCCGGCTCCTATTACAGTGAg                                                                                                                     <  2:258667/251‑1 (MQ=255)                                                                                                                                     gCCAGACCTCAGGATGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTg                                                                                                               >  2:198446/1‑251 (MQ=255)                                                                                                                                          aCCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCg                                                                                                                                                                                                    >  2:360960/1‑161 (MQ=255)                                                                                                                                          aCCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACGCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCg                                                                                                                                                                                                    <  1:360960/161‑1 (MQ=255)                                                                                                                                              cAGGCTGGAGCTTGACTGAGCCGGATCCTTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATATCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCAGCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCaa                                                                                                      <  2:154895/251‑1 (MQ=255)                                                                                                                                               cggCTGGAGCCGGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTGCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGACGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCaaa                                                                                                     <  2:21589/250‑1 (MQ=255)                                                                                                                                                   tGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGt                                                                                                 >  1:150421/1‑251 (MQ=255)                                                                                                                                                   tGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCAGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGCCCACTCAAATCgc                                                                                                 >  2:324060/1‑250 (MQ=255)                                                                                                                                                         ctgactgaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGAGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGCCCGATCAAATCGTCAAAc                                                                                            >  2:360729/1‑250 (MQ=255)                                                                                                                                                           gactgaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGATGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGCTCCTATTACCTTGAGTCGCTGACCGATCACATCGTCAAACa                                                                                           >  1:210352/1‑249 (MQ=255)                                                                                                                                                             ctgaACAAGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGcc                                                                                       >  1:213807/1‑251 (MQ=255)                                                                                                                                                               gaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACAAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATAGTCAAACAAGcc                                                                                       >  1:395040/1‑249 (MQ=255)                                                                                                                                                                  cAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGc                                                                                  >  1:295470/1‑251 (MQ=255)                                                                                                                                                                   aGGATCCATATAACAACGTTATCCGCCCCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATAATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATAGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCt                                                                                 >  2:268765/1‑251 (MQ=255)                                                                                                                                                                    ggATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCta                                                                                <  1:120154/251‑1 (MQ=255)                                                                                                                                                                      aTCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCtatt                                                                              >  1:292539/1‑251 (MQ=255)                                                                                                                                                                              aacaacGTTATCCGCCCCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATGTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGAGCCTATTACATTGAGTCGCGGACCGATCAAATCGTCAAACAAGCCAGAGCTATGATCCAAc                                                                       <  2:208573/250‑1 (MQ=255)                                                                                                                                                                                  acGTTATCCGCATCACCATCGAAGCGCTGGCTGCGACGCTGGGCGGTGCTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCTACACCCAGATCATCATCCAGGAAGAATCAGATCTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGATTCGCAGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATc                                                                  <  1:244156/251‑1 (MQ=255)                                                                                                                                                                                      tATCCGCACCTCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACCCCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTAATACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCTAACAGATcgac                                                               <  2:287996/250‑1 (MQ=255)                                                                                                                                                                                            caccacAATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCCTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTCTTATCCAACAGATCGACGAAGCCg                                                        >  2:323870/1‑251 (MQ=255)                                                                                                                                                                                              ccaccaTTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGt                                                      >  1:275813/1‑251 (MQ=255)                                                                                                                                                                                               caccaTTGAAGCGCTGGCTGCGACGCTGGGCGGTCCTCAGTCACTGCATACCAACGTCTTTGACGGAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGCAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGCCCGATCAAATCGTCAAAAAAGCCAGAGCTATTCTCCAACAGATCGCCGAAGCCgatg                                                     >  1:160671/1‑248 (MQ=255)                                                                                                                                                                                                   gttGAAGGGCTGGCTGCGACGCTGGACGGTACTCAATCACTTCATCCCACCGTATTTGACGAAGCGCTTGGTTTGCCTACAGATTTCTCAGCACGCATTGCCCGCCACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGCTCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACCAGCCAGCGCTATTATCCACCCTATCAACGAAGCCggggaaat                                                 >  2:232879/2‑245 (MQ=255)                                                                                                                                                                                                         gcgctagcaacgcccccgggCGGTCTTCAGACACTGCATACTAACGTCCTTGACAACGGGCTTGGTCTCCCTACGGATTTCTCAGCACGAATTGCCCGGAACCCCCAGCTCATCATACAGGCAGCATCAGAACTATGACGCAACGTCGCTCCACTGGCCGGATCCTCTTAAATTGAGTCGCTGAACGATCAAATCGTCAAACAAGCCAGAGCTATTATACAAAAGATCGACGCAGCAGGAGGCATGGCGa                                            >  2:374240/18‑250 (MQ=255)                                                                                                                                                                                                                ctgCGTCGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGCTCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGCTCGACGAAGCCGGTGGCATGGCGAAAGCGATc                                    >  2:378361/1‑251 (MQ=255)                                                                                                                                                                                                                  gCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCg                                   <  1:350192/250‑1 (MQ=255)                                                                                                                                                                                                                    gACGCTGAGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAg                                <  1:360729/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATg            >  1:134892/1‑250 (MQ=255)                                                                                                                                                                                                                                         aCTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCACAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGCCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGCTCAAATAGTCAACAAAGCCAGCGCTATTATCCAACAGATAGACGAAACCAGTGGCCAGGCGAACGAGATCGAAGAAGGTCTGCCAAAACGGATg            >  2:388686/1‑250 (MQ=255)                                                                                                                                                                                                                                          cTGCATACCAACGTCTTTGCCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCACAACGAATGAt          >  1:150550/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGGCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGAt          <  2:304880/249‑1 (MQ=255)                                                                                                                                                                                                                                                  cAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAgg  >  1:274664/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGG  >  NC_000913/3061670‑3062161

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