Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,041,939	A→C	10.8%	E141A (GAG→GCG)	colR →	two‑component system response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,041,939	0	A	C	10.8%	182.2 / 12.8	65	E141A (GAG→GCG)	colR	two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base C (3/4);  total (29/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 = 8.34e-01

TCCGCTCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTG  >  NC_002947/1041792‑1042087                                                                                                                                                     |                                                                                                                                                     tCCGCTCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGGGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAg                                                                                                                                                     >  1:98267/1‑149 (MQ=255)     cTCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTa                                                                                                                                                 >  7:69582/1‑149 (MQ=255)        gggggCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCA‑GGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTaac                                                                                                                                               <  1:199476/147‑3 (MQ=255)            gCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCa                                                                                                                                          <  7:208129/149‑1 (MQ=255)              tGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAgg                                                                                                                                        <  7:182809/149‑1 (MQ=255)              tGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAgg                                                                                                                                        >  3:21790/1‑149 (MQ=255)                acgacTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAggg                                                                                                                                       <  5:271495/148‑1 (MQ=255)                acgacTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAgg                                                                                                                                        <  7:245421/147‑1 (MQ=255)                  gacTATCTGCTCCAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCg                                                                                                                                    <  7:136895/149‑1 (MQ=255)                       ccTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAGGGtatcga                                                                                                                                 <  3:210404/147‑6 (MQ=255)                       ccTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTg                                                                                                                                <  7:168350/148‑1 (MQ=255)                           cTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTg                                                                                                                             >  2:288249/1‑147 (MQ=255)                              aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGCCTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAg                                                                                                                          <  1:91252/148‑1 (MQ=255)                              aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGc                                                                                                                         >  6:61797/1‑148 (MQ=255)                              aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCt                                                                                                                        >  6:197500/1‑149 (MQ=255)                               aGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCt                                                                                                                                 >  1:68540/1‑139 (MQ=255)                               aGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCt                                                                                                                                 <  2:68540/139‑1 (MQ=255)                                gCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCa                                                                                                                      <  2:98267/149‑1 (MQ=255)                                  cATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAAc                                                                                                                    <  4:72379/149‑1 (MQ=255)                                       gCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGt                                                                                                               >  8:190362/1‑149 (MQ=255)                                         ccTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCg                                                                                                             >  2:131551/1‑149 (MQ=255)                                             tCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGc                                                                                                           <  4:69027/147‑1 (MQ=255)                                                    tGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAgctg                                                                                                  >  3:110592/1‑149 (MQ=255)                                                     ggcggcGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAgctgc                                                                                                 >  2:277747/1‑149 (MQ=255)                                                      gcggcgCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAgctgct                                                                                                >  1:143330/1‑149 (MQ=255)                                                        ggcgCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTgg                                                                                              <  1:131551/149‑1 (MQ=255)                                                           gcgcATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGtatcga                                                                                                                                 >  2:109345/1‑106 (MQ=255)                                                           gcgcATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGtatcga                                                                                                                                 <  1:109345/111‑6 (MQ=255)                                                            cgcATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGt                                                                                          <  4:21790/149‑1 (MQ=255)                                                               aTCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACt                                                                                       >  8:158705/1‑149 (MQ=255)                                                                 cGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAC‑GAGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTg                                                                                      >  4:4239/1‑148 (MQ=255)                                                                  gAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGa                                                                                     <  7:190362/148‑1 (MQ=255)                                                                   aaGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTgg                                                                                              >  6:303561/1‑138 (MQ=255)                                                                   aaGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTgg                                                                                              <  5:303561/138‑1 (MQ=255)                                                                           ctgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaaaa                                                                            <  1:288249/148‑1 (MQ=255)                                                                           ctgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaaaa                                                                            <  5:61797/148‑1 (MQ=255)                                                                                ccGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGGAAAGCCCGc                                                                      <  5:207707/149‑1 (MQ=255)                                                                                       cAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACt                                                               <  5:216624/149‑1 (MQ=255)                                                                                              gTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGcgccgc                                                        >  2:50486/1‑149 (MQ=255)                                                                                               tGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaa                                                                              <  2:19186/126‑1 (MQ=255)                                                                                               tGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaa                                                                              >  1:19186/1‑126 (MQ=255)                                                                                                gggCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGa                                                      <  5:197500/149‑1 (MQ=255)                                                                                                 ggcggcGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTTCTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAg                                                     <  1:50486/149‑1 (MQ=255)                                                                                                       gcACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTg                                               <  7:158705/149‑1 (MQ=255)                                                                                                               cAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACt                                                               <  7:121491/125‑1 (MQ=255)                                                                                                               cAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACt                                                               >  8:121491/1‑125 (MQ=255)                                                                                                               cAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCTAGGTGCTGGaagaag                                        <  8:69582/148‑1 (MQ=255)                                                                                                               cAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGgaagaa                                         <  1:277747/147‑1 (MQ=255)                                                                                                                       cTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTgg                                >  2:133651/1‑149 (MQ=255)                                                                                                                         gACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTgggg                              <  3:4239/149‑1 (MQ=255)                                                                                                                               aGCTACGTCCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGcgac                          <  2:143330/147‑1 (MQ=255)                                                                                                                               aGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaaaa                                                                            >  5:24894/1‑97 (MQ=255)                                                                                                                               aGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaaaa                                                                            <  6:24894/97‑1 (MQ=255)                                                                                                                                 cTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTtgtg                                 >  4:193070/1‑138 (MQ=255)                                                                                                                                 cTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTtgtg                                 <  3:193070/138‑1 (MQ=255)                                                                                                                                  tACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTg                     >  3:285671/1‑149 (MQ=255)                                                                                                                                         tCGATACACTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACACGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGa               >  1:101590/1‑148 (MQ=255)                                                                                                                                          cGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGa               >  2:303642/1‑147 (MQ=255)                                                                                                                                           gATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATa             <  4:110592/148‑1 (MQ=255)                                                                                                                                               cccTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGAc        <  1:133651/149‑1 (MQ=255)                                                                                                                                                ccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCCACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCTGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGAGACGACTGCCCTGATAGCGa         >  2:104336/1‑147 (MQ=255)                                                                                                                                                ccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc           >  5:133633/1‑145 (MQ=255)                                                                                                                                                ccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc           <  6:133633/145‑1 (MQ=255)                                                                                                                                                ccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGa         >  4:209929/1‑147 (MQ=255)                                                                                                                                                 cTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAg      >  2:283179/1‑149 (MQ=255)                                                                                                                                                     aGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTg  <  4:285671/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TCCGCTCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTG  >  NC_002947/1041792‑1042087

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