Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,162	C→T	H255Y (CAC→TAC)	spoT →	bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,162	0	C	T	100.0%	243.9 / NA	68	H255Y (CAC→TAC)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/32);  total (36/32)

AGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGT  >  NC_000913/3822918‑3823403                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  aGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGTCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGTTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCg                                                                                                                                                                                                                                              <  2:37984/250‑1 (MQ=255) aGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCg                                                                                                                                                                                                                                              >  1:358972/1‑250 (MQ=255)         accacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGTTCATGGAc                                                                                                                                                                                                                                     >  1:227362/1‑251 (MQ=255)           cacaTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACAt                                                                                                                                                                                                                                   >  1:346798/1‑251 (MQ=255)                   aaCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTCATCAAAGAAGCGGTGACAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCtt                                                                                                                                                                                                                           >  2:491184/1‑251 (MQ=255)                    aCCGAACCCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAGGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGCGATGATCCAGAAGATTCTTTCTGACATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCACAATGGTGCTCAACGAGCCGCGTTTTTACTCGATCATGGACATCTACGCttt                                                                                                                                                                                                                          >  2:147425/1‑251 (MQ=255)                       gAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCcg                                                                                                                                                                                                                       <  1:119527/251‑1 (MQ=255)                              aatatctgtgtgttgagtctctggaTCCCAACCGTGATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCTGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAAGCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGTTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAGGGTGCTCAAAGAGCAGCGTGTTTACTCGATCATGGACAACTACTCTTTCCGCGTGATc                                                                                                                                                                                                                <  2:209464/227‑1 (MQ=255)                              aaGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGTTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTAGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATc                                                                                                                                                                                                                <  1:345544/251‑1 (MQ=255)                              aaGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATGCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCAGCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTGCCGCGTGATc                                                                                                                                                                                                                <  2:227020/251‑1 (MQ=255)                                           tGAGGAGCTGTATCCCAACCGTTATCGCGTAATCAAAGCAGTGGTGAAAGCCGAGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAACGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                                                                                                                                                                                   >  2:490137/1‑251 (MQ=255)                                                    gTATCCCAACCGTTATCGCGTAATCAAAGAAGGGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGGTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGtt                                                                                                                                                                                          <  2:386478/251‑1 (MQ=255)                                                      aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGTGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTAt                                                                                                                                                                                        >  2:196650/1‑251 (MQ=255)                                                          cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgcg                                                                                                                                                                                    >  1:482325/1‑251 (MQ=255)                                                                    cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                                                                                                                                                                           <  1:143672/250‑1 (MQ=255)                                                                                   gtggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTa                                                                                                                                                           <  1:196650/251‑1 (MQ=255)                                                                                      gtgAAAGCCGCGCGAGGCAACCGTACAGAGATGATCCAGACGATTCTTTCTGAAATCGACGGGCGTTTGCAGGAAGCGGGCATACCGTGACGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTCCGCTTTCCGCGTGATCGTCAATGATTCTGCCACCTGTTCTCGCGTGCTGGGCCAGATGCACAGCCTGTACa                                                                                                                                                         >  2:518860/1‑250 (MQ=255)                                                                                       tgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACa                                                                                                                                                         >  1:430825/1‑249 (MQ=255)                                                                                        gAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGc                                                                                                                                                      >  1:488653/1‑251 (MQ=255)                                                                                         aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGcc                                                                                                                                                     >  1:239173/1‑251 (MQ=255)                                                                                                  cgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATTCCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTAGCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCggg                                                                                                                                            <  2:358972/251‑1 (MQ=255)                                                                                                      gCAACCGTAAAGAGATGATCCAGAAGATTATTTCTGAAATCGAAGAGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTAACACCTGTTATCGCGTGCTGGGCCAGATGCACAGACTGTACAAGCCGCGTCCGGGCcgc                                                                                                                                        >  2:495847/1‑251 (MQ=255)                                                                                                              aaaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGTTCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                                                                                <  1:477229/251‑1 (MQ=255)                                                                                                                   gatgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                                            >  1:471877/1‑250 (MQ=255)                                                                                                                     tgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc                                                                                                                         >  2:416227/1‑251 (MQ=255)                                                                                                                       atCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                                       >  2:28357/1‑251 (MQ=255)                                                                                                                          cAGAAGATTCTGTCTGAAATCGAAGGGCGTTTGCCGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCTGCGTTTTTACTCGATCATGGACATCGCCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCAt                                                                                                                    <  2:502202/251‑1 (MQ=255)                                                                                                                            gaagaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTc                                                                                                                  >  1:508651/1‑251 (MQ=255)                                                                                                                              agaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCa                                                                                                                <  1:491184/251‑1 (MQ=255)                                                                                                                                atgctttctGAAATCGAAGGTGTTTTGCAGGAAGCGGTAATCCCGTGCCGCGGCAGTGGTCGCGAGAAGCATCTTTAGTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGAGTTTTTACTCGATCAGGGGCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTACCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCTCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaa                                                                                                              <  2:317213/248‑1 (MQ=255)                                                                                                                                aTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGACAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGAGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaa                                                                                                               >  1:313718/1‑250 (MQ=255)                                                                                                                                   ctttctGAAATCGAAGGGCGTTTGCAGGACGCGGGAATACCGTGCCGCGTCAGTGGACGCGCGAAGCATCTTTATTCGATTTACTGCAACATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCACTGCTTCTGACCCCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTcacaaagc                                                                                                           >  2:236266/1‑244 (MQ=255)                                                                                                                                         gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCCGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCCTTCCAAAAGCGAACgg                                                                                                     >  2:166075/1‑251 (MQ=255)                                                                                                                                          aaaTCGCAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATATACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGATATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATAGCCATTCCAAAAGCGAACGGc                                                                                                    >  2:399466/1‑251 (MQ=255)                                                                                                                                           aaTCGAAGGGCGTTGGCAGGAAGCGGGAATACCGGGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGGGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCt                                                                                                   <  1:166075/251‑1 (MQ=255)                                                                                                                                              cGAAGGGCGTGTGCAGGAAGCGGGAATACCGGGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACGCCTGTGATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGGGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTAt                                                                                                 <  1:498142/250‑1 (MQ=255)                                                                                                                                                    gCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCtt                                                                                          <  1:495626/251‑1 (MQ=255)                                                                                                                                                     cgtgtGCAGGAAGCGGGAATACCGTGCCGCGTCAGCGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGGGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCTCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGTGCCGCGTGAAAGACTATAGCGCCATTCCAAAAGCGAACGGCTATCAGTCttt                                                                                         <  1:61521/247‑1 (MQ=255)                                                                                                                                                          gCAGGAAGCGGGAATACCGGGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTCCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCCCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcac                                                                                     <  2:346798/250‑1 (MQ=255)                                                                                                                                                            aGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                                  >  1:19003/1‑251 (MQ=255)                                                                                                                                                                  cGGGAATACCGGGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatg                                                                            <  1:85775/251‑1 (MQ=255)                                                                                                                                                                    ggAATACCGTGCCGCGTCAGTTGTCGCGAGAAGCATCTTTATTCGATTTACTGAAAAATGGTGCTCAAAGAGCAGCGTTTTTACGCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGGCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatgat                                                                          <  2:287140/251‑1 (MQ=255)                                                                                                                                                                       aTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCgg                                                                       <  1:425419/251‑1 (MQ=255)                                                                                                                                                                         aCCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACGCCTCGATGATCGGcc                                                                     <  1:380230/251‑1 (MQ=255)                                                                                                                                                                          ccgggcggcgTCAGTGGTCGCGAGAATCATCTTTATTCGATTTACTTCAAAATGGAGCTCAAAGAGCAGCGATTTTACTCGAACATGGACATCTCCGCTTTCCGCGTGATCGTTAATGATTCTGACACCTTTGATCGCGTGATGGGCCAGATGCACCGCCTGTTCAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCTAACGGCTATCAGTCTTTGCACACCTCGATGGTCGGccc                                                                    <  2:139777/244‑1 (MQ=255)                                                                                                                                                                            gTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGCGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAACAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGc                                                                  >  2:65857/1‑251 (MQ=255)                                                                                                                                                                            gTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCCAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGc                                                                  >  2:255893/1‑251 (MQ=255)                                                                                                                                                                                         ggacgcgAGAAGCATCTTTATTCGCTTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCACGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGCGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAAGGTGTGCCgtt                                                     >  2:119298/4‑249 (MQ=255)                                                                                                                                                                                         ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                                     <  1:112290/251‑1 (MQ=255)                                                                                                                                                                                             gcgAGAAGCATCTTTAGTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCGGTACAAGCCGCGGCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACATCTCGATGATTGGCCCGCACGGTGTGCCGGTTGAg                                                 <  1:518860/251‑1 (MQ=255)                                                                                                                                                                                                   aGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAg                                           <  1:246155/251‑1 (MQ=255)                                                                                                                                                                                                     cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAATGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGAt                                         <  1:127863/251‑1 (MQ=255)                                                                                                                                                                                                               cGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGaa                               <  1:495847/251‑1 (MQ=255)                                                                                                                                                                                                                     aCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGACCGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGCAGATATg                         >  2:93636/1‑251 (MQ=255)                                                                                                                                                                                                                        gCAAAATGGTGCTCAAAGAGCAGGGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGGGATCGTCAATGATGCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                                            <  2:18132/149‑1 (MQ=255)                                                                                                                                                                                                                        gCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                                            >  1:18132/1‑149 (MQ=255)                                                                                                                                                                                                                          aaaaTGGTGCTCAAAGAGCAGCGTTTTTTCTCGTTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACCCCTGTTAGCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCa                    <  2:488653/251‑1 (MQ=255)                                                                                                                                                                                                                          aaaaTGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc                     >  2:15604/1‑250 (MQ=255)                                                                                                                                                                                                                          aaaaTGGAGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc                     >  2:33159/1‑250 (MQ=255)                                                                                                                                                                                                                                gTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGc              >  1:337956/1‑251 (MQ=255)                                                                                                                                                                                                                                 tGCTCAAAGAGCAGCGTTTTTACTCGCTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCCTTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCg             >  2:85119/1‑251 (MQ=255)                                                                                                                                                                                                                                 tGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                                          >  2:503142/1‑158 (MQ=255)                                                                                                                                                                                                                                 tGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                                          <  1:503142/158‑1 (MQ=255)                                                                                                                                                                                                                                  gCTCAAAGAGCAGGGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATAGCGCGATTCCAAAAGCGAACGGCTATCAGTCTGTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCgg            <  1:236266/251‑1 (MQ=255)                                                                                                                                                                                                                                        agagCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATgg      >  1:256515/1‑251 (MQ=255)                                                                                                                                                                                                                                         gagCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACATGTTATCGCATGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGAGCCGGTTGCGGTCCAGATCCGTACCGAAGATACGGACCAGATGGAGGAGATggg     >  2:23055/1‑251 (MQ=255)                                                                                                                                                                                                                                          agcagcGTTTTTACTAGATAATGGCCATATACGCTTTCCGCGTGATCGTCAATGATTCTGACACCAGTTATAGCGTGCTGGGCCCGATGAACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTACAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATgggg    >  2:523379/1‑250 (MQ=255)                                                                                                                                                                                                                                            cagcGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGgtgt  >  2:491664/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  AGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGT  >  NC_000913/3822918‑3823403

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