Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,793,968	A→T	9.7%	intergenic (‑129/+27)	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,968	0	A	T	9.7%	166.3 / 8.1	62	intergenic (‑129/+27)	pvdI/PP_4222	non‑ribosomal peptide synthetase/protein SyrP
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base T (4/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

TGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGC  >  NC_002947/4793822‑4794115                                                                                                                                                    |                                                                                                                                                    tgagtCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAg                                                                                                                                                   >  4:120714/1‑149 (MQ=255)  gagtCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGc                                                                                                                                                  >  4:283206/1‑149 (MQ=255)     tCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGc                                                                                                                                               <  3:120714/149‑1 (MQ=255)         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                                                                                                                                         <  4:84668/115‑1 (MQ=255)                                             tCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCAc                                                                                                                                         >  3:84668/1‑115 (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                                                                                                                 <  7:278857/132‑1 (MQ=255)                                                    tGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGcgcc                                                                                                                 >  8:278857/1‑132 (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                                                                                                                                          >  1:341843/1‑81 (MQ=38)                                                                             cGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCT‑TGCGGCTTGCCGCAGGcg                                                                                                                                          <  2:341843/82‑2 (MQ=38)                                                                             cGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGt                                                                       <  5:120973/149‑1 (MQ=255)                                                                              gAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTg                                                                      <  7:299251/149‑1 (MQ=255)                                                                                   gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGcc                                                                                                    >  5:88505/1‑114 (MQ=255)                                                                                   gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGcc                                                                                                    <  6:88505/114‑1 (MQ=255)                                                                                   gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCg                                                                 >  5:94674/1‑149 (MQ=255)                                                                                   gagggaggGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCg                                                                 >  7:247396/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                                                                                            >  2:280436/1‑106 (MQ=255)                                                                                                   aaCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGcc                                                                                            <  1:280436/106‑1 (MQ=255)                                                                                                    aCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGc                                                <  3:85515/149‑1 (MQ=255)                                                                                                    aCGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGc                                                <  3:166888/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                                                  <  6:234103/120‑1 (MQ=255)                                                                                                                               cgctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCgg                                                  >  5:234103/1‑120 (MQ=255)                                                                                                                                gctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATc                    >  7:127623/1‑149 (MQ=255)                                                                                                                                gctgGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATc                    <  7:291041/149‑1 (MQ=255)                                                                                                                                   ggCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACTGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAg                  >  4:220656/1‑148 (MQ=255)                                                                                                                                     cGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCAc               >  4:286109/1‑149 (MQ=255)                                                                                                                                      gAGTGCT‑TGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCg                                                   >  5:284863/1‑112 (MQ=38)                                                                                                                                      gAGTGCT‑TGCGGCTTGCCGCAGGCGCTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCg                                                   <  6:284863/112‑1 (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                          <  3:258597/124‑1 (MQ=255)                                                                                                                                                    aaGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCg                          >  4:258597/1‑124 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGAGTCTTGGGCGTTCATGCTCGGGGTCCTTACTGATCTGCAGCTCGGCCGTGCGGCCAAAGCGTTTCCAGTAGAACGAGCGGAGGGAGGGCGAATTTAACGCGGCAGCGCCATGCCCGGACGCTGCGCTGGCGAGCGCC‑TGCGGCAAGCCGCAAGCACTCGCCACAGTGGCCTCAAAGCGCCGCTTCGCGGTGCCCCTGAGCCATGGCCACGACCACTTTGCGTGGCCCGGTGAACGGTGCCCGGGCGTGCGCCGCGAGCATGTTGTCGAGCATCAGCACGTCGTTTTCCTGC  >  NC_002947/4793822‑4794115

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