Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,410,859	C→T	L103F (CTC→TTC)	yjfJ →	PspA family protein YjfJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,859	0	C	T	100.0%	142.9 / NA	45	L103F (CTC→TTC)	yjfJ	PspA family protein YjfJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/19);  total (26/19)

AGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGT  >  NC_000913/4410621‑4411101                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   agaaACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGCTTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTg                                                                                                                                                                                                                                        >  2:149763/1‑251 (MQ=255) agaaACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATCAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGACTTTCATTACCGCTg                                                                                                                                                                                                                                        >  1:378935/1‑251 (MQ=255)        cAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGCTAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGCACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGc                                                                                                                                                                                                                                  >  2:225605/1‑250 (MQ=255)              gTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGtgtt                                                                                                                                                                                                                           >  1:268578/1‑251 (MQ=255)              gTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGtgtt                                                                                                                                                                                                                           >  1:329514/1‑251 (MQ=255)                              aGCATAGTCGTGACGCTAAAGCTGAAATCGATAAAGCCGGAAAATCTCGCGTTGATCTGCGGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAggc                                                                                                                                                                                                           <  2:139630/251‑3 (MQ=255)                                   atTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCg                                                                                                                                                                                                      >  1:194501/1‑251 (MQ=255)                                        tGACGCTAAAGCTGCACTCGATAAAGCCGGCAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGACGCCCGTGCGCTGGAAGCGTTAAGCACGAACGTTAATCCGTCGTTGTTTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGACTTTCCTTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATg                                                                                                                                                                                                 >  2:429923/1‑251 (MQ=255)                                        tGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATCAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGCTTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATg                                                                                                                                                                                                 >  2:214979/1‑251 (MQ=255)                                           cGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTcc                                                                                                                                                                                                        >  2:96517/1‑241 (MQ=255)                                           cGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTcc                                                                                                                                                                                                        <  1:96517/241‑1 (MQ=255)                                            gCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGt                                                                                                                                                                                             >  2:233175/1‑251 (MQ=255)                                               aaaGCTGAACTCGATAAAGCCTGAAAATCGCGCGTTGATCTGATGGCGCGGGTTAAATTAAGTCACTATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTg                                                                                                                                                                                            <  2:54537/249‑1 (MQ=255)                                                 aGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGaa                                                                                                                                                                                         >  1:86783/1‑250 (MQ=255)                                                       aCTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAg                                                                                                                                                                                  <  1:233175/251‑1 (MQ=255)                                                          cGATAAAGCCGGAAAATCTCGAGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGtt                                                                                                                                                                                <  1:89784/250‑1 (MQ=255)                                                            aTAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTaca                                                                                                                                                                             <  1:329875/251‑1 (MQ=255)                                                              aaaGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGGGATGGCGTGGAAAAAGCAGTTacag                                                                                                                                                                            <  1:392349/250‑1 (MQ=255)                                                                 gCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCATGCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTacagcgac                                                                                                                                                                        <  2:52496/251‑1 (MQ=255)                                                                      aaaaTCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGACGCGTTAAGCAAAAACGTTAATCCGTCGTTGATTAACGACGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGCAGAGCAGGTGTTGTCGAATCTGGAAGTCACCCGTGATGGAGAGGAAAAAGCTGTTACAGCGAccgcgc                                                                                                                                                                   >  1:442194/1‑251 (MQ=255)                                                                              gcgTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATc                                                                                                                                                           <  1:283984/251‑1 (MQ=255)                                                                                           ggCGCGGGTGAAATTAAGTCACGATATGCTGTAAGTTTTACGTTAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGCATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGACGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCTCAGCGTATCGCTCAGTTTGagc                                                                                                                                              >  2:318922/1‑251 (MQ=255)                                                                                                 ggTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGGCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGATAATTTCATTACCGCTGAAGAGTAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCaaa                                                                                                                                        <  1:40570/251‑1 (MQ=255)                                                                                                  gTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAAt                                                                                                                                       <  1:341799/251‑1 (MQ=255)                                                                                                               aCGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTaa                                                                                                                           >  1:462351/1‑250 (MQ=255)                                                                                                                                              aGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTACTCCGTCGTTGGTTAACGAAGATGCTGAAGAAATCGCGCGCCATGAGACTTGCATTACCGCTGAAGAGCCAGAGTTGTCGAATCTGGAAGTCTCCCGAGATGGCGTGGAAAAAGCAGTTACAGCGGCAGCGCAGCGTATCGATCAGTTTGAGCAGCAAAAGGAAGTCGTAAAAGCCACTGAAGCCATGCAGAGAGCACAAAAgg                                                                                           >  2:297805/1‑251 (MQ=255)                                                                                                                                                   gTCTGGAAGCCCGTGCGCTGGAAGCGTAAAGCAAGAACGTTAATCCGTCGTTGATTAACGAATTTGCTGAAGAAATCGCGCTCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGTAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTTCAGCGACAGCGCTGCGTATCGCTCAGTTTGAGCATCAAATGGAATTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACATCAGGCTGTa                                                                                      <  1:226744/251‑1 (MQ=255)                                                                                                                                                        gAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaac                                                                                 >  2:365321/1‑251 (MQ=255)                                                                                                                                                                 gcgcTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGGGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCg                                                                         <  2:194501/250‑1 (MQ=255)                                                                                                                                                                                            aaTCCGTCGTTGATTAACGTAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGGAGTCTACCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGATGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTATACCTTTGACGCATATTCCACATTTTCGACAGc                                             >  2:276568/1‑251 (MQ=255)                                                                                                                                                                                               ccatcgtTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGGAACAACCTCTACCGTGGGCGCATCTTCCAGCGTTTCGACAGCGGc                                          <  1:139892/248‑1 (MQ=255)                                                                                                                                                                                                cgtcgtTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCCGCGTTTCGACAGCGGc                                          >  1:476275/1‑250 (MQ=255)                                                                                                                                                                                                    gtTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaac                                     >  2:25284/1‑250 (MQ=255)                                                                                                                                                                                                       gATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATg                                                                                                           <  1:295088/178‑1 (MQ=255)                                                                                                                                                                                                       gATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATg                                                                                                           >  2:295088/1‑178 (MQ=255)                                                                                                                                                                                                              gAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACg                           >  1:312978/1‑251 (MQ=255)                                                                                                                                                                                                               aaGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGc                          >  1:45864/1‑251 (MQ=255)                                                                                                                                                                                                                       gaagaaATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAc                  <  1:175190/251‑1 (MQ=255)                                                                                                                                                                                                                       gaagaaAGCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAc                  <  2:261626/251‑1 (MQ=255)                                                                                                                                                                                                                          gaaATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAcgcg               >  1:47695/1‑251 (MQ=255)                                                                                                                                                                                                                             aTCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCGGTAACAACCTCTACCGTTGGCGCAACTTCCAGCGATTCGACAACGGCAGACTCCTTCAAACGCCTGCAAACGCGTc             >  2:74202/1‑250 (MQ=255)                                                                                                                                                                                                                                  gcgcCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCg       <  2:45864/251‑1 (MQ=255)                                                                                                                                                                                                                                     ccTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                                                          >  2:448537/1‑181 (MQ=255)                                                                                                                                                                                                                                     ccTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                                                          <  1:448537/181‑1 (MQ=255)                                                                                                                                                                                                                                       ttGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTAAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGACGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGt  >  1:433174/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   AGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGT  >  NC_000913/4410621‑4411101

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: ‑