Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,516,959	C→T	P253P (CCG→CCA)	fecR ←	ferric citrate regulator FecR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,516,959	0	C	T	100.0%	138.3 / NA	46	P253P (CCG→CCA)	fecR	ferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/27);  total (19/27)

CGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTC  >  NC_000913/4516722‑4517205                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         cGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGGCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCGCC‑TCACCCAGTGGTTTATCGCTGa                                                                                                                                                                                                                                             <  2:127854/251‑1 (MQ=255)  gagaAGGGTAAAAAGGAAAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGaa                                                                                                                                                                                                                                            <  2:233446/251‑1 (MQ=255)             aaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGTCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGt                                                                                                                                                                                                                                 <  1:249917/251‑1 (MQ=255)              aaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTc                                                                                                                                                                                                                                <  1:35257/251‑1 (MQ=255)                   aaaTCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATGTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAGCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCg                                                                                                                                                                                                                           <  1:107414/250‑1 (MQ=255)                       cGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCa                                                                                                                                                                                                                       <  1:38541/251‑1 (MQ=255)                        gAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAg                                                                                                                                                                                                                      <  2:5289/251‑1 (MQ=255)                           aTAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGCTCGCATCGGAATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTTCGGTAACGGGTTAGCGTGGCTAT‑ACCATCACCCACTGGTTTATCTCTGTAGCTAAGGATGTCCTTCGTCCAGCtt                                                                                                                                                                                                                   >  2:148714/1‑251 (MQ=255)                             aaGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATCGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCAAGCCGTTGCGGTAACGGGTTAGAGGGGCTATCACC‑TCACCCAGTGGATTATCGCTGAAGCTAAGGATGTCCTTCGTCCAGCTTGt                                                                                                                                                                                                                 >  1:244013/1‑251 (MQ=255)                                 attattTTCCGTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTActc                                                                                                                                                                                                             <  2:94471/251‑1 (MQ=255)                                            tACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCac                                                                                                                                                                                                  >  1:146595/1‑250 (MQ=255)                                              cAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCg                                                                                                                                                                                                <  1:78578/251‑1 (MQ=255)                                                  ggtgAAATGGTTAGCCAGTACCGCTTAATAGACGGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCATCTGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGTGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGGTAGCGTGTCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGttt                                                                                                                                                                                            <  1:183061/251‑1 (MQ=255)                                                                 cAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTc                                                                                                                                                                            >  1:194891/1‑251 (MQ=255)                                                                            aaTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑Tg                                                                                                                                                                 <  2:136418/250‑1 (MQ=255)                                                                                     aaTTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGCTCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGTCGGGATCGCAGCGCCGCCCCCCGTTGCGGTAACGGGTTAGCGTGTCTATCACC‑TCCCCCAGTGGTTTTTCGCTGACGCTCAGGTTGTCCTTCGTCCAGCTTGTCCTCTCCTCATCCCCCGGTTTTACTCCTCC‑AAAcccccatgcgcttaactcccggc                                                                                                                                                        >  2:227978/1‑229 (MQ=255)                                                                                           aaCGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATTGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGTTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCAc                                                                                                                                                 <  1:230785/251‑1 (MQ=255)                                                                                           aaCGGGAAGCGTTTACGCGATACCGTTCAGGATCGCATCGGAATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGAGCAGCACGCCGTTGGGGTAACGGGATAGAGTGGGTATCACC‑TCAACCCGTTGATTATCGATGAAGCTCAGGATGCCCTTCGTCCAGCATGCACTCTCGTCCTCCAGCGGTTTCACTGCGCC‑AAACTAAGAGCCGC‑TGAACTGCCTGCTTTcgc                                                                                                                                                 >  2:192232/1‑249 (MQ=255)                                                                                              gggAAGCGTTTGCGCGATAACGTTCAGGATCGCCTCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAg                                                                                                                                              <  2:146595/251‑1 (MQ=255)                                                                                              gggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGCACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCGGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTCTCCCC‑TCACCCAGTGGTTTa                                                                                                                                                                                                                                                    >  1:121777/1‑151 (MQ=255)                                                                                               ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCTGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGGGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTACAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGc                                                                                                                                             <  1:109410/251‑1 (MQ=255)                                                                                                gAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGAGTGGCTATCAAC‑TCACCCAGTGGTTTATCGCTGAAGCCCAGGATTTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTCCGCC‑AAACTAAGAGGCGC‑AGAACTGCAGGCTTTCAcccg                                                                                                                                              >  1:226426/1‑247 (MQ=255)                                                                                                 aaGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTCTCGCTGAAGCTCAGGGTGTCCTTCGTCCAGCTTTTAATCTCGTCATCCAGCGGTTTCACTGCGCC‑CAACTCAGGGGCGC‑TGAAATGCCGGCTTTCACCAGCGt                                                                                                                                           >  2:44945/1‑251 (MQ=255)                                                                                                   gaggttGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAAAGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAACGGGTTAGCGGGGCTAGCACC‑TCACCCAGTGGTTTATCTCTGAAGCTCAGGATGGCCTTCGTCCAGCGTGAACTCTCGTCACCCAGCGGGTCCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTc                                                                                                                                         <  1:118492/247‑1 (MQ=255)                                                                                                     gTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCAc                                                                                                                                       >  1:225543/1‑251 (MQ=255)                                                                                                       ttGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACg                                                                                                                                      <  1:239175/250‑1 (MQ=255)                                                                                                               tAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGtt                                                                                                                                                                                             >  2:123255/1‑189 (MQ=255)                                                                                                               tAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGtt                                                                                                                                                                                             <  1:123255/189‑1 (MQ=255)                                                                                                                     tCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGcgcg                                                                                                                       >  1:168252/1‑251 (MQ=255)                                                                                                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCAc                                                                                                                                       <  2:113921/228‑1 (MQ=255)                                                                                                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCAc                                                                                                                                       >  1:113921/1‑228 (MQ=255)                                                                                                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTgc                                                                                                                          >  1:68293/1‑241 (MQ=255)                                                                                                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTgc                                                                                                                          <  2:68293/241‑1 (MQ=255)                                                                                                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCa                                                                                                                 <  1:41386/250‑1 (MQ=255)                                                                                                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAGCGCAGCGCAGCACAACGTTGCGGAAACGGGCTAGCGGGGCTCTCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATAGCCTTAGTCCGGCTTGTACTCTCGTCATCCAGCGGTTTCACTTCGCC‑AGCCTCAGAGTAGC‑TGAACGTCCGACTTTCACCGGCGTTCACGATGCGTTATTTcgcgtaggcgc                                                                                                                >  1:148547/1‑244 (MQ=255)                                                                                                                                             gCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCCCC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTAATCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGCGAAGCACTTcc                                                                                               >  1:53419/1‑251 (MQ=255)                                                                                                                                              cGGGAACGTCCCGCTCAGGCGCAGCCCGGCACCGGCGGGATCGCAGCTCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCTGGGGCACTGGCGAGAAGCACTTCca                                                                                              <  1:32518/251‑1 (MQ=255)                                                                                                                                                      tCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCAAGCCCTTGCGGTAACGGGTTCGCATGGCTATAACC‑TTACCCAGTGGTTTATCGCTGAAGCTCCGAATGTCCTTCGTCCAGCTTGTACTATAGCCATCCCACGGTCTCCCTGCGCCCAAACT‑A‑AAGCGCCTGGACCGCAAGATATCCCCACGGGTCACGATCCGTTTTTGCGCGGGGGCCCTGGCGTGAAAATCTTCCACAGCgaac                                                                                      >  1:134902/1‑248 (MQ=255)                                                                                                                                                       cccGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgc                                                                                      <  1:87085/250‑1 (MQ=255)                                                                                                                                                        ccGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGTTTAGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctg                                                                                    <  1:33279/251‑1 (MQ=255)                                                                                                                                                              atgcgcAGCCCGGCAACGGCGGTATCGCAGCGCTGCACGCTGTCTCGGTAACGGGTTAGCTTGGCTATCTCC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATTTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTGCACTGCGCC‑CAACGCAGAGTCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctgctgaac                                                                              <  1:158275/249‑4 (MQ=255)                                                                                                                                                                                                   cGCCGTTGCGGTAACGGTTTAGCGTTGCTATCACC‑TCACCCAGTGGTTTAGCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTGCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTa                                         <  2:58434/251‑1 (MQ=255)                                                                                                                                                                                                               aaCGGGCTCGCGTGGCTATCACC‑TCACCCCGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCAAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCCAGCTGCGTGAAATTATCCTGCTGGCGGCCGGTAAATTCTGTCCCt                             >  2:253633/1‑251 (MQ=255)                                                                                                                                                                                                                       aGCGTGGCTATCACC‑TCACCCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCTCGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTGAAGCAg                      <  1:192232/250‑1 (MQ=255)                                                                                                                                                                                                                                    cc‑tcaccCAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGa        <  2:201306/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccAGTGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTc  >  1:31261/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         CGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC‑TCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCC‑AAACTCAGAGGCGC‑TGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTC  >  NC_000913/4516722‑4517205

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