Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,041,945	C→A	9.5%	T143N (ACC→AAC)	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,945	0	C	A	9.5%	184.8 / 9.2	63	T143N (ACC→AAC)	colR	two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (24/33);  new base A (3/3);  total (27/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 = 6.43e-01

GCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTG  >  NC_002947/1041803‑1042087                                                                                                                                                |                                                                                                                                               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                                                                                                                                 <  2:68540/139‑1 (MQ=255)                    aGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCt                                                                                                                                 >  1:68540/1‑139 (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                                                                                              <  5:303561/138‑1 (MQ=255)                                                        aaGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTgg                                                                                              >  6:303561/1‑138 (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                                                                              >  1:19186/1‑126 (MQ=255)                                                                                    tGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGaa                                                                              <  2:19186/126‑1 (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           <  6:133633/145‑1 (MQ=255)                                                                                                                                     ccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc           >  5:133633/1‑145 (MQ=255)                                                                                                                                     ccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGa         >  4:209929/1‑147 (MQ=255)                                                                                                                                      cTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAg      >  2:283179/1‑149 (MQ=255)                                                                                                                                          aGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTg  <  4:285671/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               GCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAG‑GTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTG  >  NC_002947/1041803‑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: ‑