Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,820	A→G	S171S (AGT→AGC)	ecpC ←	putative fimbrial usher protein EcpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,820	0	A	G	100.0%	238.1 / NA	65	S171S (AGT→AGC)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/23);  total (42/23)

TACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGT  >  NC_000913/308689‑308954                                                                                                                                    |                                                                                                                                       gagaGCGGGCCGCCGACCCCCACTTGATGTCCGCGGAGTGAAGTACCGTTATCCAGCGAGAGTTAGCTGGTTTTTTTGCGCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAg                                                                                                                                <  2:825780/137‑1 (MQ=255)  aCAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAgg                                                                                                                               >  1:58717/1‑140 (MQ=255)  aCAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAgg                                                                                                                               >  1:970317/1‑140 (MQ=255)     gCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtg                                                                                                                            >  1:755278/1‑140 (MQ=255)      gggcccgccgACCACCACATGATGTCCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtgt                                                                                                                           <  2:101710/136‑1 (MQ=255)      cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtgt                                                                                                                           >  2:90639/1‑140 (MQ=255)      cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtgt                                                                                                                           >  2:789192/1‑140 (MQ=255)      cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtgt                                                                                                                           >  2:659886/1‑140 (MQ=255)      cGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtgt                                                                                                                           >  2:618393/1‑140 (MQ=255)      cGAGCCGTCGAGCACCACATGATGTTCGCGCAGAGCAGTAACGTTATTCAGCGACAGATAGCTGGATGCGTTGCTCCCGCTGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGgtgt                                                                                                                           >  2:1116859/1‑140 (MQ=255)       gAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGtt                                                                                                                          >  1:677687/1‑140 (MQ=255)        aGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTa                                                                                                                         >  1:1079312/1‑140 (MQ=255)        aGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTa                                                                                                                         >  1:101710/1‑140 (MQ=255)                       caTGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCg                                                                                                          >  1:376372/1‑140 (MQ=255)                       caTGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCg                                                                                                          >  2:20763/1‑140 (MQ=255)                               tCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTc                                                                                                  >  2:1117924/1‑140 (MQ=255)                               tCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTc                                                                                                  >  2:258522/1‑140 (MQ=255)                               tCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTc                                                                                                  >  1:688000/1‑140 (MQ=255)                                 gcgcGGTGCATTATCTTTATCCAGCGAAAGATAGCGGGTTTTGTTCCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGc                                                                                                <  2:635922/140‑1 (MQ=255)                                 gcgcGGTGAAGTATCTTTATTCAGCGAAAGATGGCTGGGTTTTTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGc                                                                                                <  2:970317/140‑1 (MQ=255)                                 gcgcAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGc                                                                                                >  1:1091077/1‑140 (MQ=255)                                 gcgcAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGc                                                                                                >  1:748230/1‑140 (MQ=255)                                  cgcAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCt                                                                                               >  1:635922/1‑140 (MQ=255)                                           gtAACTTTATTCAGCGAAAGATAGTTGGTTGTGTTGCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCg                                                                                      <  2:1079312/140‑1 (MQ=255)                                                 ttATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTAc                                                                                >  1:292242/1‑140 (MQ=255)                                                 ttATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTAc                                                                                >  2:138830/1‑140 (MQ=255)                                                 ttATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTAc                                                                                >  1:1122090/1‑140 (MQ=255)                                                 ttATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTAc                                                                                >  2:1039129/1‑140 (MQ=255)                                                        gCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGccc                                                                         <  1:138830/140‑1 (MQ=255)                                                           aCAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAgc                                                                      >  1:695637/1‑140 (MQ=255)                                                                     ggttgttTTCCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCt                                                            <  2:375700/137‑1 (MQ=255)                                                                         ttgtTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGAc                                                        <  2:695943/139‑1 (MQ=255)                                                                           gtTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCa                                                      >  2:797338/1‑140 (MQ=255)                                                                             tGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACt                                                    >  2:496590/1‑140 (MQ=255)                                                                              gCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTa                                                   >  2:235714/1‑140 (MQ=255)                                                                                tCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGc                                                 <  1:1117924/140‑1 (MQ=255)                                                                                     ccgTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTgcag                                            >  1:829661/1‑140 (MQ=255)                                                                                                cTGGTTGTTATAGACGCCCAAGTTATAGATCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCACAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGgc                                 <  2:688000/140‑1 (MQ=255)                                                                                                 tGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGgcg                                <  1:90639/140‑1 (MQ=255)                                                                                                  ggTTGTTATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGgcgc                               >  2:113643/1‑140 (MQ=255)                                                                                                   gttgttATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCa                              >  1:1143755/1‑140 (MQ=255)                                                                                                   gttgttATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCa                              >  2:965321/1‑140 (MQ=255)                                                                                                   gttgttATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCa                              <  1:659886/140‑1 (MQ=255)                                                                                                      gttATAGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGc                           <  2:376372/140‑1 (MQ=255)                                                                                                           aGACGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTgaga                      <  2:292242/140‑1 (MQ=255)                                                                                                             aCGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATc                    >  2:155480/1‑140 (MQ=255)                                                                                                             aCGCCCAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATc                    >  1:980290/1‑140 (MQ=255)                                                                                                                 ccAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTg                         >  2:30804/1‑131 (MQ=255)                                                                                                                 ccAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTg                         <  1:30804/131‑1 (MQ=255)                                                                                                                 ccAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGc                >  1:631084/1‑140 (MQ=255)                                                                                                                 ccAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGc                >  1:820777/1‑140 (MQ=255)                                                                                                                  cAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGc                           >  2:384377/1‑128 (MQ=255)                                                                                                                  cAAGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGc                           <  1:384377/128‑1 (MQ=255)                                                                                                                    aGTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTgc             <  1:797338/140‑1 (MQ=255)                                                                                                                     gTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTgcg            <  2:1122090/140‑1 (MQ=255)                                                                                                                     gTTATAGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTgcg            <  1:832545/140‑1 (MQ=255)                                                                                                                         taGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGtt        >  2:313033/1‑140 (MQ=255)                                                                                                                          aGCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTg       >  2:907975/1‑140 (MQ=255)                                                                                                                           gCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTgt      <  1:501192/140‑1 (MQ=255)                                                                                                                           gCTCAGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTgt      <  1:496590/140‑1 (MQ=255)                                                                                                                               aGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTgtcagt  >  1:123944/1‑140 (MQ=255)                                                                                                                               aGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTgtcagt  <  1:1139066/140‑1 (MQ=255)                                                                                                                               aGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTgtcagt  >  1:599624/1‑140 (MQ=255)                                                                                                                               aGATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTgtcagt  >  1:142263/1‑140 (MQ=255)                                                                                                                                gATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTgcg            <  2:7562/129‑1 (MQ=255)                                                                                                                                gATTGCTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTgcg            >  1:7562/1‑129 (MQ=255)                                                                                                                                    |                                                                                                                                       TACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGT  >  NC_000913/308689‑308954

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