Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,546,601	C→T	I502I (ATC→ATT)	fimD →	type I fimbriae usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,546,601	0	C	T	97.4%	130.9 / ‑3.8	39	I502I (ATC→ATT)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/14);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTG  >  NC_000913/4546351‑4546838                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              gTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGG‑CACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCTCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGCTACAACATACAGTCGAATGAATGGCTACAACATCg                                                                                                                                                                                                                                              >  1:45312/1‑251 (MQ=255)         gCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAacac                                                                                                                                                                                                                                        >  2:340574/1‑250 (MQ=255)              tAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACg                                                                                                                                                                                                                                  >  2:104108/1‑251 (MQ=255)                            ggAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGGGGGTTACCGTTATTCGCCCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGACTGGCTACAACATTGAAAAAAAGGACGGAGTTATTAAGGtt                                                                                                                                                                                                                    >  2:253791/1‑251 (MQ=255)                                           cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAgg                                                                                                                                                                                                                      <  2:322893/234‑1 (MQ=255)                                           cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAgg                                                                                                                                                                                                                      >  1:322893/1‑234 (MQ=255)                                               gggCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                                                                                                                                                                                                 >  2:111973/1‑251 (MQ=255)                                                   gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTa                                                                                                                                                                                             <  1:111973/251‑1 (MQ=255)                                                                    gACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                                                                                                                                                                             >  1:221556/1‑250 (MQ=255)                                                                             tAATTCCACACTTCCCGATGACACTCAGCATGACGGACAATCGGTGCGTTTTCTCTAGAACAAATCGGTCAATGAAGCATGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACTGATTTATTCAGGTTAAGCCGAAATGCACCGTCTATGACAACCTCGCTTATAACAAAGGCggg                                                                                                                                                                    <  2:382752/250‑1 (MQ=255)                                                                                      acTTCCCGAAGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCTAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCATGTTAAGCCGAAATTCACCGACTATTACACCCTCGCTTATAACAAACGCGGGAAATTACAAc                                                                                                                                                          <  1:200832/251‑1 (MQ=255)                                                                                          cccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa                                                                                                                                                       <  1:104108/250‑1 (MQ=255)                                                                                           ccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAcc                                                                                                                                                     >  1:407933/1‑251 (MQ=255)                                                                                              aTGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGtt                                                                                                                                                  >  2:186729/1‑251 (MQ=255)                                                                                                gACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACACCATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACAGTTAc                                                                                                                                                >  2:150230/1‑251 (MQ=255)                                                                                                   aGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACGCTCCGTTACTCa                                                                                                                                             <  2:211104/251‑1 (MQ=255)                                                                                                                 aCAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgcgc                                                                                                                                >  1:241255/1‑250 (MQ=255)                                                                                                                       ggTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCa                                                                                                                          >  1:364586/1‑250 (MQ=255)                                                                                                                               ttCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACt                                                                                                                                               >  1:344665/1‑221 (MQ=255)                                                                                                                               ttCTCTATAACAAAGCGCTCAATGAAGCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACt                                                                                                                                               <  2:344665/221‑1 (MQ=255)                                                                                                                                   ctATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTTCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAttt                                                                                                               <  1:92333/249‑1 (MQ=255)                                                                                                                                           aaTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCAGAACTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGAGCACATCAACACTGTATTTGAGTGGTAGc                                                                                                     >  1:287811/1‑251 (MQ=255)                                                                                                                                               gCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACATCATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGATTATAACAAACGCGGGACATTACAACTCACCGTTACTCAGCAACTCGGGCGCACCACAACACTGTAtttacgtggacgccatc                                                                                                 >  2:252548/1‑237 (MQ=255)                                                                                                                                               gCTCAATGAATCAGGCACGAATATGCAGTTAGGGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACAGACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                                                                 <  2:189724/251‑1 (MQ=255)                                                                                                                                                           aGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGATATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGCAATTACAACTCACCGTTACTCAGCACCTCGAGCGCCCATCAAAACTGTATTTGAGTGGTAGCCATCAAACTTCTTggg                                                                                     >  2:361817/1‑251 (MQ=255)                                                                                                                                                            ggCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                                                                    >  2:22935/1‑251 (MQ=255)                                                                                                                                                                    tcttCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGCTACAACATACAGTCGAATGAATGGCTACAACATTGACACCCAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTAACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCCTCCACCTTATTGAGGAACGAGt                                                                             >  2:293501/3‑250 (MQ=255)                                                                                                                                                                       tCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGCACGAGTAATg                                                                         >  2:405878/1‑251 (MQ=255)                                                                                                                                                                         aGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                                                                       <  1:22935/251‑1 (MQ=255)                                                                                                                                                                                  gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGCGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                                                               >  1:365023/1‑250 (MQ=255)                                                                                                                                                                                   ttcccGTTATTCGAACAGCGGATATTTTAATGTCGCTGATAAAACATACAGTCGAATGAATGGCTACAACATTTAAACACAGGACGGAGTTATTCAGAGTAAGCCGAACTTCACCGAATATTACAACCACGCTTATAACAAAAGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTACACATCAAAATTCTTGGGGAACACGTAATGTCGATgggtaat                                                             >  2:358775/4‑245 (MQ=255)                                                                                                                                                                                                                   tCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGa                              <  1:301251/250‑1 (MQ=255)                                                                                                                                                                                                                        gATACAACATACAGGCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCaa                        <  1:186729/251‑1 (MQ=255)                                                                                                                                                                                                                                aTACAGTCGAATGAATGGCTACAACATTGAAACCCAGGACGGCGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTCCAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAttt                                                                                                               >  2:198839/1‑156 (MQ=255)                                                                                                                                                                                                                                aTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAttt                                                                                                               <  1:198839/156‑1 (MQ=255)                                                                                                                                                                                                                                aTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGc                >  1:397966/1‑251 (MQ=255)                                                                                                                                                                                                                                              atgGCTACCACATTGAAACACCGGACGGAGTTCTTCCGGTTAAGCCGAAATTCACCGCCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACCCACCGTTACTCAGCAACTCGGGCGCCCATCAACACTGTATTTGAGTGGTAGCCATCAAACTTCTTGGGGAACGAGTAATGTCGGCATGCAATTCCAGGCTGGATTCAATACTGCGTTCGAAGATGTCACCTGGCCGTTCAGCTATAGCCtc  >  2:203756/1‑250 (MQ=255)                                                                                                                                                                                                                                              atgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGaa                                >  1:365218/1‑221 (MQ=255)                                                                                                                                                                                                                                              atgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGaa                                <  2:365218/221‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              GTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTG  >  NC_000913/4546351‑4546838

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