Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	305,990	G→A	I276I (ATC→ATT)	ecpD ←	fimbrial adhesin EcpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	305,990	0	G	A	97.3%	113.6 / ‑5.8	37	I276I (ATC→ATT)	ecpD	fimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (15/21);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGAT  >  NC_000913/305764‑306222                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         aTAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCTCCATCTGCGTAAAGAAGTTGCTCGAGAAGAGAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGGACGCCACGACGGCCGAGGCTAACGACGAGTTGGCGATCGGCGGGAATATATGAATTGAGGTGTTTCTTGGTCCgttt                                                                                                                                                                                                                  <  2:18530/251‑1 (MQ=255) aTAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCgttt                                                                                                                                                                                                                  >  1:254440/1‑251 (MQ=255)   aGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCTCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCgtttgt                                                                                                                                                                                                                <  1:399506/251‑1 (MQ=255)         ttGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCgtttgtt                                                                                                                                                                                                               <  2:152446/246‑1 (MQ=255)              gtTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGtat                                                                                                                                                                                                     >  1:78346/1‑251 (MQ=255)                           aTCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCtt                                                                                                                                                                                        <  2:180663/251‑1 (MQ=255)                            tCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCtt                                                                                                                                                                                        >  1:258753/1‑250 (MQ=255)                             cccGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAAGTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTAc                                                                                                                                                                                      <  1:126860/251‑1 (MQ=255)                                 aGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACt                                                                                                                                                                                  <  2:220865/251‑1 (MQ=255)                                 aGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACt                                                                                                                                                                                  <  2:79073/251‑1 (MQ=255)                                   gTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTa                                                                                                                                                                                <  2:300432/251‑1 (MQ=255)                                        ccATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCt                                                                                                                                                                           <  1:83671/251‑1 (MQ=255)                                         cATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTg                                                                                                                                                                          <  2:91287/251‑1 (MQ=255)                                              gCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCGCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCTGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGTTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGa                                                                                                                                                                         <  1:177732/247‑1 (MQ=255)                                                 tAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCt                                                                                                                                                                   <  1:400974/250‑1 (MQ=255)                                                             tCGAGAAGAACACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGc                                                                                                                                                      >  2:393646/1‑251 (MQ=255)                                                                agaagaAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATcc                                                                                                                                                                            <  1:361054/226‑1 (MQ=255)                                                                agaagaAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATcc                                                                                                                                                                            >  2:361054/1‑226 (MQ=255)                                                                                     cTGCGCTCTTCATCTCGTTGAAGGAGAAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGTCCCGACGGCCGAGGCTAACGACGAGTTGGAGATCGCCGGGCATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCg                                                                                                                              >  2:256044/1‑251 (MQ=255)                                                                                           tctTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGccc                                                                                                                           >  2:165976/1‑248 (MQ=255)                                                                                              tCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTc                                                                                                                     >  1:242582/1‑251 (MQ=255)                                                                                                                   aaTAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGAt                                                                                                 <  1:13101/250‑1 (MQ=255)                                                                                                                       ggCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGTGTCCGTTTGTTTGCAGGGtat                                                                                                                                                                                                     >  2:285627/1‑146 (MQ=255)                                                                                                                       ggCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGtat                                                                                                                                                                                                     <  1:285627/146‑1 (MQ=255)                                                                                                                          ggTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACtt                                                                                         >  2:372927/1‑251 (MQ=255)                                                                                                                                     ccGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCAGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCg                                                                                >  2:61837/1‑249 (MQ=255)                                                                                                                                           tCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGt                                                                         <  1:222217/250‑1 (MQ=255)                                                                                                                                                  aCTGCCAGTCAGACTGAACTGCAGATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGATTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGGCGGCGTTCCCGTCACCCAGAGTCGGGACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGATTGGGTATTTATCa                                                                 <  2:237709/251‑1 (MQ=255)                                                                                                                                                       cATTCAGACTTAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGTAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGt                                                            <  2:10196/251‑1 (MQ=255)                                                                                                                                                             gACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGAttt                                                      >  1:381543/1‑251 (MQ=255)                                                                                                                                                                                cGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTg                                    >  2:125438/1‑250 (MQ=255)                                                                                                                                                                                  ccccGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGAct                                 <  2:39315/251‑2 (MQ=255)                                                                                                                                                                                    ccGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGGCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATGTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGtt                               <  2:201973/251‑1 (MQ=255)                                                                                                                                                                                          gCCGAGGCTAACGACGAGTTGGCGATCGCCGGGCATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAAc                         >  1:248529/1‑251 (MQ=255)                                                                                                                                                                                          gCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAAc                         <  1:10216/251‑1 (MQ=255)                                                                                                                                                                                             gAGGCTAACGACGAGTTGGCGATCGCCGGGAATCTCTGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACAt                       >  2:343265/1‑250 (MQ=255)                                                                                                                                                                                                                 gATCGCCGGGAATATATGAATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGAt  >  2:96254/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         ATAGATCTTTGGTGTTGATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGAT  >  NC_000913/305764‑306222

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