Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,547,319	C→A	R742S (CGT→AGT)	fimD →	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,547,319	0	C	A	90.9%	312.2 / 17.9	99	R742S (CGT→AGT)	fimD	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (39/51);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACG  >  NC_000913/4547224‑4547415                                                                                                |                                                                                                 tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTgg                                                                                              <  2:1375990/100‑1 (MQ=255)  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                             <  1:727707/100‑1 (MQ=255)  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                             <  1:247702/100‑1 (MQ=255)  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                             <  1:1518695/100‑1 (MQ=255)    ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGtt                                                                                            >  1:1879359/1‑99 (MQ=255)    ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                           >  2:274294/1‑100 (MQ=255)    ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                           >  1:370389/1‑100 (MQ=255)         ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCg                                                                                      <  2:615225/100‑1 (MQ=255)         ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCg                                                                                      <  1:1371611/100‑1 (MQ=255)          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGt                                                                                     >  2:1942939/1‑100 (MQ=255)          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                     <  2:914617/100‑1 (MQ=255)          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                     <  2:867825/100‑1 (MQ=255)          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                     <  2:214098/100‑1 (MQ=255)           gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                    >  1:144119/1‑100 (MQ=255)               aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgc                                                                                >  2:382397/1‑100 (MQ=255)               aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgc                                                                                >  1:1030816/1‑100 (MQ=255)               aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgc                                                                                <  1:936148/100‑1 (MQ=255)                  gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCta                                                                             <  1:1011297/100‑2 (MQ=255)                   aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                            <  1:1808118/100‑1 (MQ=255)                   aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGCTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                            <  1:1042108/100‑1 (MQ=255)                     aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                          <  1:1441479/100‑1 (MQ=255)                      cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                         <  2:1113124/100‑1 (MQ=255)                      cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                         <  2:1776301/100‑1 (MQ=255)                      cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                         <  1:1939229/100‑1 (MQ=255)                      cGGTGGTGCCTGTTAAAGCACCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGAGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                         >  1:1546394/1‑100 (MQ=255)                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                        <  2:279890/100‑1 (MQ=255)                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                        <  2:370391/100‑1 (MQ=255)                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                        <  1:1243102/100‑1 (MQ=255)                         tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCAc                                                                      <  1:797210/100‑1 (MQ=255)                          ggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                     >  1:844116/1‑100 (MQ=255)                          ggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                     <  1:1047720/100‑1 (MQ=255)                          ggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                     <  2:514677/100‑1 (MQ=255)                           gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                    <  1:1888180/100‑1 (MQ=255)                           gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                    <  1:967336/100‑1 (MQ=255)                            tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                                   >  2:347251/1‑100 (MQ=255)                            tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCCAAAACCAGACGGGGGTGCGCACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                                   <  1:1060715/100‑1 (MQ=255)                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  <  1:18608/100‑1 (MQ=255)                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  >  2:1770003/1‑100 (MQ=255)                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  <  2:1491195/100‑1 (MQ=255)                              cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                                 >  2:319476/1‑100 (MQ=255)                              cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                                 <  2:1591504/100‑1 (MQ=255)                              cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                                 >  1:345933/1‑100 (MQ=255)                                tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAata                                                                >  2:738425/1‑99 (MQ=255)                                 gTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATc                                                              >  2:765790/1‑100 (MQ=255)                                 gTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATc                                                              >  2:81/1‑100 (MQ=255)                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCCGACGGGGGAGCGTACCGACAGGACTAGTTATGCATTGATGCATTAGGacacagca                                                                  >  1:1261679/1‑90 (MQ=255)                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                            <  2:1944951/100‑1 (MQ=255)                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                            <  2:43994/100‑1 (MQ=255)                                    aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                           >  2:859498/1‑100 (MQ=255)                                     naGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          >  1:1727774/2‑100 (MQ=255)                                     naGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          >  2:1674422/2‑100 (MQ=255)                                     aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          >  1:1728734/1‑100 (MQ=255)                                     aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          >  2:592916/1‑100 (MQ=255)                                     aaGCGCCGGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          <  1:1074290/100‑1 (MQ=255)                                        cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaaa                                                       <  1:1942926/100‑1 (MQ=255)                                          ccTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATa                                                     <  1:588922/100‑1 (MQ=255)                                           cTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATag                                                    >  1:801633/1‑100 (MQ=255)                                             ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                  <  1:1822371/100‑1 (MQ=255)                                              gcgcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGt                                                 <  2:1282619/100‑1 (MQ=255)                                                 cAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                              >  1:429320/1‑100 (MQ=255)                                                 cAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAAGAGAGTGgc                                              >  1:970912/1‑100 (MQ=255)                                                  aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcg                                             >  1:1650959/1‑100 (MQ=255)                                                   aaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcgc                                            >  2:1643545/1‑100 (MQ=255)                                                   aaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcgc                                            <  2:1904178/100‑1 (MQ=255)                                                     aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                          >  2:584912/1‑100 (MQ=255)                                                     aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                          <  2:342820/100‑1 (MQ=255)                                                       aTGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGa                                        <  1:624719/100‑1 (MQ=255)                                                        tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                       <  2:1685267/100‑1 (MQ=255)                                                         gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                      >  1:467237/1‑100 (MQ=255)                                                          cAACAGTCGAAAACCAGACTGGGGTGCGTACCGACTGGAGCTGTTATGACGTGCTGCCTTCTGCCGCTGAATACCGGGTAAATATAGTTGCTCTGTATAc                                     >  1:420518/1‑100 (MQ=255)                                                          cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGCGTGGCGGTGGATAc                                     >  1:1932860/1‑100 (MQ=255)                                                          cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                     >  2:89404/1‑100 (MQ=255)                                                          cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                     >  1:565109/1‑100 (MQ=255)                                                           aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                    <  2:293409/100‑1 (MQ=255)                                                               gTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATa                                <  1:15582/100‑1 (MQ=255)                                                                 cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                              >  2:631952/1‑100 (MQ=255)                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                             >  1:341067/1‑100 (MQ=255)                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                             >  2:1377805/1‑100 (MQ=255)                                                                    aaaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTg                           <  2:711370/100‑1 (MQ=255)                                                                    aaaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTg                           >  1:1767037/1‑100 (MQ=255)                                                                    aaaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTAGGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTg                           <  1:955434/100‑1 (MQ=255)                                                                     aaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                          <  2:1940795/100‑1 (MQ=255)                                                                       ccAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                        <  2:844117/100‑1 (MQ=255)                                                                       ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                        <  2:1363848/100‑1 (MQ=255)                                                                        cAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTg                       >  2:935179/1‑100 (MQ=255)                                                                         aGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGa                      <  2:1852650/100‑1 (MQ=255)                                                                           aCGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATa                    <  2:1879371/100‑1 (MQ=255)                                                                              ggggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACg                 >  1:946561/1‑100 (MQ=255)                                                                                ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc               <  1:1004946/100‑1 (MQ=255)                                                                                 gTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCg              <  1:1263358/100‑1 (MQ=255)                                                                                  tGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGa             >  1:149068/1‑100 (MQ=255)                                                                                    cGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAtt           <  2:345935/100‑1 (MQ=255)                                                                                    cGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCCGGATACCAATACCCTGGCTGATAACGTCGAtt           <  1:1667708/100‑1 (MQ=255)                                                                                    cGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAACAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAtt           >  2:55646/1‑100 (MQ=255)                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt          >  2:650489/1‑100 (MQ=255)                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt          <  1:1643535/100‑1 (MQ=255)                                                                                      tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa         <  1:564584/100‑1 (MQ=255)                                                                                       aCCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAg        <  1:382395/100‑1 (MQ=255)                                                                                             tGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGagatcg  >  2:512256/1‑95 (MQ=255)                                                                                                |                                                                                                 TGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACG  >  NC_000913/4547224‑4547415

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