Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,826,466	G→A	S153F (TCC→TTC)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,466	0	G	A	100.0%	216.8 / NA	60	S153F (TCC→TTC)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/29);  total (31/29)

GTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTG  >  NC_000913/1826187‑1826751                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                              gtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGt                                                                                                                                                                                                                                                                           >  1:178341/1‑300 (MQ=255) gtTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCAAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGt                                                                                                                                                                                                                                                                           >  1:177683/1‑300 (MQ=255)      cACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATc                                                                                                                                                                                                                                                                      <  1:756238/300‑1 (MQ=255)        cGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTACCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCaa                                                                                                                                                                                                                                                                    >  1:47093/1‑300 (MQ=255)          aTACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTACAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAgg                                                                                                                                                                                                                                                                  >  2:767700/1‑300 (MQ=255)           tACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGt                                                                                                                                                                                                                                                                 <  2:164263/300‑1 (MQ=255)             ccGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGc                                                                                                                                                                                                                                                               >  1:616993/1‑300 (MQ=255)                  acTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAGTTTATAGTCCGAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCAAgcg                                                                                                                                                                                                                                                          >  2:696059/1‑300 (MQ=255)                   ctctCACCAGATAGCAGCGCAGTAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTACAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCAAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAgcgc                                                                                                                                                                                                                                                         >  2:286613/1‑300 (MQ=255)                   ctctCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTCTGTAGATCAAGGTGCCAgcgc                                                                                                                                                                                                                                                         >  2:423488/1‑300 (MQ=255)                       cACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAg                                                                                                                                                                                                                                                     >  1:491500/1‑300 (MQ=255)                       cACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCAAGTCAGAAGTTTCTCGTCTAAGGGGATGTCGCGTTGCGGTAGTACACCGAACTGCGGATACAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAg                                                                                                                                                                                                                                                     >  1:491492/1‑300 (MQ=255)                         ccAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGAt                                                                                                                                                                                                                                                   >  1:232546/1‑300 (MQ=255)                                       agcaATTACGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGTGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGCCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTc                                                                                                                                                                                                                                     <  2:375474/300‑1 (MQ=255)                                                                    ggCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACATGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAATAAGCCACGAATTGCGGGATGGAGGTCAAGGTGCCAGCGCGCAGGTTCTGTGCCCTGGGCGGAAAAATCTTCCATGCGCTGTGCCAggg                                                                                                                                                                                                        <  2:81380/300‑3 (MQ=255)                                                                    ggCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGNTGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGtt                                                                                                                                                                                                        <  1:556648/300‑1 (MQ=255)                                                                    ggCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCATGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGtt                                                                                                                                                                                                        <  2:779965/300‑1 (MQ=255)                                                                    ggCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCGGGTCATAAAAATCTTCCAGGCACTGTTCCAGtt                                                                                                                                                                                                        <  1:157839/300‑1 (MQ=255)                                                                    ggCGAAGATCGTTTTGCCCAAAGGGCAACGCTGTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTTGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAACTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTGGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGtt                                                                                                                                                                                                        <  1:129778/300‑1 (MQ=255)                                                                                                           gTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCt                                                                                                                                                                 >  1:455007/1‑300 (MQ=255)                                                                                                            ttCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTACGCTTTCAAAAAATAGCTg                                                                                                                                                                >  2:67056/1‑300 (MQ=255)                                                                                                              ccAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCATTGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGcc                                                                                                                                                              <  2:63043/300‑1 (MQ=255)                                                                                                                   gTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGCCAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGa                                                                                                                                                         >  1:180521/1‑300 (MQ=255)                                                                                                                   gTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGa                                                                                                                                                         >  2:472721/1‑300 (MQ=255)                                                                                                                        ggCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGcc                                                                                                                                                    >  1:172808/1‑300 (MQ=255)                                                                                                                                   cAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGAt                                                                                                                                         <  1:318286/300‑1 (MQ=255)                                                                                                                                                  tGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTACGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGc                                                                                                                          <  1:317345/300‑1 (MQ=255)                                                                                                                                                   gAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCATGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                                                                                                         <  2:291317/300‑1 (MQ=255)                                                                                                                                                   gAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                                                                                                         >  1:440398/1‑300 (MQ=255)                                                                                                                                                         gttAAACTTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGTCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCATGGAATGTCGCGTTGCGGTAATACACCGAATTGCGGATGCAAGAAGCCNCGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAAATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACg                                                                                                                   <  2:181785/298‑1 (MQ=255)                                                                                                                                                            aaaCGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGttt                                                                                                                <  2:462953/300‑1 (MQ=255)                                                                                                                                                                           ttCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGTGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTACGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGATATGGCAATAACGTTTCCCTTGCTTCAgcgc                                                                                                 >  2:420845/1‑300 (MQ=255)                                                                                                                                                                             ccTGATGGAACACCAGCGCCTACAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCAAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAAAAGCCACGAATTGCGGTATGTAGATCAAGGTGTCAGTGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACCTTCGATTCATATCGCTATGGCAATAACGTTTACCTTGCTTCCGCGCag                                                                                               >  2:242669/1‑300 (MQ=255)                                                                                                                                                                                gATGGAACACCAGCGCCTCCAGCCCCGCCACACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTGTCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCTTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCaggag                                                                                            <  2:46884/300‑1 (MQ=255)                                                                                                                                                                                 aTGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTNCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTAAGCGCaggagg                                                                                           >  1:578736/1‑300 (MQ=255)                                                                                                                                                                                       cacCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTccc                                                                                     >  2:177832/1‑300 (MQ=255)                                                                                                                                                                                             cgcCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGTGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCNGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCATGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt                                                                               <  2:699182/300‑1 (MQ=255)                                                                                                                                                                                               ccTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCAAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCa                                                                             >  2:15533/1‑300 (MQ=255)                                                                                                                                                                                                cTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAc                                                                            >  1:714599/1‑300 (MQ=255)                                                                                                                                                                                                       cccGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAAc                                                                     <  1:276867/300‑1 (MQ=255)                                                                                                                                                                                                                      cACGTCAGAAATTTCTCGTCCCATGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                      <  2:338145/300‑1 (MQ=255)                                                                                                                                                                                                                      cACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACggg                                                      <  1:458205/300‑1 (MQ=255)                                                                                                                                                                                                                                   tctcGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTgag                                         <  1:87435/300‑1 (MQ=255)                                                                                                                                                                                                                                     tcGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGat                                       <  1:683022/300‑1 (MQ=255)                                                                                                                                                                                                                                       gTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGatat                                     <  1:759592/300‑1 (MQ=255)                                                                                                                                                                                                                                       gTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGatat                                     <  2:528210/300‑1 (MQ=255)                                                                                                                                                                                                                                        tCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATc                                    >  1:646992/1‑300 (MQ=255)                                                                                                                                                                                                                                            aGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCgcgc                                <  2:506347/300‑1 (MQ=255)                                                                                                                                                                                                                                            aGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATNCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCgcgc                                >  1:561746/1‑300 (MQ=255)                                                                                                                                                                                                                                             gggAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCAGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGACGTGAGATATCGCGcc                               >  1:273106/1‑300 (MQ=255)                                                                                                                                                                                                                                                 aTGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTNCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGc                           <  1:547269/300‑1 (MQ=255)                                                                                                                                                                                                                                                         ttGCGGTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTc                   >  1:703138/1‑300 (MQ=255)                                                                                                                                                                                                                                                              gTAATCCACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCACCCGCCAACGTAACGGGATCTCGCCGTGAGGTATCGCGCCAAGCAACGCGTCCCGCa              >  1:612094/1‑300 (MQ=255)                                                                                                                                                                                                                                                              gTAATACACCGAACTGCGGATGCAAGAAGCCNCGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGNCAGGTTTCTCCGCTTTCAACAAATAGGTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGGCCAGCa              >  2:561871/1‑300 (MQ=255)                                                                                                                                                                                                                                                              gTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACTGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCa              <  2:363956/300‑1 (MQ=255)                                                                                                                                                                                                                                                              gTAATACACCGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCa              >  1:410422/1‑300 (MQ=255)                                                                                                                                                                                                                                                                       cGAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAg     <  1:538274/300‑1 (MQ=255)                                                                                                                                                                                                                                                                        gAACTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCCGCAACCGCCCACGTAACGGGATCTCGCCGTGAGATATCGCGCAACGCAACGCGTCCAGCATCTACACAgg    >  2:516114/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        gAACTGCGGAGGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAgg    <  2:789624/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          aCTGCGGATGCAAGAAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAACAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTg  <  1:25957/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                              GTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTG  >  NC_000913/1826187‑1826751

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