Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,116,922	A→C	7.2%	L158R (CTC→CGC)	PP_2733 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,116,922	0	A	C	7.2%	156.1 / 4.6	56	L158R (CTC→CGC)	PP_2733	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base C (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

GGGCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGAT  >  NC_002947/3116775‑3117070                                                                                                                                                     |                                                                                                                                                       gggCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAg                                                                                                                                                       <  4:226779/149‑1 (MQ=255)   gCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGgc                                                                                                                                                     <  1:127563/149‑1 (MQ=255)          ttGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCTCGCGgcgcggcgc                                                                                                                                              <  3:256481/149‑5 (MQ=255)             tgAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCcagcag                                                                                                                                           >  7:155478/1‑149 (MQ=255)             tgAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCcagca                                                                                                                                            <  6:240787/148‑1 (MQ=255)             tgAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCcag                                                                                                                                              >  5:240787/1‑146 (MQ=255)                    gagCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAgcaggggc‑agg                                                                                                                                   >  8:157266/1‑149 (MQ=255)                    gagCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAgcagggac‑agg                                                                                                                                   >  3:31106/1‑149 (MQ=255)                      gCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGc                                                                                                                                 <  3:317026/149‑1 (MQ=255)                       ccTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGc                                                                                                                                 >  6:236386/1‑148 (MQ=255)                             gCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGcc                                                                                                                          >  4:190434/1‑149 (MQ=255)                              cccAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGcc                                                                                                                          <  7:157266/148‑1 (MQ=255)                               ccAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCa                                                                                                                        >  7:251387/1‑149 (MQ=255)                               ccAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCa                                                                                                                        <  8:155478/149‑1 (MQ=255)                                cAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCa                                                                                                                        >  5:89062/1‑148 (MQ=255)                                      ggTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑Gc                                                                                                                <  3:190434/149‑1 (MQ=255)                                             gtggtgTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAg                                                                                                           >  8:18314/1‑147 (MQ=255)                                              tggtgTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAg                                                                                                        >  7:247474/1‑149 (MQ=255)                                                     tgGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCg                                                                                                 <  3:58605/149‑1 (MQ=255)                                                       gccggccgTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCg                                                                                                 >  4:58605/1‑147 (MQ=255)                                                            ccgTTCACCTGGCGCATGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGAC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCa                                                                                          <  4:31106/149‑1 (MQ=255)                                                              gTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAgc                                                                                        >  3:216858/1‑149 (MQ=255)                                                                 cACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAgcgc                                                                                      <  8:247474/148‑1 (MQ=255)                                                                        cgcAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATg                                                                              <  8:251387/149‑1 (MQ=255)                                                                                cTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggcaaggca                                                                      <  6:89062/149‑1 (MQ=255)                                                                                  gggCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAac                                                                    <  4:216858/149‑1 (MQ=255)                                                                                   ggCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAac                                                                    <  5:236386/148‑1 (MQ=255)                                                                                    gCACGATGTTGGGAGCAGGTGCGTTACTGCTCCGAAGCGATGTCGCCGGTCAGCGCCGCCCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCAATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGCCAAGGCAaca                                                                   >  1:53110/1‑149 (MQ=255)                                                                                     cACGCTGTTGGGAGCAGGTGCGTTTCTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACAcc                                                                 <  1:84617/149‑1 (MQ=255)                                                                                          tGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGc                                                              >  5:131225/1‑147 (MQ=255)                                                                                              gggAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGc                                                          >  7:251541/1‑147 (MQ=255)                                                                                              gggAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGc                                                          <  8:251541/147‑1 (MQ=255)                                                                                               ggAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTg                                                       >  1:62563/1‑149 (MQ=255)                                                                                                gAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTgg                                                      >  8:174047/1‑149 (MQ=255)                                                                                                       tGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGGCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAAGACCTGCCAGCCTGGcgccgc                                                <  2:53110/148‑1 (MQ=255)                                                                                                               tGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAgg                                       >  6:82968/1‑149 (MQ=255)                                                                                                                 cTCCG‑GGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGcc                                     >  5:241156/1‑149 (MQ=255)                                                                                                                    cG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCaga                                  <  6:131225/149‑1 (MQ=255)                                                                                                                     g‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCagag                                 >  3:209330/1‑149 (MQ=255)                                                                                                                     g‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCagag                                 >  6:26554/1‑149 (MQ=255)                                                                                                                     g‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCACCGCCTACATAGGCCagag                                 >  5:355026/1‑149 (MQ=255)                                                                                                                       aGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAgg                                <  5:82968/149‑1 (MQ=255)                                                                                                                            tGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGCGAGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGgc                              >  1:68250/1‑147 (MQ=255)                                                                                                                            tGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGgc                              <  2:68250/146‑1 (MQ=255)                                                                                                                                gCCGGTCAGCGCAGCGCCTCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAAt                       >  8:79555/1‑149 (MQ=255)                                                                                                                                 ccGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGTCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATc                      <  3:235898/149‑1 (MQ=255)                                                                                                                                    gTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGcc                   <  7:174047/149‑1 (MQ=255)                                                                                                                                     tCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGcc                                     <  7:293111/130‑1 (MQ=255)                                                                                                                                     tCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGcc                                     >  8:293111/1‑130 (MQ=255)                                                                                                                                       aGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCACCGCCTACATAGGCCAGAGGGCCGCCAATCGcccc                 <  6:355026/148‑1 (MQ=255)                                                                                                                                         tggcgcgcctcgCGGCGCGGCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCt                                                        >  6:18280/4‑107 (MQ=255)                                                                                                                                         tggcgcgcctcgCGGCGCGGCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCt                                                        <  5:18280/104‑1 (MQ=255)                                                                                                                                          ggcgcgcctcgCGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCcag              >  1:276205/3‑148 (MQ=255)                                                                                                                                           ccgcgccgcgAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCcagca            <  1:312372/149‑1 (MQ=255)                                                                                                                                                  gcgAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACGCGCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggcaagg                                                                        <  1:243457/83‑1 (MQ=255)                                                                                                                                                  gcgAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATggcaagg                                                                        >  2:243457/1‑82 (MQ=255)                                                                                                                                                     aGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGAt  <  7:79555/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GGGCAGGCGTTGTGAAAAAGAGCCTTCTGCCCAGCAAGGTTTGCGTGGTGTGTGGCCGGCCGTTCACCTGGCGCAAGCGCTGGGCACGCTGTTGGGAGCAGGTGCGTTACTGCTCCG‑AGCGATGTCGCCGGTCAGCGCCGCGCCGCGAGGCGCGCCAGCAGGGGC‑AGGGCCAAAGCCCATAACG‑GCATCAGTAGCAGGCCGGTTTCCAGCGCGCCCAATGGCAAGGCAACACCTGCCAGCCTGGCGCCGCCTACATAGGCCAGAGGGCCGCCAATCGCCCCCAGCAGTGCGCCGAT  >  NC_002947/3116775‑3117070

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