Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,293	T→C	H41R (CAC→CGC)	fecD ←	ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,512,293	0	T	C	100.0%	131.2 / NA	47	H41R (CAC→CGC)	fecD	ferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/22);  total (25/22)

CGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGC  >  NC_000913/4512052‑4512539                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       cGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCCGTCGGTACTCCCTCAATACATAATAATGCTCGCGTCCGGccc                                                                                                                                                                                                                                               >  2:239488/1‑250 (MQ=255)  cccATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGAAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTg                                                                                                                                                                                                                                              <  2:60553/250‑1 (MQ=255)        ccgcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTc                                                                                                                                                                                                                                        >  1:67757/1‑251 (MQ=255)         cgcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCg                                                                                                                                                                                                                                       <  1:265706/251‑1 (MQ=255)           cccgcAAAGGCCAGCAGAGGACGCACACGCCCGGGCAGTGACGGAAAAAGAAGTAGAGCCCCCCCAGAGGCCGGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCGGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGt                                                                                                                                                                                                                                     <  2:67757/251‑1 (MQ=255)            ccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCCGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAACCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAACACATAATAATGCTCCCGTCCGGCCTTCCAGTCGGTc                                                                                                                                                                                                                                    >  2:87561/1‑251 (MQ=255)                             ggCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacgg                                                                                                                                                                                                                   >  1:149937/1‑251 (MQ=255)                              gcagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGACTATCCGGTAATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggc                                                                                                                                                                                                                  >  1:118847/1‑251 (MQ=255)                                 gcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcacg                                                                                                                                                                                                               >  1:15786/1‑251 (MQ=255)                                  caCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcacgg                                                                                                                                                                                                              >  1:214652/1‑251 (MQ=255)                                      aTCACGGGCAGTGACGGCATAAGCAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACCTAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCCGCAGCGCGCGCCACGGCACgggggt                                                                                                                                                                                                          >  2:162916/1‑249 (MQ=255)                                              cAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCa                                                                                                                                                                                                                                           >  2:212048/1‑210 (MQ=255)                                              cAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCa                                                                                                                                                                                                                                           <  1:212048/210‑1 (MQ=255)                                                    cGGCATAAGAGGTAGAGCCCCCACAGAGGCCAGGCGGTCGGCATGGTTAACACCGGGAATACCCGGTGATGCCAGAGGTTTGCGCACAAGCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCGCGCGTCCGGCCTGCCAGTCGGGCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGa                                                                                                                                                                                            <  2:43152/251‑1 (MQ=255)                                                          aagaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATa                                                                                                                                                                                      <  1:114583/251‑1 (MQ=255)                                                           agaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATa                                                                                                                                                                                      >  1:42213/1‑250 (MQ=255)                                                                                            gcATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaa                                                                                                                                                    <  1:132997/251‑1 (MQ=255)                                                                                             cATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaa                                                                                                                                                   <  1:162916/251‑1 (MQ=255)                                                                                                 gtaaaccccGAGATTACCCGGTGATCCCTGCGGGTTGCGCCAAACCCCCTGTCCCAGCACCCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATa                                                                                                                                                <  2:218689/243‑1 (MQ=255)                                                                                                   tAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                                                             >  1:157538/1‑251 (MQ=255)                                                                                                             aaTATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAttt                                                                                                                                   >  1:66426/1‑251 (MQ=255)                                                                                                              atatCCGGTGAGCCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCATCGCAAtttt                                                                                                                                  <  2:42213/251‑1 (MQ=255)                                                                                                               tatCCGGTGATGCCAGAGGGTTGCGCACACTCCCCTGTATCAGCACGCCCGCCCCGGCGACGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTCCTCCATCCATCCATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCCCGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGCGCTGAAAATAACCCACGCACTTTTc                                                                                                                                 >  2:269351/1‑251 (MQ=255)                                                                                                                   cGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATACTAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATcc                                                                                                                             >  2:196482/1‑251 (MQ=255)                                                                                                                    ggTGATGCCAGAGGGTTGCGCACCATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCt                                                                                                                            >  2:174670/1‑251 (MQ=255)                                                                                                                           ccAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAATGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACCTAATAATGCTCGCGTCAGGACTGCAAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGCTCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAACAACGCAATTGTCATCCTCGcctcc                                                                                                                     >  1:117079/1‑251 (MQ=255)                                                                                                                             agagGGTTGCGCCCACTCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGCAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                                                                                   >  1:204312/1‑251 (MQ=255)                                                                                                                             agagGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGCTCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGCAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                                                                                   >  2:226145/1‑251 (MQ=255)                                                                                                                               agGGTTGCGCACAATCCCCTGTACCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                                                                 <  1:124243/251‑1 (MQ=255)                                                                                                                                        cacaAGCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGAc                                                                                                        <  1:174670/251‑1 (MQ=255)                                                                                                                                                   gTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGc                                                                                             >  1:180205/1‑251 (MQ=255)                                                                                                                                                            cgcccgccACGGCGAGGGCTGCACCGCAAAACAGTGCCAGCAGCAAGAGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATACCAGCGCACAGCCCGCAAGGGCAAGGGCGATGAAAATAACCAGCGCAATTTTCATCATCAACTCCTCAAAACCCAGCCAAAGAAAGGGCTGCCAATc                                                                                     >  2:265093/1‑250 (MQ=255)                                                                                                                                                            cgcccgcaACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGAAAGCGCCGCAGTCGGTACTCCATCAATACATAATAATGCCCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACCGCGCACAGCCCGCCAGGGCAAGGGGGATGAAAATAACAAGCGCAATTTTCCTCCTCGCCTCCCCACAAGACAGCCAAAGCAAGGGCTGCCCATcc                                                                                    >  2:182766/1‑250 (MQ=255)                                                                                                                                                                      ggCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCAc                                                                          >  1:203689/1‑251 (MQ=255)                                                                                                                                                                                gCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCg                                                                <  2:148243/251‑1 (MQ=255)                                                                                                                                                                                                      gcaATCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGACAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCTCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCa                                          <  1:55537/251‑1 (MQ=255)                                                                                                                                                                                                      gcaAGCGCGGCATTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCGGCCAATCAGCGCCAGc                                                                            <  2:184426/217‑1 (MQ=255)                                                                                                                                                                                                      gcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGc                                                                            >  1:184426/1‑217 (MQ=255)                                                                                                                                                                                                              ggCAGTCGGTACTCCATCAATCCATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCACGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATACCCAGCGCAATTTTCATCCTCGCCTCCTCACACGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGt                                  >  2:125286/1‑251 (MQ=255)                                                                                                                                                                                                              ggCAGGCGGTACTCCATCAATCCATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGTGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGTGAAGGCCAgttcgcgt                                  <  2:247829/251‑8 (MQ=255)                                                                                                                                                                                                                           ccATCAATACATAATAATGCTCGCGCCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAATCAAGGGCTGCCAATCAGCGCCAGCACTTCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCt                      <  1:179774/250‑1 (MQ=255)                                                                                                                                                                                                                           ccATCAACACATAATAACGCTCGCGTCCGGCCTGCCAGTCGGTCACCAGCGCGCGCCACGGCACGGGGATCACCCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGTGGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCTCAAGCCAGACAAAGCAAGGGCTGACCATCAGCGCCAGCACTGCGCCTGCGTGCGGATCTCCGGGGAAGGCCAGTGCGCGTGTGAGTACATCt                      <  2:192365/250‑1 (MQ=255)                                                                                                                                                                                                                               cAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCca                  <  1:101267/250‑1 (MQ=255)                                                                                                                                                                                                                                 aTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCcagca               <  1:239488/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCATAATAATGCTCGCGTCCGGCCTGCCCGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCagcagc             <  1:197489/251‑1 (MQ=255)                                                                                                                                                                                                                                        ataatGCTCGCGTCCGGCCGGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCTCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCcatcagcatcag        <  2:251109/251‑1 (MQ=255)                                                                                                                                                                                                                                              cTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGc  >  1:210442/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       CGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGC  >  NC_000913/4512052‑4512539

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