Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	111,300	T→C	L86P (CTG→CCG)	mutT →	dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	111,300	0	T	C	100.0%	288.9 / NA	78	L86P (CTG→CCG)	mutT	dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/32);  total (46/32)

GGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGG  >  NC_000913/111152‑111444                                                                                                                                                     |                                                                                                                                                 ggcggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGGGCGGGAACTTCAGGGAAAAGTCGGGATTAACCCCCAAAATTTTTCGCTTTTTGAAAAACTGGAATATGAAATCCCGGACAGGCATATAACACTGTGGCTTTGGCCgg                                                                                                                                                >  1:39106/1‑151 (MQ=255)  gcggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGt                                                                                                                                               <  2:840532/151‑1 (MQ=255)  gcggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGt                                                                                                                                               <  1:477079/151‑1 (MQ=255)      tAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTAACCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGa                                                                                                                                            >  1:327683/1‑150 (MQ=255)         aaTTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTTAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGaa                                                                                                                                           >  2:569733/1‑148 (MQ=255)         aaTTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGaa                                                                                                                                           >  1:1236950/1‑148 (MQ=255)              aaaTGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAAATTCAGGAATAAGTCGGGATTAACCCCCAACATTTTTCGCTATTTGAAAAAATGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTggg                                                                                                                                    >  2:506664/1‑150 (MQ=255)                aTGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAg                                                                                                                                 <  2:1096407/151‑1 (MQ=255)                  gggTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAgg                                                                                                                                <  1:365689/150‑1 (MQ=255)                    gTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGa                                                                                                                             <  1:1141834/151‑1 (MQ=255)                    gTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGa                                                                                                                             <  2:1146244/151‑1 (MQ=255)                    gTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGGTTTGGCCGGTCGAACGCTGGGAAgggg                                                                                                                              >  1:1027277/1‑150 (MQ=255)                        aaCGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGACATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCg                                                                                                                         <  2:1193369/151‑1 (MQ=255)                        aaCGCCGGAACAGGCGGTGGGGCGTGAACTTTAGGAAGAAGTCGGGATTACCCCCCAACATTTTTTGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGGAGGGGAGCCg                                                                                                                         >  1:345049/1‑151 (MQ=255)                                  cAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAaag                                                                                                               <  1:403465/151‑1 (MQ=255)                                      cGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAgg                                                                                                           <  1:756223/151‑1 (MQ=255)                                       ggtggtGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAggg                                                                                                          >  2:865775/1‑151 (MQ=255)                                              cGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAGGAAGGGCAACCCg                                                                                                   >  1:450656/1‑151 (MQ=255)                                              cGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGAAAAGAAGGGCAACCCg                                                                                                   >  2:400024/1‑151 (MQ=255)                                                 gAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTg                                                                                                <  1:685123/151‑1 (MQ=255)                                                 gAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTg                                                                                                >  2:481838/1‑151 (MQ=255)                                                            agaagTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCg                                                                                     >  1:479210/1‑151 (MQ=255)                                                               agTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTg                                                                                  >  1:1052749/1‑151 (MQ=255)                                                                  cGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAACTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtc                                                                               <  2:5877/151‑1 (MQ=255)                                                                   gggATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtcg                                                                              >  1:142134/1‑151 (MQ=255)                                                                    ggATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtcgg                                                                             >  1:404298/1‑151 (MQ=255)                                                                    ggATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGATggtcgg                                                                             >  1:836826/1‑151 (MQ=255)                                                                     gATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtcagt                                                                            >  2:511173/1‑151 (MQ=255)                                                                      aTTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtcggtc                                                                           >  2:408562/1‑151 (MQ=255)                                                                      aTTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtcggtc                                                                           <  1:1090559/151‑1 (MQ=255)                                                                         aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTa                                                                        <  2:148240/151‑1 (MQ=255)                                                                          ccccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTa                                                                        >  2:452172/1‑150 (MQ=255)                                                                            ccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACGCTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATg                                                                     >  1:776699/1‑151 (MQ=255)                                                                            ccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATg                                                                     >  2:608204/1‑151 (MQ=255)                                                                            ccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGACGCTGGTCGGTCTTAATg                                                                     >  2:167268/1‑151 (MQ=255)                                                                                    tttttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCgatgat                                                             >  2:334969/1‑151 (MQ=255)                                                                                       ttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGAtttt                                                          >  1:202597/1‑151 (MQ=255)                                                                                         cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTcc                                                        <  1:813883/151‑1 (MQ=255)                                                                                         cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTcc                                                        >  2:514495/1‑151 (MQ=255)                                                                                             aTTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgcc                                                     >  1:222289/1‑150 (MQ=255)                                                                                             aTTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGACTTAATGCCGATGATTTTCCgcc                                                     >  1:587923/1‑150 (MQ=255)                                                                                                 gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcca                                                <  1:85775/151‑1 (MQ=255)                                                                                                 gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcca                                                <  2:654335/151‑1 (MQ=255)                                                                                                 gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcc                                                 <  1:14517/150‑1 (MQ=255)                                                                                                 gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcc                                                 <  2:164734/150‑1 (MQ=255)                                                                                                 gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcc                                                 >  1:626142/1‑150 (MQ=255)                                                                                                   aaaaCTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAAt                                              <  2:94502/151‑1 (MQ=255)                                                                                                     aaCTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATg                                             >  2:952843/1‑150 (MQ=255)                                                                                                     aaCGGGAATATGAATTCCCGGACATGCATATAACGCTGTGGTTTTGGCCGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGa                                            <  2:20586/151‑1 (MQ=255)                                                                                                           aaTATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGt                                      <  1:260218/151‑1 (MQ=255)                                                                                                            atatGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGt                                      >  1:1001147/1‑150 (MQ=255)                                                                                                            atatGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTa                                     <  2:43448/151‑1 (MQ=255)                                                                                                              atGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAt                                   >  1:1205204/1‑151 (MQ=255)                                                                                                                gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAtt                                  <  2:124810/150‑1 (MQ=255)                                                                                                                gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAtt                                  <  1:1234097/150‑1 (MQ=255)                                                                                                                gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAtt                                  >  2:208693/1‑150 (MQ=255)                                                                                                                gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTg                                 >  2:1078588/1‑151 (MQ=255)                                                                                                                 aaTTCCCGGACAGGCATATAACCCTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTg                                 >  1:275114/1‑150 (MQ=255)                                                                                                                  aTTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGc                                >  2:467272/1‑150 (MQ=255)                                                                                                                  aTTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCg                               >  2:332785/1‑151 (MQ=255)                                                                                                                      ccGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGc                           >  1:783146/1‑151 (MQ=255)                                                                                                                        ggACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCtt                         >  2:957526/1‑151 (MQ=255)                                                                                                                         gACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCtt                         >  2:582502/1‑150 (MQ=255)                                                                                                                          aCAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTaa                       <  1:340466/151‑1 (MQ=255)                                                                                                                             ggCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAAc                     <  1:707552/150‑1 (MQ=255)                                                                                                                                atatAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCt                 <  2:126026/151‑1 (MQ=255)                                                                                                                                atatAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCt                 >  2:562379/1‑151 (MQ=255)                                                                                                                                atatAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCt                 >  1:50449/1‑151 (MQ=255)                                                                                                                                   taACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGt               >  2:932898/1‑150 (MQ=255)                                                                                                                                    aaCACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGAGGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAg             <  1:985271/151‑1 (MQ=255)                                                                                                                                     acacTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAgg            >  2:905612/1‑151 (MQ=255)                                                                                                                                      cacTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGt           <  1:211793/151‑1 (MQ=255)                                                                                                                                         tgtgGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTc          >  1:228885/1‑149 (MQ=255)                                                                                                                                          gtgGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCa         <  2:1083778/149‑1 (MQ=255)                                                                                                                                          gtgGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAAGGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGa       <  1:52993/151‑1 (MQ=255)                                                                                                                                             gTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATaa    <  1:700394/151‑1 (MQ=255)                                                                                                                                             gTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATa     >  2:583281/1‑150 (MQ=255)                                                                                                                                               tttGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAgg  >  2:63540/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGG  >  NC_000913/111152‑111444

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