Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,547,373	C→A	L760M (CTG→ATG)	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,547,373	0	C	A	95.5%	302.9 / 7.1	89	L760M (CTG→ATG)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (33/52);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

AACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCAC  >  NC_000913/4547238‑4547507                                                                                                                                         |                                                                                                                                       aaCGATACGGTGGTGCTTGTTAAGGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  2:923037/140‑1 (MQ=255) aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  1:958428/140‑1 (MQ=255) aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  2:118062/140‑1 (MQ=255) aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  2:418833/140‑1 (MQ=255) aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  1:23144/140‑1 (MQ=255) aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  2:606695/140‑1 (MQ=255) aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGa                                                                                                                                    <  1:1203439/140‑1 (MQ=255)     aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATAcc                                                                                                                                >  1:920508/1‑140 (MQ=255)       aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGTATACCaa                                                                                                                              >  1:60822/1‑140 (MQ=255)        cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAAt                                                                                                                             >  2:242007/1‑140 (MQ=255)           tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCTCT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATGcc                                                                                                                          <  1:1113958/140‑1 (MQ=255)           tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATAcc                                                                                                                          >  2:654794/1‑140 (MQ=255)           tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATAcc                                                                                                                          <  1:467298/140‑1 (MQ=255)           tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATAcc                                                                                                                          <  1:31761/140‑1 (MQ=255)             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCt                                                                                                                        >  1:676466/1‑140 (MQ=255)             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCt                                                                                                                        >  2:611876/1‑140 (MQ=255)                   gTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGa                                                                                                                  >  1:448960/1‑140 (MQ=255)                    ttAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGAt                                                                                                                 <  2:766169/140‑1 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                                                                                                                     >  1:611134/1‑132 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                                                                                                                     <  2:611134/132‑1 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACg                                                                                                             <  2:570088/140‑1 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACg                                                                                                             <  2:1277883/140‑1 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACg                                                                                                             <  1:184792/140‑1 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACg                                                                                                             >  2:453987/1‑140 (MQ=255)                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACg                                                                                                             <  1:1074576/140‑1 (MQ=255)                            ccTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCg                                                                                                          >  2:1106616/1‑140 (MQ=255)                             cTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGGACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGAt                                                                                                        >  1:345091/1‑140 (MQ=255)                                 cgcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAg                                                                                                    >  2:182616/1‑140 (MQ=255)                                 cgcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAg                                                                                                    >  2:1236067/1‑140 (MQ=255)                                         aTGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAcgcg                                                                                            <  2:1005195/140‑1 (MQ=255)                                               aaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCt                                                                                      <  2:330653/140‑1 (MQ=255)                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATTCCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACg                                                                                  <  2:533100/140‑1 (MQ=255)                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACg                                                                                  <  1:301198/140‑1 (MQ=255)                                                          cAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATTTCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTcc                                                                           <  1:1021934/140‑1 (MQ=255)                                                          cAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTcc                                                                           <  1:999729/140‑1 (MQ=255)                                                          cAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTcc                                                                           <  2:676466/140‑1 (MQ=255)                                                          cAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTcc                                                                           <  1:527611/140‑1 (MQ=255)                                                          cAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTcc                                                                           <  1:222170/140‑1 (MQ=255)                                                          cAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTcc                                                                           <  1:1269294/140‑1 (MQ=255)                                                                  ggTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTg                                                                   <  1:1007848/140‑1 (MQ=255)                                                                       gTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCg                                                              <  1:888636/140‑1 (MQ=255)                                                                         aCCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGAt                                                            >  2:620614/1‑140 (MQ=255)                                                                         aCCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGAt                                                            <  2:1279992/140‑1 (MQ=255)                                                                           cGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCg                                                          >  2:1167173/1‑140 (MQ=255)                                                                           cGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCg                                                          >  2:201413/1‑140 (MQ=255)                                                                           cGACTGGCGTGGTTATGCCGCGCTGCCCTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGACTGATAACGCCGATTTAGATAACGCTGATGCTAAAGATGATCCCACTCGGGGGGCGATCg                                                          >  2:216024/1‑140 (MQ=255)                                                                            gACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGt                                                         <  2:469776/140‑1 (MQ=255)                                                                             aCTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTg                                                        >  2:1000769/1‑140 (MQ=255)                                                                             aCTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTg                                                        >  1:205747/1‑140 (MQ=255)                                                                             aCTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTg                                                        >  2:1113363/1‑140 (MQ=255)                                                                                ggCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGTATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGa                                                     <  2:60822/140‑1 (MQ=255)                                                                                  cGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGc                                                   <  2:490708/140‑1 (MQ=255)                                                                                    tGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCag                                                 <  1:344451/140‑1 (MQ=255)                                                                                               tgctgcCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcg                                      >  2:555518/1‑140 (MQ=255)                                                                                               tgctgcCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcg                                      >  2:404721/1‑140 (MQ=255)                                                                                                gctgcCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgc                                     >  2:967157/1‑140 (MQ=255)                                                                                                gctgcCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgc                                     >  1:823142/1‑140 (MQ=255)                                                                                                 ctgcCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgcg                                    >  1:527952/1‑140 (MQ=255)                                                                                                 ctgcCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGG‑TACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgcgc                                   >  2:1250346/1‑140 (MQ=255)                                                                                                    ccTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGt                                 <  1:654794/140‑1 (MQ=255)                                                                                                    ccTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGt                                 <  1:999224/140‑1 (MQ=255)                                                                                                    ccTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGt                                 <  2:345091/140‑1 (MQ=255)                                                                                                      ttATGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTg                               >  1:1169368/1‑140 (MQ=255)                                                                                                         tGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGa                            <  2:396874/140‑1 (MQ=255)                                                                                                         tGCCACT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGa                            <  1:242007/140‑1 (MQ=255)                                                                                                           ccACT‑GAATATCGGGAAAATAGAGTGGCGATGG‑TACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATaa                         <  1:1250346/140‑1 (MQ=255)                                                                                                              cT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATaaaa                       <  1:354606/140‑1 (MQ=255)                                                                                                              cT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATaaaa                       >  1:855636/1‑140 (MQ=255)                                                                                                              cT‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATaaaa                       >  1:562805/1‑140 (MQ=255)                                                                                                               t‑GAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAAc                      <  1:1134292/140‑1 (MQ=255)                                                                                                                 gAATATCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACt                     >  2:1082860/1‑140 (MQ=255)                                                                                                                    tatCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCt                  >  2:1294903/1‑140 (MQ=255)                                                                                                                    tatCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCt                  >  1:785978/1‑140 (MQ=255)                                                                                                                    tatCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCt                  >  2:563500/1‑140 (MQ=255)                                                                                                                      tCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCa                <  2:716936/140‑1 (MQ=255)                                                                                                                      tCGGGAAAATAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCa                >  1:97872/1‑140 (MQ=255)                                                                                                                            aaaTAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGc          >  1:36625/1‑140 (MQ=255)                                                                                                                             aaTAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCt         <  1:182616/140‑1 (MQ=255)                                                                                                                              aTAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTg        <  2:221210/140‑1 (MQ=255)                                                                                                                               tAGAGTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGa       <  2:413822/140‑1 (MQ=255)                                                                                                                                 gagTGGCGATGGATACCAATACCCTGGCTGATATCGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc     <  2:1314677/140‑1 (MQ=255)                                                                                                                                 gagTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc     <  1:1236067/140‑1 (MQ=255)                                                                                                                                 gagTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc     <  1:788200/140‑1 (MQ=255)                                                                                                                                 gagTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc     <  2:1093280/140‑1 (MQ=255)                                                                                                                                 gagTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGAcc     <  1:1000769/140‑1 (MQ=255)                                                                                                                                   gTGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCca   <  2:920508/140‑1 (MQ=255)                                                                                                                                    tGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcac  <  1:1106616/140‑1 (MQ=255)                                                                                                                                    tGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcac  >  2:587376/1‑140 (MQ=255)                                                                                                                                    tGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcac  <  2:176649/140‑1 (MQ=255)                                                                                                                                    tGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCca   <  1:973456/139‑1 (MQ=255)                                                                                                                                    tGGCGATGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCca   >  2:973456/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       AACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACT‑GAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCAC  >  NC_000913/4547238‑4547507

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑