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%	297.5 / NA	79	T650T (ACG→ACA)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/47);  total (32/47)

CGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTC  >  NC_000913/3062651‑3062946                                                                                                                                                     |                                                                                                                                                    cgAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACt                                                                                                                                                   >  1:395997/1‑151 (MQ=255)      gTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATcc                                                                                                                                              >  2:529815/1‑151 (MQ=255)           cGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGaa                                                                                                                                         <  1:321467/151‑1 (MQ=255)            gCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGaa                                                                                                                                         >  1:664222/1‑150 (MQ=255)             cgagagcCTATTCCGATCGCGGTTTCGACGTAGATTTAAGCACGATGTTCTCTCCACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAAc                                                                                                                                        <  2:380496/145‑1 (MQ=255)             ccAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACt                                                                                                                                       <  1:475032/151‑1 (MQ=255)              cAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTg                                                                                                                                      <  1:61870/151‑1 (MQ=255)              cAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTg                                                                                                                                      >  1:405166/1‑151 (MQ=255)                           gATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTg                                                                                                                          >  2:206822/1‑150 (MQ=255)                              ctcGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaa                                                                                                                       >  2:676068/1‑150 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGaa                                                                                                                                <  2:124489/139‑1 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaaa                                                                                                                      <  1:137168/149‑1 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaa                                                                                                                       <  1:366100/148‑1 (MQ=255)                                cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGa                                                                                                                         <  1:637721/146‑1 (MQ=255)                                 ggTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaaaaa                                                                                                                    >  1:455347/1‑150 (MQ=255)                                   tttCGACGTAGATTTAAGCCCGCTGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATgg                                                                                                                 >  1:312449/1‑151 (MQ=255)                                      cGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGa                                                                                                              <  1:620504/151‑1 (MQ=255)                                       gACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAc                                                                                                             >  2:209133/1‑151 (MQ=255)                                        aCGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcg                                                                                                            >  2:732765/1‑151 (MQ=255)                                          gTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcgcg                                                                                                          <  1:262114/151‑1 (MQ=255)                                           tAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcgcg                                                                                                          >  1:357783/1‑150 (MQ=255)                                            aGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGa                                                                                                         >  2:129232/1‑150 (MQ=255)                                                ttAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGat                                                                                                     <  2:395976/150‑1 (MQ=255)                                                ttAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGat                                                                                                     >  2:43058/1‑150 (MQ=255)                                                 tAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGatat                                                                                                   >  1:446408/1‑151 (MQ=255)                                                  aaGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGatat                                                                                                   <  1:206830/150‑1 (MQ=255)                                                   aGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCt                                                                                                 <  2:627587/151‑1 (MQ=255)                                                     cccGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGc                                                                                               >  2:757219/1‑151 (MQ=255)                                                       cGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGt                                                                                             <  1:1493/151‑1 (MQ=255)                                                        gATGTTCTCTTCACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCACGGAACTGGTGGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTg                                                                                                <  1:328182/147‑1 (MQ=255)                                                           gTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTc                                                                                         >  1:182845/1‑151 (MQ=255)                                                                 tACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCggg                                                                                    >  2:447797/1‑150 (MQ=255)                                                                  acacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGt                                                                                   <  1:526612/150‑1 (MQ=255)                                                                  acacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGt                                                                                   <  1:458311/150‑1 (MQ=255)                                                                   cacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCggg                                                                                    <  2:764549/148‑1 (MQ=255)                                                                    acCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGc                                                                                >  2:287140/1‑151 (MQ=255)                                                                     ccTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTgg                                                                                 <  2:625217/149‑1 (MQ=255)                                                                        gAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTc                                                                             <  1:93066/150‑1 (MQ=255)                                                                           gagaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTc                                                                         >  1:133388/1‑151 (MQ=255)                                                                            agaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTcc                                                                        >  1:179175/1‑151 (MQ=255)                                                                            agaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTcc                                                                        <  1:604312/151‑1 (MQ=255)                                                                                cgcccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCGCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgccg                                                                    >  1:172386/1‑151 (MQ=255)                                                                                 gcccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCcgccgc                                                                   >  2:618282/1‑151 (MQ=255)                                                                                 gcccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCcgccgc                                                                   <  1:402689/151‑1 (MQ=255)                                                                                    cgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAgg                                                                <  2:522381/151‑1 (MQ=255)                                                                                      ccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGAt                                                              <  2:763160/151‑1 (MQ=255)                                                                                        tGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGAtt                                                             >  2:647303/1‑150 (MQ=255)                                                                                         ggCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGAtt                                                             <  2:667600/149‑1 (MQ=255)                                                                                               aGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCtt                                                      <  2:649346/150‑1 (MQ=255)                                                                                                gAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTcc                                                    <  1:351707/151‑1 (MQ=255)                                                                                                gAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTcc                                                    >  2:583577/1‑151 (MQ=255)                                                                                                 aaaaCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCt                                                   <  2:201018/151‑1 (MQ=255)                                                                                                      gacgTTCACGTAGTGGGCGCATCCTCCCTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAag                                              <  1:454236/151‑1 (MQ=255)                                                                                                       acgTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAaga                                             <  2:730281/151‑1 (MQ=255)                                                                                                        cgTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAaga                                             <  2:102392/150‑1 (MQ=255)                                                                                                         gTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAgc                                           <  2:305996/151‑1 (MQ=255)                                                                                                           tCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAgcgc                                         >  1:466554/1‑151 (MQ=255)                                                                                                            cACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGcgcg                                        <  1:593149/151‑1 (MQ=255)                                                                                                               gTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCgcggcg                                     >  1:495020/1‑151 (MQ=255)                                                                                                                 aGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGt                                    <  1:693174/150‑1 (MQ=255)                                                                                                                    gggCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggcg                                <  2:406363/151‑1 (MQ=255)                                                                                                                    gggCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggca                                <  2:425144/151‑2 (MQ=255)                                                                                                                    gggCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggc                                 <  2:474313/150‑1 (MQ=255)                                                                                                                      gcgcATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTggcggc                              >  1:644434/1‑151 (MQ=255)                                                                                                                          aTCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAtt                          >  1:466568/1‑151 (MQ=255)                                                                                                                            ccTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAttt                         <  2:81435/150‑1 (MQ=255)                                                                                                                            ccTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTa                        <  1:552232/151‑1 (MQ=255)                                                                                                                               cACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATg                      <  2:385838/150‑1 (MQ=255)                                                                                                                                   ggCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTcc                  >  1:263238/1‑150 (MQ=255)                                                                                                                                      tGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGt              >  1:354775/1‑151 (MQ=255)                                                                                                                                            tCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCt        >  1:341466/1‑151 (MQ=255)                                                                                                                                             cATAAACCACTGATCCCGGAGCTGGTCGCAGCGTTGAGAAAATGGGGACGAGAAGATATCTGCGTGGTCGTGATTGGCGTCATTCCGACGCAGGATTTCGCTTTCCTGTAAGAGCGCGGCGTGGCGGCGATTTATGGTCCTGGTACCCCt        <  2:571239/150‑1 (MQ=255)                                                                                                                                             cATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCt        <  1:310642/150‑1 (MQ=255)                                                                                                                                             cATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTa       <  2:344870/151‑1 (MQ=255)                                                                                                                                              aTAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTAt      <  1:529841/151‑1 (MQ=255)                                                                                                                                                aaaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATg     <  1:486654/150‑1 (MQ=255)                                                                                                                                                aaaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGc    <  1:28212/151‑1 (MQ=255)                                                                                                                                                 aaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCt   <  2:590426/151‑1 (MQ=255)                                                                                                                                                  aaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTc  >  2:196593/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTC  >  NC_000913/3062651‑3062946

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