Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,814,697	C→T	F277F (TTC→TTT)	katE →	catalase HPII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,814,697	0	C	T	100.0%	207.6 / NA	62	F277F (TTC→TTT)	katE	catalase HPII
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/30);  total (32/30)

CAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAAC  >  NC_000913/1814466‑1814936                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                cAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCga                                                                                                                                                                                                                                >  2:568319/1‑249 (MQ=255)     aaCACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCgaaggga                                                                                                                                                                                                                           <  1:573946/250‑1 (MQ=255)       cacGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCgaagggaagg                                                                                                                                                                                                                        <  1:367828/251‑1 (MQ=255)        acGCCAATCTTCTTTATCCAGGTTGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGc                                                                                                                                                                                                                       <  1:613377/251‑1 (MQ=255)         cGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGc                                                                                                                                                                                                                       >  2:564002/1‑250 (MQ=255)            cAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACg                                                                                                                                                                                                                   >  1:373662/1‑251 (MQ=255)              aTCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTAATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGtt                                                                                                                                                                                                                 >  2:297698/1‑251 (MQ=255)                  tcttTTTCCAGGATGCCCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGTCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGTCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGttt                                                                                                                                                                                                                <  1:623852/248‑1 (MQ=255)                   cttTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTACAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAATGCCCACGATACTTTCTGGGATTATGTTTCTGTGCAACATGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGAGGCATCCCCCGCAGTTACCGCACAATGGAAGGCGTCGGTATTCACACCTTTCGCCTGATTAATGCGGAAGGGAAGGCAACGTTTGt                                                                                                                                                                                                              >  1:261156/1‑249 (MQ=255)                     ttATCCAGGATGCGCATAAATTCCCTGATGTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTa                                                                                                                                                                                                             <  1:628848/248‑1 (MQ=255)                     ttATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGTGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACg                                                                                                                                                                                                           <  1:438809/250‑1 (MQ=255)                     ttATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATTTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTAc                                                                                                                                                                                                            <  2:631162/249‑1 (MQ=255)                       aTCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGttt                                                                                                                                                                                                        >  1:402509/1‑251 (MQ=255)                         ccAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTcc                                                                                                                                                                                                      >  1:190213/1‑251 (MQ=255)                           aGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTCTTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCAc                                                                                                                                                                                                    >  2:504390/1‑251 (MQ=255)                            ggATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGTCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGAGCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACt                                                                                                                                                                                                   <  1:461698/251‑1 (MQ=255)                            ggATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACt                                                                                                                                                                                                   >  2:347508/1‑251 (MQ=255)                              aTGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTgg                                                                                                                                                                                                 >  1:316533/1‑251 (MQ=255)                                  gcATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGACGGCAACGTTTGTACGTTTCCACTGGAAAc                                                                                                                                                                                             >  2:235235/1‑251 (MQ=255)                                                  tttGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGc                                                                                                                                                                             >  2:166695/1‑251 (MQ=255)                                                         aTGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTc                                                                                                                                                                      >  1:85888/1‑251 (MQ=255)                                                              gTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCCCTGGAAACCACTGGCAGGTAAAGCCCCACTCGtttt                                                                                                                                                                 >  2:51463/1‑250 (MQ=255)                                                                aaaaCCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTggg                                                                                                                                                               <  1:347508/251‑1 (MQ=255)                                                                 aaaCCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCTCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGTCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTGTCCACTGGAATCCACTGGCAGGTAAAGCCTCACTCGTTTGGGa                                                                                                                                                              <  1:583382/251‑1 (MQ=255)                                                                  aaCCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGAt                                                                                                                                                             <  1:568319/251‑1 (MQ=255)                                                                           cGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACaa                                                                                                                                                    <  2:225018/251‑1 (MQ=255)                                                                            gCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCAc                                                                                                                                                                        <  1:75952/230‑1 (MQ=255)                                                                            gCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCAc                                                                                                                                                                        >  2:75952/1‑230 (MQ=255)                                                                               cTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCCCTCGTTTGGGATGAAGCACAAAAAc                                                                                                                                                >  1:607526/1‑251 (MQ=255)                                                                                tGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACt                                                                                                                                               >  2:207829/1‑251 (MQ=255)                                                                                  ggCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCa                                                                                                                                             <  1:49454/251‑1 (MQ=255)                                                                                              aGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCATCATGGAAGGCTTCGGTATTCACTCCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTc                                                                                                                                                                      <  2:639412/214‑1 (MQ=255)                                                                                              aGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCATCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTc                                                                                                                                                                      >  1:639412/1‑214 (MQ=255)                                                                                              aGGGCAAAGAGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCTTTATCGGATCGCGGCATCCCCCGCCGTTACCACGACATATAAGGCTTCGATATGCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCGCAGATGGTACGTTTCCACTGAATACCACTAAAAGGTAAAGCCTCACTCGTTTGGGGGGAAACAGACAAACTCACCGGACGTGaac                                                                                                                                 >  1:610922/1‑249 (MQ=255)                                                                                                          ccACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAcc                                                                                                                     >  1:581437/1‑251 (MQ=255)                                                                                                                        gggATTATGTTCCTCTGCAACCTGAAACTCTGCCCAACGTGATGTGGGGGATGTCGGATCGCGGCACCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCCCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAATCTCACCGGACGTGCCCCGGACTTCCTCCGCCGCGAGTtg                                                                                                          <  2:373662/248‑1 (MQ=255)                                                                                                                              aTGTTTCTCTGCAACCTGACACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGTAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGATGTGAACCGGACTTCCACCGCcgcgccttgagggcagcc                                                                                                 >  1:56488/1‑237 (MQ=255)                                                                                                                                          aaCCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGATGTGGGAAGCCATTGAAGCAggg                                                                                     >  2:47635/1‑250 (MQ=255)                                                                                                                                          aaCCTGAAACTCGGCACAACGTGATGGGGGCGAGGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACCCCTTTCGCCTGATTAAGGCCGAAGGGAAGGCAACGTTGGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGGGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGc                                                                                     <  2:316533/251‑1 (MQ=255)                                                                                                                                                       gCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAAt                                                                        >  1:264141/1‑251 (MQ=255)                                                                                                                                                         acaACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATAc                                                                      <  1:236994/251‑1 (MQ=255)                                                                                                                                                              gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGAt                                                                                                                                                             >  2:627663/1‑159 (MQ=255)                                                                                                                                                              gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGAt                                                                                                                                                             <  1:627663/159‑1 (MQ=255)                                                                                                                                                                gATGTGGGCGATGTCGGATCGCGCCACCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCCCACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCCCTGGCAGGTAAAGCCTCACCCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTg                                                                <  1:58228/250‑1 (MQ=255)                                                                                                                                                                   gtgGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCACCGTTTGTACGTTTCCACTGGAACCCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCt                                                            <  1:408860/251‑1 (MQ=255)                                                                                                                                                                     gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGTTATTCACACCTTTCCCCTGATTAATGCCGAAGGGAAGTCAACTTTGGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCtt                                                           <  2:537674/250‑1 (MQ=255)                                                                                                                                                                          aTGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGtt                                                     <  1:95141/251‑1 (MQ=255)                                                                                                                                                                           tGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTg                                                    >  1:382838/1‑251 (MQ=255)                                                                                                                                                                            gTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACTTTTGTACTTTTCCACTGGAACCCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGa                                                   <  2:190213/251‑1 (MQ=255)                                                                                                                                                                            gTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTg                                                    <  1:47635/250‑1 (MQ=255)                                                                                                                                                                             tCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCCCTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGa                                                   <  1:564002/250‑1 (MQ=255)                                                                                                                                                                                   cgcgGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCCGTTGATTCCTg                                             >  2:343028/1‑250 (MQ=255)                                                                                                                                                                                                    ttACCGCACCATGGAAGGCTTCGGTATTCACACCTTTAGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGAAAGCCTCACTTGTTTGGGATGAAGCACAAAAACTCAACGGAAGTGACCCCGACTTCCACCGCCGCGAGTTGTGGGGAGCCATTGCAGCCAGCGATTTTCCGGCATACGAACTGGGCTTCCCGTTGATTTCTGGAGAAGATGAATTCAAGt                           >  2:485151/1‑251 (MQ=255)                                                                                                                                                                                                      aCCGCACCATGGAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGCAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTc                         >  1:436936/1‑251 (MQ=255)                                                                                                                                                                                                        cGCACCATGGAAGGCTTGGGTATTCACCCCTTTCGCCTGATTAATGCCGAAGGGAAGGAACCGTTTGTACGTTTCCATTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTGGGGATGAAGCACAAAACCTCACCGGACGTGACCAGGACTTCCACCGCCGCGAGGTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCGTCCAGTTGAGTCCTGAAGAAGATGAATTCAAGTTCGa                       <  1:209361/251‑1 (MQ=255)                                                                                                                                                                                                                ggAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtc               >  1:168443/1‑251 (MQ=255)                                                                                                                                                                                                                ggAAGGCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAt                <  1:207829/250‑1 (MQ=255)                                                                                                                                                                                                                     gCTTCGGTATTCACCACTTTCGCCGGATTAAGGCGGAAGGGAAGGAAACGTTTGTGCGTTGCCATGGGACACAATCGGCAGGTAAAGCCGCCGTCTTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGATTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAGTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTtctc          <  2:613375/251‑1 (MQ=255)                                                                                                                                                                                                                     gCTTCGGTATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGGTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTACTGAAGAAGATGAATTCAAGTTCGACTTCGATCTtctc          >  1:467592/1‑251 (MQ=255)                                                                                                                                                                                                                         cgttATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGAt       <  1:604685/247‑1 (MQ=255)                                                                                                                                                                                                                            tATTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATcca    >  1:57324/1‑250 (MQ=255)                                                                                                                                                                                                                             aTTCACACCTTTCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATccaac  >  1:58036/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                CAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCAATTCCACAAGGGCAAAGTGCCCACGATACTTTCTGGGATTATGTTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAAC  >  NC_000913/1814466‑1814936

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