Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,793,961	C→T	10.3%	intergenic (‑122/+34)	pvdI ← / ← PP_4222	non‑ribosomal peptide synthetase/protein SyrP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,793,961	0	C	T	10.3%	158.6 / 5.9	58	intergenic (‑122/+34)	pvdI/PP_4222	non‑ribosomal peptide synthetase/protein SyrP
Reads supporting (aligned to +/- strand):  ref base C (29/23);  new base T (4/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGAGTGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTT  >  NC_002947/4793815‑4794109                                                                                                                                                   |                                                                                                                                                      gTTTGAGTGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTg                                                                                                                                                     >  1:32045/1‑149 (MQ=255)  tttGAGTGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGc                                                                                                                                                    >  6:100777/1‑149 (MQ=255)  tttGAGTGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGc                                                                                                                                                    <  2:295300/149‑1 (MQ=255)    tgagtgagTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGGTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCg                                                                                                                                                   <  5:246384/148‑1 (MQ=255)        tgagtCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGTGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAg                                                                                                                                              >  5:242950/1‑149 (MQ=255)                        aTGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGTGCTTGCGGCTTGCC‑GCAGGCGCTCGCca                                                                                                                             <  5:281697/149‑1 (MQ=39)                              ggggTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTg                                                                                                                        <  3:67972/148‑1 (MQ=255)                                         tGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAgcgc                                                                                                            >  4:181000/1‑149 (MQ=255)                                         tGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAgcgc                                                                                                            >  7:111462/1‑149 (MQ=255)                                         tGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGC‑GCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGcgcc                                                                                                           >  4:128790/1‑149 (MQ=255)                                          gATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCGCAAAGcgcc                                                                                                           >  7:252939/1‑149 (MQ=255)                                               gCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTc                                                                                                      >  6:135110/1‑149 (MQ=255)                                                 gctttggCCGCACGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgc                                                                                                    >  4:81233/6‑149 (MQ=255)                                                 aGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgc                                                                                                    >  2:26822/1‑149 (MQ=255)                                                  gCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgcg                                                                                                   <  7:4036/149‑1 (MQ=255)                                                   cTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCgg                                                                                                  >  6:308533/1‑149 (MQ=255)                                                    tCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGt                                                                                                 >  5:294109/1‑149 (MQ=255)                                                      ggCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTg                                                                                                <  5:100777/148‑1 (MQ=255)                                                       gCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGcc                                                                                              >  1:177031/1‑149 (MQ=255)                                                              ggCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgcg                                                                                                   <  6:87059/137‑1 (MQ=255)                                                              ggCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTcgcg                                                                                                   >  5:87059/1‑137 (MQ=255)                                                                ccAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGC‑GCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCAt                                                                                    <  3:128790/149‑1 (MQ=255)                                                                           ccAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGAc                                                                           <  8:111462/148‑1 (MQ=255)                                                                            cAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCa                                                                         >  4:245010/1‑149 (MQ=255)                                                                             aGTAGAACGAGCGGAGGGAGGGCGAATTTAAAGCGGCAGCGCCATGCCCGGACGCTGCGATGGCGAGCGCCTGAGGCAAGCC‑GCAAGCAATAGCCACAGTGGCCTCAAAGCGACGCTTCGCGGTGCCACTGAGCAATGGCCAAGACCAc                                                                        >  6:216791/1‑149 (MQ=255)                                                                               tAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACtt                                                                      <  3:181000/149‑1 (MQ=255)                                                                                  aaCGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGc                                                                   >  1:80982/1‑149 (MQ=255)                                                                                      aGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTgg                                                               <  2:177031/149‑1 (MQ=255)                                                                                      aGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTgg                                                               >  7:71714/1‑149 (MQ=255)                                                                                       gCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGc                                                              >  7:161953/1‑149 (MQ=255)                                                                                        cGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGcc                                                             <  6:294109/149‑1 (MQ=255)                                                                                           agggagggCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCgg                                                          <  6:242950/149‑1 (MQ=255)                                                                                              gagggCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTg                                                        >  1:166242/1‑148 (MQ=255)                                                                                                  gCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATg                                                                                   >  7:131685/1‑117 (MQ=255)                                                                                                  gCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATg                                                                                   <  8:131685/117‑1 (MQ=255)                                                                                                   cGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGt                                                  >  1:241047/1‑149 (MQ=255)                                                                                                    gAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACg                                                    <  5:288829/146‑1 (MQ=255)                                                                                                    gAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACg                                                    >  6:288829/1‑146 (MQ=255)                                                                                                          aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGTCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCggg                                           >  3:63139/1‑149 (MQ=255)                                                                                                          aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                            <  2:80982/148‑1 (MQ=255)                                                                                                          aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                            <  5:308533/148‑1 (MQ=255)                                                                                                              cgGCAGCGCCATGCCCGGACGCTGCGCTGTCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTg                                       <  4:63139/149‑1 (MQ=255)                                                                                                               ggCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTgc                                      <  8:71714/149‑1 (MQ=255)                                                                                                                gCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCCGGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTgc                                      >  1:67442/1‑149 (MQ=255)                                                                                                                gCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACttt                                                                     >  4:168321/1‑117 (MQ=255)                                                                                                                       cATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAg                              >  1:41631/1‑149 (MQ=255)                                                                                                                           cccGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCAtg                          <  2:32045/149‑1 (MQ=255)                                                                                                                                 cgctgcgctgGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCg                    <  5:135110/149‑1 (MQ=255)                                                                                                                                  gctgcgctgGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGa                   <  8:161953/149‑1 (MQ=255)                                                                                                                                  gccactgtgGCGAGTGCTTGCGGCTTGCC‑GCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTc                     >  3:254486/8‑147 (MQ=38)                                                                                                                                    cactgtgGCGAGTGCTTGCGGCTTGCC‑GCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAg                  >  3:242919/6‑148 (MQ=38)                                                                                                                                          ggCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGc           <  3:245010/149‑1 (MQ=255)                                                                                                                                           gCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCa          >  4:94625/1‑149 (MQ=255)                                                                                                                                             gAGTGCTTGCGGCTTGCC‑GCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGcgccgc                                 <  2:116583/124‑1 (MQ=38)                                                                                                                                             gAGTGCTTGCGGCTTGCC‑GCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcg        >  3:27512/1‑149 (MQ=255)                                                                                                                                             gAGTGCTTGCGGCTTGCC‑GCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcg        >  1:2278/1‑149 (MQ=255)                                                                                                                                             gAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGc                                                <  5:283532/109‑1 (MQ=255)                                                                                                                                             gAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGc                                                >  6:283532/1‑109 (MQ=255)                                                                                                                                                   cTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGttt  >  2:134426/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTTTGAGTGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCTGCGGCAAGCC‑GCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTT  >  NC_002947/4793815‑4794109

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