Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,132,041	G→A	5.5%	A53V (GCA→GTA)	fimD ←	type I pili subunit FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,132,041	0	G	A	5.5%	224.3 / 2.1	73	A53V (GCA→GTA)	fimD	type I pili subunit FimD
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base A (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ATGGCTGCAGCTCACCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGTGGG  >  NC_002947/2131895‑2132187                                                                                                                                                     |                                                                                                                                                   aTGGCTGCAGCTCACCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCgctgct                                                                                                                                                  >  5:199855/1‑149 (MQ=255) aTGGCTGCAGCTCACCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCgctgct                                                                                                                                                  >  1:139220/1‑149 (MQ=255)    gCTGCAGCTCACCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTg                                                                                                                                                <  2:135647/148‑1 (MQ=255)        cAGCTCACCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCgg                                                                                                                                             <  3:61229/147‑1 (MQ=255)              aCCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAg                                                                                                                                     >  1:107626/1‑149 (MQ=255)              aCCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAg                                                                                                                                     >  1:70658/1‑149 (MQ=255)                cgggcGGGCCTGGGTGGAAATCCAGCTCATTCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCCCTCTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAg                                                                                                                                     <  1:101546/145‑1 (MQ=255)                  tCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCAc                                                                                                                                 >  6:311500/1‑149 (MQ=255)                    cGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACtt                                                                                                                               >  5:159400/1‑149 (MQ=255)                     gggCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTg                                                                                                                              <  7:95351/149‑1 (MQ=255)                              tggAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAgg                                                                                                                       >  1:174729/1‑147 (MQ=255)                              tggAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAgg                                                                                                                       >  4:285442/1‑147 (MQ=255)                                gAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGGGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGaa                                                                                                                     <  6:83855/147‑1 (MQ=255)                                gAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGaa                                                                                                                     <  1:366353/147‑1 (MQ=255)                                 aaaTCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGcc                                                                                                                  <  6:199855/149‑1 (MQ=255)                                  aaTCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCg                                                                                                                 >  5:41917/1‑149 (MQ=255)                                      cAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGctc                                                                                                              <  4:348127/148‑1 (MQ=255)                                       aGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTg                                                                                                            <  3:285442/149‑1 (MQ=255)                                       aGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTg                                                                                                            <  4:113336/149‑1 (MQ=255)                                          tCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTACTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGaa                                                                                                         >  1:40734/1‑149 (MQ=255)                                          tCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTACTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGaa                                                                                                         >  1:46363/1‑149 (MQ=255)                                            aTGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAAcc                                                                                                       >  5:21538/1‑149 (MQ=255)                                                tGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCg                                                                                                     >  8:207945/1‑147 (MQ=255)                                                    cGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCt                                                                                               <  6:159400/149‑1 (MQ=255)                                                    cGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCt                                                                                               <  3:282191/149‑1 (MQ=255)                                                      gCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGc                                                                                             >  7:170456/1‑149 (MQ=255)                                                      gCAGGTTCAAGATTGACCTGTACCTCTACACGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTAGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGc                                                                                             >  6:293261/1‑149 (MQ=255)                                                          gTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTgg                                                                                         <  3:209380/149‑1 (MQ=255)                                                           ttCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTgg                                                                                         >  1:121746/1‑148 (MQ=255)                                                             cAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAg                                                                                      >  2:357978/1‑149 (MQ=255)                                                              aaGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAg                                                                                      <  4:173847/148‑1 (MQ=255)                                                               aGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTa                                                                                    <  1:329855/149‑1 (MQ=255)                                                                     aCCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGgcgc                                                                              >  4:192470/1‑149 (MQ=255)                                                                       cTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCt                                                                             >  6:125830/1‑148 (MQ=255)                                                                       cTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCt                                                                             >  8:332879/1‑148 (MQ=255)                                                                          tACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTgaga                                                                         <  3:192470/149‑1 (MQ=255)                                                                          tACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTgaga                                                                         <  2:70658/149‑1 (MQ=255)                                                                          tACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTgaga                                                                         <  8:170456/149‑1 (MQ=255)                                                                             ctctACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGc                                                                       >  2:253491/1‑148 (MQ=255)                                                                                tACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAg                                                                                                                        >  4:128013/1‑96 (MQ=255)                                                                                tACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAg                                                                                                                        <  3:128013/96‑1 (MQ=255)                                                                                  cccGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAgg                                                                   <  6:41917/147‑1 (MQ=255)                                                                                   ccGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTgg                                                                >  2:173546/1‑149 (MQ=255)                                                                                   ccGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTgg                                                                <  2:121746/149‑1 (MQ=255)                                                                                    cGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGa                                                               <  6:21538/149‑1 (MQ=255)                                                                                    cGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGa                                                               <  2:107626/149‑1 (MQ=255)                                                                                    cGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGa                                                               <  5:311500/149‑1 (MQ=255)                                                                                       tAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCt                                                            <  5:64469/149‑1 (MQ=255)                                                                                         gCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTa                                                           <  2:139220/148‑1 (MQ=255)                                                                                           ggCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTAcc                                                         >  5:185433/1‑148 (MQ=255)                                                                                            gCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTg                                                                            >  2:29361/1‑128 (MQ=255)                                                                                            gCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTg                                                                            <  1:29361/128‑1 (MQ=255)                                                                                            gCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCAGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCca                                                       <  7:280859/149‑1 (MQ=255)                                                                                              cAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCcagc                                                     <  7:207945/149‑1 (MQ=255)                                                                                                 gcgcGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAgg                                                  >  1:125525/1‑149 (MQ=255)                                                                                                   gcgAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGca                                                >  5:154979/1‑149 (MQ=255)                                                                                                                     tGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAAtgcgctgcgc                                >  8:251423/1‑149 (MQ=255)                                                                                                                       gTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGcgctgcgctg                              >  8:201574/1‑149 (MQ=255)                                                                                                                           aCATGCTCAGCGACCAGGCACTCGCTCCTTGCGGTAGGCCAGGCACTTGCCGCAGGAAGACCCTCTGGAACCCGGATTCTGCCTGGGAGTAGGGCGCTCAGAGCAAGGTGGAGATACCCAGCAGGCACACCAATGCGCTGCGCTggagg                          >  3:33934/1‑149 (MQ=255)                                                                                                                            cATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTggcggc                         <  7:332879/149‑1 (MQ=255)                                                                                                                            cATGCTCAGCGACCAGGCACTCGCTACTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAACAATGCGCTGCGCTggcggc                         <  2:40734/149‑1 (MQ=255)                                                                                                                            cATGCTCAGCGACCAGGCACTCGCTACTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAACAATGCGCTGCGCTggcggc                         <  2:46363/149‑1 (MQ=255)                                                                                                                                cTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCa                      <  1:173546/148‑1 (MQ=255)                                                                                                                                cTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTATGCCTGGGAGTAGGGCGCAGAGAGCAAGCTGGAGCTACCCAGCAGGCACACCAATGCGCTGCGCTGGCGGCGCa                      >  6:131182/1‑148 (MQ=255)                                                                                                                                 tCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCa                      >  4:125976/1‑147 (MQ=255)                                                                                                                                 tCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCa                      >  8:296110/1‑147 (MQ=255)                                                                                                                                    gCGATCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGa                                                               <  2:289667/103‑1 (MQ=255)                                                                                                                                    gCGATCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGa                                                               >  1:289667/1‑103 (MQ=255)                                                                                                                                         cAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGcac                                               >  2:23962/1‑114 (MQ=255)                                                                                                                                         cAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGcac                                               <  1:23962/114‑1 (MQ=255)                                                                                                                                         cAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTgcgc              >  3:3741/1‑147 (MQ=255)                                                                                                                                         cAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTgcgc              >  5:515/1‑147 (MQ=255)                                                                                                                                               cTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAg      <  6:515/149‑1 (MQ=255)                                                                                                                                                tCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGt     >  3:33877/1‑149 (MQ=255)                                                                                                                                                tCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGT‑GAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGTg    >  4:176800/1‑149 (MQ=255)                                                                                                                                                  gctgctTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGt     <  7:251423/147‑1 (MQ=255)                                                                                                                                                  gctgctTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGt     <  2:174729/147‑1 (MQ=255)                                                                                                                                                  gctgctTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGTgg   >  8:235890/1‑149 (MQ=255)                                                                                                                                                   ctgctTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGTggg  <  2:63593/149‑1 (MQ=255)                                                                                                                                                   ctgctTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGTggg  <  1:253491/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ATGGCTGCAGCTCACCGTCCGGGCCTGGCTGGAAATCCAGCTCATGCTGGCCGGCAGGTTCAAGATTGACCTGTACCTCTACCCGGTAGCGGCCAGGCGCGAGCCC‑‑CTGATCGCTGGTCAACATGCTCAGCGACCAGGCACTCGCTGCTTGCGGTTGGCCAGGCACTTGGCGCAGGAAGCCGCTCTGGAACCCGGAGTCTGCCTGGGAGTAGGGCGCTGAGAGCAAGGTGGAGCTACCCAGCAGGCACAGCAATGCGCTGCGCTGGCGGCGCATAGTGCGCCAGCGTAGTGGG  >  NC_002947/2131895‑2132187

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