Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,780,121	A→G	8.6%	V2398A (GTG→GCG) ‡	pvdJ ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,780,121	0	A	G	8.6%	156.2 / 4.5	58	V2398A (GTG→GCG) ‡	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (2/3);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTGCGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCG  >  NC_002947/4779973‑4780269                                                                                                                                                      |                                                                                                                                                      cACCTGCGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGcaca                                                                                                                                                       >  3:154632/1‑149 (MQ=255)         cATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGcc                                                                                                                                               <  1:97763/149‑1 (MQ=255)         cATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGCAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGcc                                                                                                                                               >  2:310826/1‑149 (MQ=255)             caccaccacGGCACTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAGCGGTCTGCCCATCGGCGCACGTGTGCGccgat                                                                                                                                            <  1:231622/149‑4 (MQ=255)               ccaccacGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATTGGCGCACGTGTGCGccgatggg                                                                                                                                         <  8:202055/149‑7 (MQ=255)                        cgcTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCgg                                                                                                                                <  7:202774/149‑1 (MQ=255)                         gcTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATAGGCGCACACGTGCGCCAGGTCGCCACTGCgg                                                                                                                                >  8:213389/1‑148 (MQ=255)                           tgtgCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTCCGGtc                                                                                                                              >  5:348005/1‑147 (MQ=255)                               cacaCCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAgg                                                                                                                          <  1:212420/148‑1 (MQ=255)                                 cacCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGgcgc                                                                                                                       <  3:310534/149‑1 (MQ=255)                                    cTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACGTGTGCGccga                                                                                                                                             <  6:60821/124‑3 (MQ=255)                                    cTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACGTGTGCGccga                                                                                                                                             >  5:60821/1‑122 (MQ=255)                                    cTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCa                                                                                                                    >  2:368337/1‑149 (MQ=255)                                    cTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCa                                                                                                                    >  3:337754/1‑149 (MQ=255)                                         tGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCgg                                                                                                               <  4:154632/149‑1 (MQ=255)                                           ggCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCggg                                                                                                              <  7:273001/148‑1 (MQ=255)                                                gcagGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTg                                                                                                        <  3:66985/149‑1 (MQ=255)                                                 cagGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTgg                                                                                                       <  2:43833/149‑1 (MQ=255)                                                  agGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGCAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGc                                                                                                      <  1:310826/149‑1 (MQ=255)                                                    gcgcGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCg                                                                                                     >  8:342757/1‑148 (MQ=255)                                                       cgattcgattTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACg                                                                                                 <  7:342757/149‑1 (MQ=255)                                                        gattcgattTCGCCCAGTTCTATGCGGTAGCCGCGCATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACg                                                                                                 >  2:12785/1‑148 (MQ=255)                                                           tcgattTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATc                                                                                             >  6:58108/1‑149 (MQ=255)                                                            cgattTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTAGAGGCGCTCACCCGGGCGGCTGGCGAACGGATc                                                                                             <  5:58108/148‑1 (MQ=255)                                                                     ccAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAg                                                                                   >  2:296358/1‑149 (MQ=255)                                                                     ccAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAg                                                                                   >  1:166227/1‑149 (MQ=255)                                                                        gTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTc                                                                                <  5:64649/149‑1 (MQ=255)                                                                            tATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATg                                                                            >  3:333551/1‑149 (MQ=255)                                                                                  gTAGCCGCGGATCTTCACCTGGATGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCAATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACACGGGCGGCTGGCGAACGGATCGGGGACACCGCTCAATGCCCTgc                                                                      >  8:145359/1‑149 (MQ=255)                                                                                  gTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTgc                                                                      >  5:161997/1‑149 (MQ=255)                                                                                        gcgGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCg                                                               >  8:180411/1‑149 (MQ=255)                                                                                                tCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATcccc                                                        >  8:228371/1‑148 (MQ=255)                                                                                                 cACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGTCAAAGCTCAATGCCGTGCGCA‑CCGGATc                                                           <  2:85675/144‑1 (MQ=255)                                                                                                 cACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATcccc                                                        >  5:348354/1‑147 (MQ=255)                                                                                                 cACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATccc                                                         <  7:154120/146‑1 (MQ=255)                                                                                                 cACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATc                                                           <  6:149151/144‑1 (MQ=255)                                                                                                  aCCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGTTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAg                                                     <  1:368337/149‑1 (MQ=255)                                                                                                  aCCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCaa                                                      <  1:12785/148‑1 (MQ=255)                                                                                                  aCCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAg                                                     <  2:41202/149‑1 (MQ=255)                                                                                                    cTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAgg                                                    >  2:33873/1‑148 (MQ=255)                                                                                                        aaGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAAcc                                               <  4:337754/149‑1 (MQ=255)                                                                                                             ggCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGcc                                          <  1:33873/149‑1 (MQ=255)                                                                                                                gcgcCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGa                                       >  5:6000/1‑149 (MQ=255)                                                                                                                    cccACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACAcc                                   <  8:70443/149‑1 (MQ=255)                                                                                                                      cacaTACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACAGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTa                                 >  1:182514/1‑149 (MQ=255)                                                                                                                            cTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCgg                           >  1:214238/1‑149 (MQ=255)                                                                                                                             tCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGc                          >  1:214694/1‑149 (MQ=255)                                                                                                                                    gTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGAGAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTATCCGCGCCCGACACCTACGCCGGCCACATg                    <  1:296358/149‑1 (MQ=255)                                                                                                                                    gTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATg                    <  6:348354/149‑1 (MQ=255)                                                                                                                                     tCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATg                    >  4:285075/1‑148 (MQ=255)                                                                                                                                       tGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCaa                 >  1:343817/1‑149 (MQ=255)                                                                                                                                         cccATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACt               >  8:324625/1‑149 (MQ=255)                                                                                                                                            aTCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGc            >  7:11398/1‑149 (MQ=255)                                                                                                                                             tcgtcgCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGcc           <  4:333551/149‑1 (MQ=255)                                                                                                                                             ctggCGCACGTGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGcc           >  3:14018/3‑149 (MQ=255)                                                                                                                                               ggCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCATTTTGTATGTCGCATTGGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCg          >  1:122654/1‑149 (MQ=255)                                                                                                                                                gcgcACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGt        >  6:301373/1‑149 (MQ=255)                                                                                                                                                 cgcACACGTGCGGCAGGTCGACACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTg       <  6:348005/149‑1 (MQ=255)                                                                                                                                                 cgcACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTg       <  7:213389/149‑1 (MQ=255)                                                                                                                                                      acGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCg  >  8:26948/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CACCTGCGCATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCG  >  NC_002947/4779973‑4780269

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑