Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,780,122	C→T	8.6%	V2398M (GTG→ATG) ‡	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,122	0	C	T	8.6%	162.9 / 4.8	58	V2398M (GTG→ATG) ‡	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base T (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 = 8.12e-01

CATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCGG  >  NC_002947/4779981‑4780270                                                                                                                                               |                                                                                                                                                      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                                                                                                                                              >  5:60821/1‑122 (MQ=255)                            cTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACGTGTGCGccga                                                                                                                                              <  6:60821/124‑3 (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)                                                                                                                                               cGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCgg  >  7:302570/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      CATCCACCACCACGGCGCTGTGCACACCTTCATGGGCCAGCAGGCGCGATTCGATTTCGCCCAGTTCTATGCGGTAGCCGCGGATCTTCACCTGGAAGTCGGCGCGCCCCACATACTCCAG‑GGTCTGCCCATCGGCGCACACGTGCGCCAGGTCGCCACTGCGGTACAGGCGCTCACCCGGGCGGCTGGCGAACGGATCGGGGACAAAGCTCAATGCCGTGCGCA‑CCGGATCCCCAAGGTAACCGCGCCCGACACCTACGCCGGCCACATGCAACTCGCCGGTGGCCCGG  >  NC_002947/4779981‑4780270

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