Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,793,976	A→G	7.8%	intergenic (‑137/+19)	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,976	0	A	G	7.8%	180.7 / 4.7	64	intergenic (‑137/+19)	pvdI/PP_4222	non‑ribosomal peptide synthetase/protein SyrP
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (3/2);  total (34/30)
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.78e-01

GGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGCCAGGG  >  NC_002947/4793830‑4794120                                                                                                                                                    |                                                                                                                                                 gggCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCA‑TGCGGCAAGCCGCAAg                                                                                                                                                 <  5:71170/148‑1 (MQ=255)                 ggTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGcc                                                                                                                                >  6:125605/1‑149 (MQ=255)                     cTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCa                                                                                                                             >  7:339625/1‑148 (MQ=255)                       tACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAg                                                                                                                          >  2:8760/1‑149 (MQ=255)                         cTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAgc                                                                                                                         <  8:339625/148‑1 (MQ=255)                            aTCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgccg                                                                                                                     >  7:37689/1‑149 (MQ=255)                            aTCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgccg                                                                                                                     >  8:171295/1‑149 (MQ=255)                            aTCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgccg                                                                                                                     >  8:291041/1‑149 (MQ=255)                                     tCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCAc                                                                                                                                              >  3:84668/1‑115 (MQ=255)                                     tCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCAc                                                                                                                                              <  4:84668/115‑1 (MQ=255)                                       ggCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGc                                                                                                                                 <  4:315979/126‑1 (MQ=255)                                       ggCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGc                                                                                                                                 >  3:315979/1‑126 (MQ=255)                                       ggCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGt                                                                                                            <  4:231471/147‑1 (MQ=255)                                            tGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgcc                                                                                                                      >  8:278857/1‑132 (MQ=255)                                            tGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgcc                                                                                                                      <  7:278857/132‑1 (MQ=255)                                              cGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGa                                                                                                    <  5:125605/148‑1 (MQ=255)                                              cGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGa                                                                                                    <  8:37689/148‑1 (MQ=255)                                                gCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGc                                                                                                  >  6:120973/1‑148 (MQ=255)                                                  cAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCAt                                                                                               <  3:283206/149‑1 (MQ=255)                                                      gCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGcc                                                                                           <  2:5060/149‑1 (MQ=255)                                                      gCGTTTCCAGTAGAACGAGCGGAGGGAGGGCAAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGc                                                                                            >  2:66107/1‑148 (MQ=255)                                                                   aaCGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGc                                                                              >  4:166888/1‑149 (MQ=255)                                                                    aCGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCg                                                                             <  1:8760/149‑1 (MQ=255)                                                                     cGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCT‑TGCGGCTTGCCGCAGGcg                                                                                                                                               <  2:341843/82‑2 (MQ=38)                                                                     cGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCT‑TGCGGCTTGCCGCAGGcg                                                                                                                                               >  1:341843/1‑81 (MQ=38)                                                                     cGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGt                                                                            <  5:120973/149‑1 (MQ=255)                                                                      gAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTg                                                                           <  7:299251/149‑1 (MQ=255)                                                                           gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGcc                                                                                                         <  6:88505/114‑1 (MQ=255)                                                                           gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGcc                                                                                                         >  5:88505/1‑114 (MQ=255)                                                                           gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCg                                                                      >  7:247396/1‑149 (MQ=255)                                                                           gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCg                                                                      >  5:94674/1‑149 (MQ=255)                                                                                   gCGAATTTAACGCAGCAGCGCCATGCCCGGACGCTGCGCTCGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACgg                                                              >  4:225430/1‑149 (MQ=255)                                                                                    cGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCa                                                                                    >  8:301464/1‑126 (MQ=255)                                                                                    cGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGt                                                             <  7:171295/149‑1 (MQ=255)                                                                                      aaTTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGc                                                           >  5:1531/1‑149 (MQ=255)                                                                                        tttAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCCGTGCGGCAAGCCGCAAGCAGTCGCCACAGTGGCCACAAAGCGCCGGTTCGCGGTGCCCCTGCGCCATGGCCCCGACCACTTTGCGTGGCCCGGTGAACGGTGcc                                                          >  1:258095/1‑149 (MQ=255)                                                                                           aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGcc                                                                                                 <  1:280436/106‑1 (MQ=255)                                                                                           aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGcc                                                                                                 >  2:280436/1‑106 (MQ=255)                                                                                            aCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGc                                                     <  3:166888/149‑1 (MQ=255)                                                                                            aCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGc                                                     <  3:85515/149‑1 (MQ=255)                                                                                              gcgGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGt                                                   >  5:139000/1‑149 (MQ=255)                                                                                               cgGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGACGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTg                                                  >  2:95404/1‑149 (MQ=255)                                                                                                      gcCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCcgcg                                           <  1:66107/149‑1 (MQ=255)                                                                                                           gCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGCCCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCAt                                      >  5:365154/1‑149 (MQ=255)                                                                                                                 aCGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTc                                <  2:258095/149‑1 (MQ=255)                                                                                                                       cgctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                                       >  5:234103/1‑120 (MQ=255)                                                                                                                       cgctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                                       <  6:234103/120‑1 (MQ=255)                                                                                                                        gctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATc                         <  7:291041/149‑1 (MQ=255)                                                                                                                        gctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATc                         >  7:127623/1‑149 (MQ=255)                                                                                                                           ggCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACTGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAg                       >  4:220656/1‑148 (MQ=255)                                                                                                                             cGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAc                    >  4:286109/1‑149 (MQ=255)                                                                                                                              gAGTGCT‑TGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCg                                                        <  6:284863/112‑1 (MQ=38)                                                                                                                              gAGTGCT‑TGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCg                                                        >  5:284863/1‑112 (MQ=38)                                                                                                                              gAGTGCT‑TGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcg                   >  6:179116/1‑149 (MQ=255)                                                                                                                              gAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcg                   >  4:146120/1‑149 (MQ=255)                                                                                                                                 cgcC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAcgtcg                <  8:327496/149‑1 (MQ=255)                                                                                                                                    cttGCGGCTTGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGtt              >  5:105698/3‑149 (MQ=255)                                                                                                                                        cGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTc           <  8:247396/148‑1 (MQ=255)                                                                                                                                           cAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGc       <  6:139000/149‑1 (MQ=255)                                                                                                                                            aaGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCg                               >  4:258597/1‑124 (MQ=255)                                                                                                                                            aaGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCg                               <  3:258597/124‑1 (MQ=255)                                                                                                                                                cGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGCCAggg  <  3:286109/149‑1 (MQ=255)                                                                                                                                                cGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGCCAggg  <  3:220656/149‑1 (MQ=255)                                                                                                                                                 gCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGCCAggg  >  5:173014/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGCCAGGG  >  NC_002947/4793830‑4794120

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