Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,697	C→T	F334F (TTC→TTT)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,637,697	0	C	T	94.3%	87.7 / ‑4.4	35	F334F (TTC→TTT)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (13/20);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACG  >  NC_000913/4637452‑4637938                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  aaaaaCTATATTGAGCAGTATGTTTATACGTTAACCCATTAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGTGGAAACTTTACGCGAAGGTACGCCTCCGGAAGTGGTGGCTCGTTTTACCTACAACATTCTGACTCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAATTCGTTCTGACTGTTGTTGATGTGTCGGCATTATTgcgccg                                                                                                                                                                                                                                              <  2:265707/251‑1 (MQ=255)  aaaaCTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTcgccgc                                                                                                                                                                                                                                             <  1:189215/251‑1 (MQ=255)     aCTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGGGAATCATAAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTc                                                                                                                                                                                                                                          >  1:377362/1‑251 (MQ=255)             gAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAg                                                                                                                                                                                                                                   <  1:247743/250‑1 (MQ=255)                      gTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCgg                                                                                                                                                                                                                         >  1:211660/1‑250 (MQ=255)                         tatgcgttAACCCATGAGCTAAAAAGCCCACTGGCGGCGATGCGTGTCGCGGCGGAAATTTTTCGCGAAGGTCCGCCACCTGAAGTGGTGGCTCGTTTTACCTACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGGTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTGGATGTGGCTGCATTATTTCGCCGCGTCAGCGaggctctatccgtgc                                                                                                                                                                                                                       <  2:157051/244‑15 (MQ=255)                         tACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGTAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGGAGAGGGGGCGCGTTTTACTGGCAACATTCTGACGCAAAAAGCGCGTGTGCAGGCGTTGGTAGGAACCGTACTACGCCAGGCCAGACTGGGGAACCGTCAGGAAGTCGTTCTGACTGCTGATGATCTGGCAGCCTTTTTTCGCCGCGTCAGGGATGCGCGCAcccggcc                                                                                                                                                                                                                      >  2:59037/1‑246 (MQ=255)                           tgcgttAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATAGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTTGCCGCGTCAGCGAAGCGCGCACCGTGCAGt                                                                                                                                                                                                                    >  1:247325/5‑251 (MQ=255)                                           cTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGTCTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGGTGGCAGaaaaaaacat                                                                                                                                                                                                    >  1:236118/1‑248 (MQ=255)                                             aaaaaGCCCACGTGCGGTGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGGCCGCCACCGGTAGTGGTGGCGCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGaaaaaaacat                                                                                                                                                                                                    <  2:120724/249‑4 (MQ=255)                                                aagcccgcggggggcgATTCGGGGCGCGGCGGAAATTTTACGCGAATGTCCGCCACTGGAAGTGGTGGCTCGTTTTACCGACTACATTCCGACGCAAAAGTCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGATTTCGTTCGGACTGTTGTTGATGGGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCCCCGGTCAGTTGGCAGAAAAAAACATCActc                                                                                                                                                                                               <  2:363832/239‑1 (MQ=255)                                                aaGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGCGGTGGCTCGTTTTACCAACAACATTCTGACGCAAAATGCGCGTCTGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTCTTTCGCCGCGTCAGCGAAACGCGCACCGTGCAGTTGGCAGAAAAAAACATCActc                                                                                                                                                                                               >  1:54977/1‑251 (MQ=255)                                                      ccTGGCGGCGATTCGTGGCGCAGCGGAAACTTTACGCGAAGGTCCGACACCGGAAGTGGTGGCTCGGTTTACCGACAACAGTCTGACGCAAAATGCGTGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGATTCGTCAGGAAGCCGTTCTGACTGTTGTTGATGTGGCGGCAGTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCTGTTGGCAGaaacatactgcacaatatagt                                                                                                                                                                                         <  2:377362/250‑19 (MQ=255)                                                         ggcggcGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGtt                                                                                                                                                                                       <  1:354095/250‑1 (MQ=255)                                                                    gTGGAGCTGCGGAAATTTTACTCGAATGTACGCCACCGGAGTTTGTGGCTCGTTTTACCGACAACTTTCTGACGCAAAAGGCGCGTATGCAGGCATTGGTAGAAACGTTTCTACGCCAGGCAAGACTGTAGAATCGTCAGGATGTCGTTCTGACTGTTGTTGATTTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTTCATGTTTTGCCCACTGAg                                                                                                                                                                           <  2:355470/251‑1 (MQ=255)                                                                         gcggcgGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGGTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTa                                                                                                                                                                       >  1:186163/1‑250 (MQ=255)                                                                            gcggAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTaatgt                                                                                                                                                                   <  1:348469/251‑4 (MQ=255)                                                                                 aaTTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTtgctg                                                                                                                                                              <  1:122151/251‑1 (MQ=255)                                                                                                 ccgccaccgGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACt                                                                                                                                              <  2:142323/251‑1 (MQ=255)                                                                                                   gccaccgGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTgg                                                                                                                                            >  1:274589/5‑251 (MQ=255)                                                                                                           aaGTGGTGTCTCGGTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGGTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGgcg                                                                                                                                    <  1:171859/251‑1 (MQ=255)                                                                                                           aaGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGgcg                                                                                                                                    >  1:276361/1‑251 (MQ=255)                                                                                                               ggtggCTCGTTTTACTGACAACATTCTGCCGCAAAATGCGCGGTTGCAGGCATTGGTAGAAACGTTACTATGCCATGCAAGACGGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGATGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACGCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTTCTTGAGCAGGCGCTgg                                                                                                                                <  2:198925/251‑1 (MQ=255)                                                                                                                      agTTTTACCGACAACATTCAGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGTAAGACCGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGTGAAGCGCTCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACCGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAAttt                                                                                                                         <  1:371520/250‑1 (MQ=255)                                                                                                                                     aTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGTAGCAGGCGCTGGGGAATTTACTGGATAACGCCAt                                                                                                          <  1:120509/251‑1 (MQ=255)                                                                                                                                       tCTGCCGCAAAATGCGCGGATGCAGGTATTGGCTGAAACGTTCCTACGCTATGCAAGACTGGAGACTCTTCATGAAGTTGTTCTGACTGCTGCTGATGTGGAGGCATTATCCCGCCGCGTCTGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAACCACTCTGCAGGTTTCGCCCACCGAGGTTAACGCTGCTGCTGAACCTGCGTTACTGGCTCAGGCGCTGGGAAATTTAATGGATAACGTCATc                                                                                                         <  2:99140/250‑1 (MQ=255)                                                                                                                                        cTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGCCTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTTCCCGCGTCTGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTCCATGTTATGCCAACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTCCTGGTTAACCCCATCg                                                                                                        >  1:22941/1‑250 (MQ=255)                                                                                                                                                aaaTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAccc                                                                                               >  1:47727/1‑251 (MQ=255)                                                                                                                                                  aTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAccccc                                                                                             <  1:3141/251‑1 (MQ=255)                                                                                                                                                             cAGGCATTGGGAGCAACGTTCCTACTCCAGGCAAGCCTGTAGAATCGTCCGGAAGTCTTTCTGCCTGTTGTTGATGTGGCGGCATTAGTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTGCTTGATAACGCCATCGATTTTACCCCCGAGATCGGtc                                                                                   <  2:315553/250‑2 (MQ=255)                                                                                                                                                                gCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAt                                                                               <  1:209287/251‑1 (MQ=255)                                                                                                                                                                      gTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCt                                                                         <  1:197893/251‑1 (MQ=255)                                                                                                                                                                                     cGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGa                                                          <  1:41988/251‑1 (MQ=255)                                                                                                                                                                                       ccAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTgggtc                                                        >  1:104248/1‑248 (MQ=255)                                                                                                                                                                                                                          tGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGt                     >  1:202727/1‑251 (MQ=255)                                                                                                                                                                                                                                             ggcATTATTTCGCCGCGGCAGCGAATCGATCACCGTGCAGTTGGCAGAAAAAAATATCCCTTTACATGTTATGCCCTCTGAGGTTAATTTTGATGCTGAACCGGCTTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCTCCCTTAAGGTGCTGTAGACCGGTAGTGGTATTCCTGACTAcg  <  2:63590/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  AAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACG  >  NC_000913/4637452‑4637938

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