Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,799	G→A	T650T (ACG→ACA)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,799	0	G	A	100.0%	315.9 / NA	81	T650T (ACG→ACA)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/42);  total (39/42)

CGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTC  >  NC_000913/3062651‑3062946                                                                                                                                                     |                                                                                                                                                    cgAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACt                                                                                                                                                   >  2:845713/1‑151 (MQ=255)   aaaGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGa                                                                                                                                                 <  2:91310/151‑1 (MQ=255)     aGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATc                                                                                                                                               >  2:127333/1‑151 (MQ=255)      gTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATcc                                                                                                                                              >  1:595037/1‑151 (MQ=255)           cGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGaa                                                                                                                                         <  2:659093/151‑1 (MQ=255)           cGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGa                                                                                                                                          >  2:697885/1‑150 (MQ=255)               aGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTgg                                                                                                                                     <  1:849574/151‑1 (MQ=255)                gcgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGt                                                                                                                                    <  1:17319/151‑1 (MQ=255)                 cgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTc                                                                                                                                   >  1:432892/1‑151 (MQ=255)                  gcCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCg                                                                                                                                  <  2:833223/151‑1 (MQ=255)                     tATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGa                                                                                                                                 >  2:173367/1‑149 (MQ=255)                               tcGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaa                                                                                                                       >  1:276287/1‑149 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAgcgc                                                                                                                            <  1:36899/143‑1 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTg                                                                                                                          <  1:890309/145‑1 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaaaa                                                                                                                     >  1:592170/1‑150 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaaa                                                                                                                      <  2:818718/149‑1 (MQ=255)                                 ggTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAAt                                                                                                                   >  1:29409/1‑151 (MQ=255)                                 ggTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAAt                                                                                                                   >  2:532411/1‑151 (MQ=255)                                     ttgACGTAGATTTAACCCTGATGTTCTATGCCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGAAACTGGTCGAAGCGCTGAAAAAATgg                                                                                                                 <  1:684538/147‑1 (MQ=255)                                     tCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATgggg                                                                                                               >  2:1020655/1‑151 (MQ=255)                                      cGACTTAGATTTATGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGGTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAgc                                                                                                                              <  2:450942/135‑1 (MQ=255)                                        aCGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcg                                                                                                            >  1:437078/1‑151 (MQ=255)                                          gTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcgcg                                                                                                          <  1:820209/151‑1 (MQ=255)                                          gTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcgcg                                                                                                          <  2:725252/151‑1 (MQ=255)                                              aTTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGa                                                                                                         >  1:327055/1‑148 (MQ=255)                                                 tAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGatat                                                                                                   >  1:222819/1‑151 (MQ=255)                                                   aGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCt                                                                                                 >  1:639528/1‑151 (MQ=255)                                                    gCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTg                                                                                                >  1:733457/1‑151 (MQ=255)                                                     cccGATGTTCTCTACACCAGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGAATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTg                                                                                                <  2:822633/150‑1 (MQ=255)                                                      ccGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCg                                                                                              >  2:870071/1‑151 (MQ=255)                                                       cGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGACCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGt                                                                                             <  1:283816/151‑1 (MQ=255)                                                        gATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTg                                                                                            >  2:315018/1‑151 (MQ=255)                                                         aTGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTgg                                                                                           >  1:782104/1‑151 (MQ=255)                                                           gTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTc                                                                                         >  2:813085/1‑151 (MQ=255)                                                                ctACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTcgcg                                                                                      <  2:559500/149‑1 (MQ=255)                                                                  acacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGt                                                                                   <  2:53508/150‑1 (MQ=255)                                                                       tGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTc                                                                             >  1:507466/1‑151 (MQ=255)                                                                       tGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTc                                                                             >  2:881467/1‑151 (MQ=255)                                                                        gAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTc                                                                             >  1:443170/1‑150 (MQ=255)                                                                           gagaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTc                                                                         <  2:227794/151‑1 (MQ=255)                                                                              aTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgc                                                                      >  2:531183/1‑151 (MQ=255)                                                                              aTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgc                                                                      <  2:738934/151‑1 (MQ=255)                                                                               tCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgcc                                                                     >  1:254548/1‑151 (MQ=255)                                                                                cgcccgccTGGCCGTAGACAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgccg                                                                    >  1:635974/1‑151 (MQ=255)                                                                                cgcccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgccg                                                                    >  2:112668/1‑151 (MQ=255)                                                                                 gcccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCcgccgc                                                                   <  1:198809/151‑1 (MQ=255)                                                                                  cccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCa                                                                  <  2:76474/151‑1 (MQ=255)                                                                                     gccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAgg                                                                >  1:456013/1‑150 (MQ=255)                                                                                     gccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGa                                                               <  2:839890/151‑1 (MQ=255)                                                                                        tGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGAt                                                              >  1:925064/1‑149 (MQ=255)                                                                                         ggCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTa                                                            <  1:134060/150‑1 (MQ=255)                                                                                               aGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCtt                                                      <  1:866029/150‑1 (MQ=255)                                                                                                 aaaaCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCt                                                   <  1:544789/151‑1 (MQ=255)                                                                                                   aaCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGc                                                 >  2:818812/1‑151 (MQ=255)                                                                                                   aaCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGc                                                 <  1:676299/151‑1 (MQ=255)                                                                                                       acgTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCa                                                >  1:406637/1‑148 (MQ=255)                                                                                                         gTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAgc                                           <  2:254433/151‑1 (MQ=255)                                                                                                         gTTCACGTAGTGGGCGCATCCTCACCGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAgc                                           <  1:427332/151‑1 (MQ=255)                                                                                                               gTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCgcggcg                                     >  2:1039036/1‑151 (MQ=255)                                                                                                               gTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCgcggcg                                     >  1:664335/1‑151 (MQ=255)                                                                                                                 aGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGt                                    <  2:338484/150‑1 (MQ=255)                                                                                                                 aGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTg                                   >  1:533521/1‑151 (MQ=255)                                                                                                                    gggCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggc                                 <  1:364306/150‑1 (MQ=255)                                                                                                                     ggCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggcgg                               <  1:870104/151‑1 (MQ=255)                                                                                                                      gcgcATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGGAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggcggc                              <  1:782689/151‑1 (MQ=255)                                                                                                                          aTCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAtt                          <  2:606718/151‑1 (MQ=255)                                                                                                                          aTCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAtt                          <  2:1996/151‑1 (MQ=255)                                                                                                                               cACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTAt                       <  2:651297/149‑1 (MQ=255)                                                                                                                                  tGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTcc                  >  2:979904/1‑151 (MQ=255)                                                                                                                                    gCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAg                <  1:201625/151‑1 (MQ=255)                                                                                                                                      tGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAATCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAgg               >  1:182493/1‑150 (MQ=255)                                                                                                                                       gCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGt              >  2:214216/1‑150 (MQ=255)                                                                                                                                        ccGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTac            <  1:389413/151‑1 (MQ=255)                                                                                                                                        ccGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTa             <  2:426264/150‑1 (MQ=255)                                                                                                                                         cGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTaca           <  2:763528/151‑1 (MQ=255)                                                                                                                                         cGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTaca           <  2:704700/151‑1 (MQ=255)                                                                                                                                           gTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACAcc         <  2:456002/151‑1 (MQ=255)                                                                                                                                              aTAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTAt      <  1:233856/151‑1 (MQ=255)                                                                                                                                                 aaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGc    >  1:874861/1‑150 (MQ=255)                                                                                                                                                 aaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCt   >  1:191250/1‑151 (MQ=255)                                                                                                                                                  aaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTc  <  2:449734/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTC  >  NC_000913/3062651‑3062946

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