Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,297,902	A→T	6.2%	A740A (GCT→GCA)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,297,902	0	A	T	6.2%	198.8 / 2.9	64	A740A (GCT→GCA)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base T (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTGCGGG  >  NC_002947/2297762‑2298046                                                                                                                                             |                                                                                                                                                 tcgtTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCaa                                                                                                                                          <  3:87112/149‑1 (MQ=255)       acggCCGTGAGCGCGTCCTCTTGTAGTACGCCGGCCAGGTCGTGCAGGCCTTGCTGCAACAGCTGCTCGTAGTTGGCCAGTGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCGGCAATTGCTg                                                                                                                                    <  7:209473/147‑1 (MQ=255)                 gcgcTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCgg                                                                                                                          <  7:34223/149‑1 (MQ=255)                  cgcTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCggg                                                                                                                         >  3:244456/1‑149 (MQ=255)                       ccTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCcgc                                                                                                                     <  6:42138/148‑1 (MQ=255)                           tGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTg                                                                                                                                             <  1:239596/120‑1 (MQ=255)                           tGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTg                                                                                                                                             >  2:239596/1‑120 (MQ=255)                             gtcgtaCGCCGGCCAGGTCGTTCAGGCCTTGCGGCAACAGCTGCGCGGCGTTGGCCAGGGCTTGGGGCTGTTGC‑CCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGt                                                                                                              <  1:287009/145‑1 (MQ=255)                             gtagtaCGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGt                                                                                                              >  5:92437/1‑149 (MQ=255)                              tagtaCGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGtt                                                                                                             <  6:222114/149‑1 (MQ=255)                                gtaCGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAgg                                                                                                                             <  6:20789/131‑1 (MQ=255)                                gtaCGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAgg                                                                                                                             >  5:20789/1‑131 (MQ=255)                                gtaCGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGAGGGCCGCGTCCCGTTg                                                                                                            >  8:127534/1‑148 (MQ=255)                                gtaCGCCGGCCAGGTCGTGCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGc                                                                                                           <  4:244456/149‑1 (MQ=255)                                 taCGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCa                                                                                                          <  5:113683/149‑1 (MQ=255)                                   cGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAgg                                                                                                        <  7:143852/149‑1 (MQ=255)                                   cGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAg                                                                                                         <  6:92437/148‑1 (MQ=255)                                     ccGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGc                                                                                                      <  8:286643/149‑1 (MQ=255)                                          cAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTcc                                                                                                                 >  5:223156/1‑133 (MQ=255)                                          cAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTcc                                                                                                                 <  6:223156/133‑1 (MQ=255)                                          cAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAg                                                                                                   >  5:37438/1‑147 (MQ=255)                                          cAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCa                                                                                                 <  6:37438/149‑1 (MQ=255)                                           aGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAg                                                                                                >  1:167666/1‑149 (MQ=255)                                               cGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGgcg                                                                                            >  2:323772/1‑149 (MQ=255)                                                    aGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTg                                                                                       <  2:22012/149‑1 (MQ=255)                                                      gCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGc                                                                                     >  6:108421/1‑149 (MQ=255)                                                         ttGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGtt                                                                                  >  2:6103/1‑149 (MQ=255)                                                              gcAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGttt                                                                                 <  5:3719/145‑1 (MQ=255)                                                              gcAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGttt                                                                                 >  6:3719/1‑145 (MQ=255)                                                                aaCAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTc                                                                           <  5:108421/149‑1 (MQ=255)                                                                   aGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTg                                                                        >  2:212362/1‑149 (MQ=255)                                                                    gctgctCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTgg                                                                       >  2:134827/1‑149 (MQ=255)                                                                        ctCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTc                                                                           <  4:278599/141‑1 (MQ=255)                                                                        ctCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTc                                                                           >  3:278599/1‑141 (MQ=255)                                                                        ctCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGa                                                                   >  2:209554/1‑149 (MQ=255)                                                                           gtcgtTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGAttt                                                                <  8:7999/149‑1 (MQ=255)                                                                            tcgtTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATttc                                                               >  6:164126/1‑149 (MQ=255)                                                                              gtTGGCCAGGTCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATttct                                                              <  2:99747/148‑1 (MQ=255)                                                                               ttGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATttctt                                                             >  2:3102/1‑148 (MQ=255)                                                                                 ggCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCa                                                           <  1:3102/148‑1 (MQ=255)                                                                                                  gTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAgctgg                                         >  2:60451/1‑149 (MQ=255)                                                                                                   ttGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAgctggc                                        >  7:59832/1‑149 (MQ=255)                                                                                                   ttGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAgctgg                                         <  7:127534/148‑1 (MQ=255)                                                                                                       tCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTAATCCAGGCGACGCAGCTGGCTGTc                                    >  8:173488/1‑149 (MQ=255)                                                                                                         cAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGa                                  <  8:59832/149‑1 (MQ=255)                                                                                                              ggCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAg                             <  1:6103/149‑1 (MQ=255)                                                                                                              ggCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAg                             <  2:167666/149‑1 (MQ=255)                                                                                                                cgcTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGgc                           <  7:173488/149‑1 (MQ=255)                                                                                                                 gcTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGgcg                          >  6:78068/1‑149 (MQ=255)                                                                                                                 gcTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCAAGGCGACGCAGCTGGCTGTCGAGCCAGccg                          >  5:297809/1‑149 (MQ=255)                                                                                                                    gggCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTa                       >  2:248443/1‑149 (MQ=255)                                                                                                                       ccTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAgg                                                                                                                             <  4:79252/44‑1 (MQ=255)                                                                                                                       ccTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAgg                                                                                                                             >  3:79252/1‑44 (MQ=255)                                                                                                                       ccTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATttct                                                              >  6:88646/1‑107 (MQ=255)                                                                                                                       ccTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATttct                                                              <  5:88646/107‑1 (MQ=255)                                                                                                                           ttgttgGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGc                  <  5:164126/147‑1 (MQ=255)                                                                                                                              ttgGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTgt             >  1:143355/1‑149 (MQ=255)                                                                                                                              ttgGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTgt             >  3:160767/1‑149 (MQ=255)                                                                                                                                gggCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTgtcg           >  4:64486/1‑149 (MQ=255)                                                                                                                                     tcgtGGGCTGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTg      >  4:110623/1‑149 (MQ=255)                                                                                                                                     tcgtGGGCTGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTg      >  4:106072/1‑149 (MQ=255)                                                                                                                                         gggCTGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTGCggg  <  3:110623/149‑1 (MQ=255)                                                                                                                                         gggCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTGCggg  <  1:323772/149‑1 (MQ=255)                                                                                                                                          ggCTGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTGCggg  <  3:106072/148‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                 TCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGATTTCTTCATCCAGGCGACGCAGCTGGCTGTCGAGCCAGGCGCTACGGGCCTTGTCGTCCTGCGGG  >  NC_002947/2297762‑2298046

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: ‑