Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,826,054	T→C	E290E (GAA→GAG)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,054	0	T	C	100.0%	294.0 / NA	79	E290E (GAA→GAG)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/38);  total (41/38)

TTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCG  >  NC_000913/1825906‑1826201                                                                                                                                                     |                                                                                                                                                    ttATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTc                                                                                                                                                   >  2:150721/1‑151 (MQ=255)  tataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTcc                                                                                                                                                  >  1:168292/1‑151 (MQ=255)   ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCct                                                                                                                                                 <  1:1012510/151‑1 (MQ=255)   ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCct                                                                                                                                                 >  1:635706/1‑151 (MQ=255)       gAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTcc                                                                                                                                             >  1:785708/1‑151 (MQ=255)             aaCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTg                                                                                                                                       >  1:756463/1‑151 (MQ=255)               cTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTgcg                                                                                                                                     >  2:766989/1‑151 (MQ=255)                 tAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTgcgc                                                                                                                                    <  2:827267/150‑1 (MQ=255)                  aaTACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGcc                                                                                                                                   >  1:52032/1‑150 (MQ=255)                   aTACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAg                                                                                                                                 <  2:103829/151‑1 (MQ=255)                     aCCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGc                                                                                                                                >  1:94717/1‑150 (MQ=255)                       ccGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAAt                                                                                                                             <  1:223872/151‑1 (MQ=255)                       ccGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAAt                                                                                                                             <  2:677136/151‑1 (MQ=255)                        cGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTGTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATg                                                                                                                            <  1:1110880/151‑1 (MQ=255)                         gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGt                                                                                                                           >  2:61192/1‑151 (MQ=255)                             aaTGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCt                                                                                                                       <  2:1028481/151‑1 (MQ=255)                                gATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCttt                                                                                                                     <  2:94718/150‑1 (MQ=255)                                gATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtc                                                                                                                    >  2:1195199/1‑151 (MQ=255)                                 aTTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct                                                                                                                   <  2:955397/151‑1 (MQ=255)                                 aTTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct                                                                                                                   >  2:597361/1‑151 (MQ=255)                                   ttCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtctc                                                                                                                  >  2:270084/1‑150 (MQ=255)                                       gCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCaa                                                                                                                <  1:397113/148‑1 (MQ=255)                                       gCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACg                                                                                                             >  1:1183692/1‑151 (MQ=255)                                        cGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACgg                                                                                                            <  2:141362/151‑1 (MQ=255)                                        cGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACgg                                                                                                            <  1:506583/151‑1 (MQ=255)                                            tAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCAt                                                                                                         >  1:1055110/1‑150 (MQ=255)                                              aGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATaaa                                                                                                      <  1:912870/151‑1 (MQ=255)                                                 aTTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAAtt                                                                                                   >  1:129761/1‑151 (MQ=255)                                                    aGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTc                                                                                                  >  1:235290/1‑149 (MQ=255)                                                    aGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAg                                                                                                <  1:483661/151‑1 (MQ=255)                                                       ttATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCgt                                                                                             >  1:372514/1‑151 (MQ=255)                                                                    cATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCata                                                                                <  2:331507/151‑1 (MQ=255)                                                                       tAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGc                                                                             >  1:1119003/1‑151 (MQ=255)                                                                        aaTCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGc                                                                             >  1:372462/1‑150 (MQ=255)                                                                        aaTCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGc                                                                             <  2:384767/150‑1 (MQ=255)                                                                            ccGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTc                                                                        >  2:22358/1‑151 (MQ=255)                                                                            ccGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTc                                                                        <  2:1007618/151‑1 (MQ=255)                                                                             cgcgcgTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTAACTTGCCATATGCATTTCg                                                                       <  2:577406/151‑1 (MQ=255)                                                                               cgcgTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGaa                                                                     <  1:1195168/151‑1 (MQ=255)                                                                                gcgTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGa                                                                      >  1:942149/1‑149 (MQ=255)                                                                                gcgTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGa                                                                      >  2:577306/1‑149 (MQ=255)                                                                                  gTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAgg                                                                   >  2:58159/1‑150 (MQ=255)                                                                                    aaTCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGAc                                                                 <  2:304072/150‑1 (MQ=255)                                                                                    aaTCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACg                                                                <  2:1055135/151‑1 (MQ=255)                                                                                       cccAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggc                                                              <  2:1140205/150‑1 (MQ=255)                                                                                          aaCGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagt                                                          >  1:1003419/1‑151 (MQ=255)                                                                                           tcGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATACAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtg                                                         >  1:226699/2‑151 (MQ=255)                                                                                              cTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcg                                                      >  1:811023/1‑151 (MQ=255)                                                                                                aCCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcgag                                                    <  2:736502/151‑1 (MQ=255)                                                                                                aCCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcg                                                      >  1:638199/1‑149 (MQ=255)                                                                                                  cAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcgagt                                                   >  2:444992/1‑150 (MQ=255)                                                                                                  cAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTa                                                  <  2:1158765/151‑1 (MQ=255)                                                                                                   aaCGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcg                                                             <  1:1190360/139‑1 (MQ=255)                                                                                                   aaCGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcg                                                             >  2:1190391/1‑139 (MQ=255)                                                                                                   aaCGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTaa                                                 <  1:97727/151‑1 (MQ=255)                                                                                                    aCGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAAt                                                <  1:114249/151‑1 (MQ=255)                                                                                                     cGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAAtt                                               <  1:215416/151‑1 (MQ=255)                                                                                                         ttAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCg                                           >  2:513518/1‑151 (MQ=255)                                                                                                             ggAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac                                       <  2:693821/151‑1 (MQ=255)                                                                                                              gAATAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAacc                                      >  1:97065/1‑151 (MQ=255)                                                                                                               aaTAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAacca                                     <  1:756815/151‑1 (MQ=255)                                                                                                                aTAACACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAaccac                                    >  2:949496/1‑151 (MQ=255)                                                                                                                  aaCACATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCAccc                                  <  1:270081/151‑1 (MQ=255)                                                                                                                   acacATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCAccc                                  >  1:174766/1‑150 (MQ=255)                                                                                                                   acacATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCg                                 >  2:1010090/1‑151 (MQ=255)                                                                                                                   acacATACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCg                                 <  2:992889/151‑1 (MQ=255)                                                                                                                    cacaTACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCg                                 >  1:1070163/1‑150 (MQ=255)                                                                                                                     acaTACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt                               >  1:565396/1‑151 (MQ=255)                                                                                                                     acaTACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt                               >  2:842850/1‑151 (MQ=255)                                                                                                                      caTACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATa                              <  1:1126802/151‑1 (MQ=255)                                                                                                                       aTACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATaa                             <  1:929021/151‑1 (MQ=255)                                                                                                                       aTACTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATaa                             <  2:12331/151‑1 (MQ=255)                                                                                                                         aCTCTACATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACg                           >  1:310542/1‑151 (MQ=255)                                                                                                                              aCATCATGGGTTACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAg                        <  1:330309/149‑1 (MQ=255)                                                                                                                                         tACGGTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTtctc           >  2:739740/1‑151 (MQ=255)                                                                                                                                            ggTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACg        <  1:1152656/151‑1 (MQ=255)                                                                                                                                            ggTAAAACGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACg        >  2:304938/1‑151 (MQ=255)                                                                                                                                                 aaCGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTTTCACGATAcc   >  2:270258/1‑151 (MQ=255)                                                                                                                                                  aCGCTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCg  <  1:803653/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCG  >  NC_000913/1825906‑1826201

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