Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	496,096	T→A	8.2%	intergenic (‑63/‑33)	PP_0408 ← / → PP_0409	hypothetical protein/two‑component system sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	496,096	0	T	A	8.2%	172.9 / 3.7	60	intergenic (‑63/‑33)	PP_0408/PP_0409	hypothetical protein/two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base A (2/3);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGC  >  NC_002947/495951‑496243                                                                                                                                                  |                                                                                                                                                     tCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCa                                                                                                                                                   >  8:68714/1‑149 (MQ=255)     cATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGc                                                                                                                                               >  7:123918/1‑149 (MQ=255)           cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCCCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac                                                                                                                                         <  8:276512/149‑3 (MQ=255)           cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac                                                                                                                                         <  7:160584/149‑3 (MQ=255)            ggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACTGCCATGCTGCATGGCCGTCttt                                                                                                                                        <  1:192042/149‑1 (MQ=255)                                ggACCGAGGCCATGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGAt                                                                                                                    <  7:80878/149‑1 (MQ=255)                                 gACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                   <  2:26058/149‑1 (MQ=255)                                   ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                   >  5:116423/1‑147 (MQ=255)                                    cGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTCCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAg                                                                                                                >  5:20376/1‑149 (MQ=255)                                       ggCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGc                                                                                                               >  4:113674/1‑147 (MQ=255)                                        gCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAg                                                                                                            <  3:113674/149‑1 (MQ=255)                                         ccAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGt                                                                                                           <  7:68714/149‑1 (MQ=255)                                           aGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                         <  8:123918/149‑1 (MQ=255)                                            ggCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                            <  4:269954/129‑1 (MQ=255)                                            ggCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                            >  3:269954/1‑129 (MQ=255)                                               agcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCg                                                                                                     >  5:39621/1‑149 (MQ=255)                                                gcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGCCAGCCTAAAAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTGGTTGCCAAGTCGATCGTATAGATGAAGCGAGTGCGTCGc                                                                                                    >  6:22387/1‑149 (MQ=255)                                                        ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTgg                                                                                             <  6:39621/148‑1 (MQ=255)                                                          aGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTa                                                                                          >  7:75066/1‑149 (MQ=255)                                                                    gATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCtt                                                                                <  8:75066/149‑1 (MQ=255)                                                                     aTAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTc                                                                                                                                <  1:247358/100‑1 (MQ=255)                                                                     aTAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTc                                                                                                                                >  2:247358/1‑100 (MQ=255)                                                                      tAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTtgt                                                                              >  7:124374/1‑149 (MQ=255)                                                                        aaGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTtgtg                                                                             <  6:116423/148‑1 (MQ=255)                                                                           tgtgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTg                                                                         >  2:305813/1‑149 (MQ=255)                                                                               gACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                            <  8:74457/94‑1 (MQ=255)                                                                               gACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                            >  7:74457/1‑94 (MQ=255)                                                                               gACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGccc                                                                      <  6:64722/148‑1 (MQ=255)                                                                               gACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCt                                                                     <  7:81631/149‑1 (MQ=255)                                                                                aCATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTtgt                                                                              <  3:38264/139‑1 (MQ=255)                                                                                aCATGAGTCAACCCAAGGCAAGAATGAAACCGGCAGCCTAATAGTATTTGCACGTGACGGCCATGCTGCATGCCCGTCTTTGCCAACTCAAGCGTATAGATGAATCGAGTCCGTCGCCTGTTGGTTTTCGGCTGCTtgt                                                                              >  4:38264/1‑139 (MQ=255)                                                                                 cATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCtt                                                                    >  7:217795/1‑148 (MQ=255)                                                                                  aTGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGa                                                                  >  5:164414/1‑149 (MQ=255)                                                                                    gAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑tt                                                               >  7:93323/1‑149 (MQ=255)                                                                                    gAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑tt                                                               >  3:138410/1‑149 (MQ=255)                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTg                                                          >  2:74464/1‑149 (MQ=255)                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTg                                                          <  8:124374/149‑1 (MQ=255)                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTg                                                          >  2:24036/1‑149 (MQ=255)                                                                                                ggCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAg                                                   <  1:74464/149‑1 (MQ=255)                                                                                                 gCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGc                                                  >  2:80118/1‑149 (MQ=255)                                                                                                    aGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTg                                                          <  7:234754/138‑1 (MQ=255)                                                                                                    aGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTg                                                          >  8:234754/1‑138 (MQ=255)                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTg                                                              <  5:171109/118‑1 (MQ=255)                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTg                                                              >  6:171109/1‑118 (MQ=255)                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgt                               >  3:246735/1‑149 (MQ=255)                                                                                                                          aGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCt                         >  3:235044/1‑149 (MQ=255)                                                                                                                          aGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCt                         >  3:230861/1‑149 (MQ=255)                                                                                                                              tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAg                     >  2:85209/1‑149 (MQ=255)                                                                                                                                tGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                    <  4:230861/148‑1 (MQ=255)                                                                                                                                tGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                    <  4:235044/148‑1 (MQ=255)                                                                                                                                  cccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGgc                                           <  7:48352/123‑1 (MQ=255)                                                                                                                                  cccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGgc                                           >  8:48352/1‑123 (MQ=255)                                                                                                                                  cccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGc                          <  1:75237/140‑1 (MQ=255)                                                                                                                                  cccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGc                          >  2:75237/1‑140 (MQ=255)                                                                                                                                    aaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                    <  1:58551/142‑1 (MQ=255)                                                                                                                                    aaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                    >  2:58551/4‑144 (MQ=255)                                                                                                                                       gACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcag            >  2:162524/1‑149 (MQ=255)                                                                                                                                          ggCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcag         <  1:85209/149‑1 (MQ=255)                                                                                                                                              aTGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGt     >  1:85046/1‑149 (MQ=255)                                                                                                                                              aTGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGt     >  1:7084/1‑149 (MQ=255)                                                                                                                                                 ctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGc  <  4:138410/149‑1 (MQ=255)                                                                                                                                                 ctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGc  <  8:93323/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGA‑TTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGC  >  NC_002947/495951‑496243

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