Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,261,617	C→T	10.1%	A735A (GCC→GCT) ‡	fimV →	motility protein FimV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,261,617	0	C	T	10.1%	160.4 / 6.3	59	A735A (GCC→GCT) ‡	fimV	motility protein FimV
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (4/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CCGACGAACTGATAATCGCTGAGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGTTTGATCCGG  >  NC_002947/2261478‑2261765                                                                                                                                            |                                                                                                                                                      ccGACGAACTGATAATCGCTGAGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAg                                                                                                                                              <  1:263215/149‑1 (MQ=255)  ctactaaCTGATAATCGCTGAGCCTGTGGTGGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGc                                                                                                                                             <  3:69321/144‑1 (MQ=255)        aCTGATAATCGCTGAGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCg                                                                                                                                       <  5:134476/149‑1 (MQ=255)          tGATAATCGCTGAGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTg                                                                                                                                     >  3:289306/1‑149 (MQ=255)            aTAATCGCTGAGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGa                                                                                                                                   <  2:68812/149‑1 (MQ=255)                      aGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTg                                                                                                                           <  4:74340/147‑1 (MQ=255)                         cTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTc                                                                                                                      <  1:255218/149‑1 (MQ=255)                         cTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTc                                                                                                                      <  4:256204/149‑1 (MQ=255)                           gtgCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGa                                                                                                                    >  5:73622/1‑149 (MQ=255)                                 gATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGa                                                                                                              >  3:6963/1‑149 (MQ=255)                                     aGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGt                                                                                                           <  4:279646/148‑1 (MQ=255)                                      gCCCGTACATGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCATCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCa                                                                                                         <  5:246947/149‑1 (MQ=255)                                      gCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGATCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCa                                                                                                         <  5:251392/149‑1 (MQ=255)                                         cGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCAGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGa                                                                                                                    <  7:248657/135‑1 (MQ=255)                                         cGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCAGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGa                                                                                                                    >  8:248657/1‑135 (MQ=255)                                          gTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGc                                                                                                       >  3:293073/1‑147 (MQ=255)                                            aCAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGc                                                                                                   >  5:247550/1‑149 (MQ=255)                                              aGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGCCCTGTCTGATTTCGACCTCGACGTCAGCGAGCCg                                                                                                 >  5:196611/1‑149 (MQ=255)                                                gaagaaGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGt                                                                                               >  2:8700/1‑149 (MQ=255)                                                gaagaaGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGt                                                                                               >  8:155921/1‑149 (MQ=255)                                                gaagaaGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGt                                                                                               >  3:315765/1‑149 (MQ=255)                                                   gaagCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGc                                                                                            <  6:73622/149‑1 (MQ=255)                                                          aGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑gg                                                                                      <  5:55066/147‑1 (MQ=255)                                                            ccGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGc                                                                                  <  7:155921/149‑1 (MQ=255)                                                               gaagaaCAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGc                                                                                  <  2:83505/146‑1 (MQ=255)                                                               gaagaaCAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGc                                                                                  >  1:83505/1‑146 (MQ=255)                                                                aagaaCAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCt                                                                              <  3:254681/149‑1 (MQ=255)                                                                  gaaCAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGa                                                                            <  4:293073/149‑1 (MQ=255)                                                                   aaCAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAg                                                                           <  1:8700/149‑1 (MQ=255)                                                                     cAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGAc                                                                                                             >  7:220796/1‑114 (MQ=255)                                                                     cAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGAc                                                                                                             <  8:220796/114‑1 (MQ=255)                                                                             tGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGa                                                                   <  4:315765/147‑1 (MQ=255)                                                                              cgtatctaTGCCGAGCTGAAAGCCGACGCCGACGCCGATTTCGATGCCATGCTGGCTCAGGACGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGc                                                                 <  8:136628/141‑1 (MQ=255)                                                                               ccGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAg                                                                                                    >  3:73580/1‑113 (MQ=255)                                                                               ccGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAg                                                                                                    <  4:73580/113‑1 (MQ=255)                                                                                cGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGc                                                                 >  7:269061/1‑146 (MQ=255)                                                                                 gTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGa                                                             <  8:269061/149‑1 (MQ=255)                                                                                      tGCCGAGCTGCAAGCCGACGCCGAAGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGCCCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCATGATGCGCTGGAGc                                                                 >  7:136628/1‑140 (MQ=255)                                                                                                           cgacgccgaTTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTg                                                                     >  4:207151/1‑115 (MQ=255)                                                                                                           cgacgccgaTTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTg                                                                     <  3:207151/115‑1 (MQ=255)                                                                                                             acgccgaTTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCATGCGGTGGACCTGTCTGATTTCGACATCGACGTCAGCGAGCCGGTTGGAAC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTGACCGCCGCACTGGCt                                 >  6:276465/1‑149 (MQ=255)                                                                                                                 cgaTTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCt                             >  2:42800/1‑149 (MQ=255)                                                                                                                 cgaTTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCt                             >  8:194864/1‑149 (MQ=255)                                                                                                                  gaTTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCtt                            <  6:196611/149‑1 (MQ=255)                                                                                                                       cGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGa                         >  2:77519/1‑147 (MQ=255)                                                                                                                       cGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACa                       >  4:68501/1‑149 (MQ=255)                                                                                                                         aGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGc                     >  7:92055/1‑149 (MQ=255)                                                                                                                             cATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAAc                 >  8:34588/1‑149 (MQ=255)                                                                                                                               tGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCg               <  7:34588/149‑1 (MQ=255)                                                                                                                               ctgaggCTCGGCCTGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCAccgc                                           >  1:16363/5‑121 (MQ=255)                                                                                                                               ctgaggCTCGGCCTGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCAccgc                                           <  2:16363/117‑1 (MQ=255)                                                                                                                               ctgaggCTCGGCCTGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCAccgc                                           <  2:11255/117‑1 (MQ=255)                                                                                                                               ctgaggCTCGGCCTGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCAccgc                                           >  1:11255/5‑121 (MQ=255)                                                                                                                                 cTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGa             <  8:92055/149‑1 (MQ=255)                                                                                                                                      tCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGttt         >  6:289661/1‑148 (MQ=255)                                                                                                                                      gcAGGCTGAGGCTCGGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGttt         >  2:203452/2‑148 (MQ=255)                                                                                                                                       cAGGCTGAGGCTCGGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGTTTGa       >  3:253468/1‑149 (MQ=255)                                                                                                                                        aGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGTTTg        <  1:77519/147‑1 (MQ=255)                                                                                                                                          gCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGTTTGATcc    <  4:289306/149‑1 (MQ=255)                                                                                                                                            cGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCCAGGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGTTTGATCCg   >  1:96220/1‑149 (MQ=255)                                                                                                                                            cGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGGGTTTGATCCgg  >  1:78672/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                      CCGACGAACTGATAATCGCTGAGCCTGTGCTCGATGAGCCCGTACAGGAAGAAGCTCAGCCGGAAGAACAGCCCTTTGCCGTGGATGCCGAGCTGCAAGCCGACGCCGACGCCGATTTCGAGGCCATGCTGGCTCAGGCCGAGCCTCAGCCTGCCGTGGACCTGTCTGATTTCGACCTCGACGTCAGCGAGCCGGTTGCCCC‑GGCTGCCCCTGAGGCGCTGGAGCAAGATCCGCTGGATGTCACCGCCGAACTGGCTGCCTTCGACAGCGAACCGGAGTTTGATCCGG  >  NC_002947/2261478‑2261765

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