Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,137	A→G	V138A (GTA→GCA)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,137	0	A	G	100.0%	259.6 / NA	74	V138A (GTA→GCA)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (46/28);  total (46/28)

AGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGC  >  NC_000913/3815893‑3816355                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                           cgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCCGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGAATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAgcc                                                                                                                                                                                                                      >  1:83033/2‑250 (MQ=255)    caacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGcc                                                                                                                                                                                                                    >  1:254615/1‑250 (MQ=255)      acaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCCTTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCAc                                                                                                                                                                                                                  >  2:356388/1‑250 (MQ=255)             gAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAc                                                                                                                                                                                                          <  1:262726/251‑1 (MQ=255)               gTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                                                                                        <  2:306623/251‑1 (MQ=255)               gTAGCCCTTCATGGGTGAACCGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCACGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                                                                                        >  2:222902/1‑251 (MQ=255)                 aGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGt                                                                                                                                                                                                      <  1:211593/251‑1 (MQ=255)                    tctTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCTCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAAt                                                                                                                                                                                                   >  1:136764/1‑251 (MQ=255)                     ctTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTaa                                                                                                                                                                                                    <  2:358316/249‑1 (MQ=255)                      ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCGACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCg                                                                                                                                                                                                 >  1:164877/1‑251 (MQ=255)                      ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAAt                                                                                                                                                                                                   <  1:50134/249‑1 (MQ=255)                      ttCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCCCACAGGCACCCggcat                                                                                                                                                                                                   >  2:50134/1‑245 (MQ=255)                         aTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGaag                                                                                                                                                                                              <  1:213632/251‑3 (MQ=255)                            ggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCgg                                                                                                                                                                                           >  2:354385/1‑251 (MQ=255)                                         ccTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACCATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGaagc                                                                                                                                                                                             >  1:41527/1‑233 (MQ=255)                                         ccTTCTGCCGTCCCCGTCACTGCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCGGTAGCGCATCTGCCATCGCCACACAGGCACCCGTatgcgaagc                                                                                                                                                                                             <  2:41527/236‑9 (MQ=255)                                                            cTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCAAAGCGGTACGCGTGCCACCATCAGCCTGACGCACGTCg                                                                                                                                                           >  2:207745/1‑251 (MQ=255)                                                                   gatgCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCAGCAGCCTGAAGCACGTCGCAGTCCa                                                                                                                                                    <  2:438059/251‑1 (MQ=255)                                                                     tgCGCCCTTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCTTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAg                                                                                                                                                   <  1:29982/250‑1 (MQ=255)                                                                        gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                                                                                                >  2:90046/1‑250 (MQ=255)                                                                        gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTa                                                                                                                                               <  1:180941/251‑1 (MQ=255)                                                                              cTTCGGTCATCCCTAAGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAACATCCCTTTCATTGGATTGGTTTTCCGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCCCCATCAGCCCGACGCACGTCCCCGTCCAGCGTAATGGTg                                                                                                                                         >  1:359794/1‑251 (MQ=255)                                                                                 cGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAc                                                                                                                                      <  1:32142/251‑1 (MQ=255)                                                                                       tcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGCTCGCAAATTGCTTCACCGTTTACAATAACGACAGAAACCGCGGCTACCATCCCTTTCCTTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                >  1:431867/1‑251 (MQ=255)                                                                                         aCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                >  1:98349/1‑249 (MQ=255)                                                                                              gTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCTCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCtt                                                                                                                         <  2:39345/251‑1 (MQ=255)                                                                                                      cGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCTACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTTCCACCATCAGCCTGAAGCACGTCGCAGGCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATc                                                                                                                  <  2:261527/250‑1 (MQ=255)                                                                                                              gCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACGCACCCAGCGCTTTCAAATCTACTGCcg                                                                                                          <  2:43264/250‑1 (MQ=255)                                                                                                                   tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTCATGGTGAACTCACCCAGCGCTTTCAAATCTCCTGCCGCGCGa                                                                                                     >  2:75994/2‑250 (MQ=255)                                                                                                                   tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                                                                                                >  1:41135/2‑175 (MQ=255)                                                                                                                   tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                                                                                                <  2:41135/174‑1 (MQ=255)                                                                                                                   tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAACGCTTTCAAATCTACTGCCGCGCGa                                                                                                     >  1:261153/2‑250 (MQ=255)                                                                                                                      gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAagtg                                                                                                 >  2:85361/1‑249 (MQ=255)                                                                                                                             aCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGCAGTGCACGGGc                                                                                          >  2:200923/1‑251 (MQ=255)                                                                                                                                tATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                       >  1:243691/1‑251 (MQ=255)                                                                                                                                      aGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCCGCGGCTACCATCCCTTTCATTGGATTGGTTGTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGatg                                                                                 <  2:302192/251‑3 (MQ=255)                                                                                                                                        aTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCt                                                                               <  1:396126/251‑1 (MQ=255)                                                                                                                                                  gCTTCACCGTTTACAATACAGACAGAAACCCCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCAGTTTTCCAACAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAt                                                                     >  1:159969/1‑251 (MQ=255)                                                                                                                                                   cTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAtt                                                                    >  1:151756/1‑251 (MQ=255)                                                                                                                                                      caccGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATtgtg                                                                 <  1:184819/249‑1 (MQ=255)                                                                                                                                                       accGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTATAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGAGGTACGCGTGCCACCATCGGCCTGACGCACGTCGCAGTCCAGAGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCCAAGTGCACGGTCGAACAGCCGCTGGTTTTCCAttacgc                                                                >  1:430867/2‑247 (MQ=255)                                                                                                                                                             tacaATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCCGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAc                                                          >  2:352252/2‑251 (MQ=255)                                                                                                                                                              acaATACCGACAGAAACCGCGGTTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCATACGCTGGATTTCCATTGTGCGTCCAcc                                                         >  1:139371/1‑251 (MQ=255)                                                                                                                                                              acaATACCGACAGAAACCGAGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGCGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAcc                                                         >  2:49218/1‑251 (MQ=255)                                                                                                                                                               caATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCCCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAccc                                                        >  2:193784/1‑251 (MQ=255)                                                                                                                                                                 ataccGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGGTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTg                                                      <  2:90908/248‑1 (MQ=255)                                                                                                                                                                      gACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTAAGCTTGCCGTTTTCCACCAGCTTATGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAATCACACAGCGCTTTCAAATCTACTGCCGCGAGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCttc                                                  >  1:279953/1‑249 (MQ=255)                                                                                                                                                                       aCAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                <  1:85361/251‑1 (MQ=255)                                                                                                                                                                       aCAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                >  1:51288/1‑251 (MQ=255)                                                                                                                                                                         aGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCTCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCtt                                              <  2:167034/251‑1 (MQ=255)                                                                                                                                                                          gAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCttt                                             >  2:301654/1‑251 (MQ=255)                                                                                                                                                                                  ggCTACCATCCCTTTCATTGGATGGGTTTTCAGCTTTCCGTTTTCCACCAGCATCTGTAGCCCATATGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCTAGCGTAATGGAGAACTCACCCAGCGCTTTCAAATCTACCGCTGCACGAGGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTACACCCTGCTTACCTTTCGCTGCtt                                     >  2:130871/1‑251 (MQ=255)                                                                                                                                                                                   gCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCACGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTc                                    >  1:399401/1‑251 (MQ=255)                                                                                                                                                                                      aCCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTTTGTAGCGCATCTGCCAGCGCCACTCAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGCAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACg                                 >  2:59546/1‑251 (MQ=255)                                                                                                                                                                                      aCCATCCCTTTCATTGGATTGGTGTTCAGCGTGCCGTTTTCCTCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGGGCCACCATCAGCCTGAAGCACGTCGCAGTCCGGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCCCGAAGGGCACGTGCGATCAGACGCTGGATGGCCATTGTGCGTCCACCCTGCTTTCCTTTCGCTGCATCACg                                 <  2:431867/251‑1 (MQ=255)                                                                                                                                                                                        cATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACg                                 >  2:134325/1‑249 (MQ=255)                                                                                                                                                                                         aTCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGc                              >  1:208247/1‑251 (MQ=255)                                                                                                                                                                                               atcattgcaTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACCCGCGCCACCATAAGCCTGAAGCACGTCGCAGTCCCGCGTATTGGTGAACTCAACCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCGTTCTGCGTCCACACTGCTTACCTTCCGCTGCTCCACGCGCGTTACg                        >  2:275140/9‑251 (MQ=255)                                                                                                                                                                                                 cattggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACgg                       <  1:356388/246‑1 (MQ=255)                                                                                                                                                                                                 cattggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGCATCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACTTCGCAGTCCCGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCTCACGTAGTGCACGGGCGCTCAGACGCTGGCTTTCCATTGGGCGTCCACCCTGCTTACCTTTCGCTGCTTCCCGATCGTTACggg                      >  2:232991/5‑251 (MQ=255)                                                                                                                                                                                                  attggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGAATCTGCCAGCGCCAAACATGCACACGTAATCGAAGCAGTACGCGTGCCACCATCCGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACATTTCGCTGCTACACGAGCGTTACGGgt                     >  2:182279/4‑251 (MQ=255)                                                                                                                                                                                                     ggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGtgtg                  <  1:67286/251‑1 (MQ=255)                                                                                                                                                                                                      gATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTgg                 >  1:289358/1‑251 (MQ=255)                                                                                                                                                                                                        ttGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGt               >  1:42037/1‑251 (MQ=255)                                                                                                                                                                                                          ggTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCCGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGt               >  2:278577/1‑249 (MQ=255)                                                                                                                                                                                                              ttCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTACTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAAc          <  2:419077/250‑1 (MQ=255)                                                                                                                                                                                                                cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                       >  1:188296/1‑171 (MQ=255)                                                                                                                                                                                                                cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                       <  2:188296/171‑1 (MQ=255)                                                                                                                                                                                                                cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcg       >  1:176700/1‑249 (MQ=255)                                                                                                                                                                                                                 aGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcgg      >  1:127309/1‑248 (MQ=255)                                                                                                                                                                                                                  gcgtGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCTGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTTAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCGCCCTGCTTACCGTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGatcgcggc     <  2:210769/248‑8 (MQ=255)                                                                                                                                                                                                                     tGCCGTTTTCCACCATCTTCGGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAcc                                                                                                                                 <  1:261137/124‑1 (MQ=255)                                                                                                                                                                                                                     tGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGgcagc  >  1:190173/1‑251 (MQ=255)                                                                                                                                                                                                                     tGCCGTTTTCCACCAGATTCTGTAGCGCATCTGCCAGCGCCACGACGGCACCCGTCATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAcc                                                                                                                                 >  2:261137/1‑124 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                           AGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGC  >  NC_000913/3815893‑3816355

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