Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,832,054	A→T	7.9%	intergenic (+29/+123)	pvdS → / ← pvdY	alternative sigma factor/hydroxyproline acetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,832,054	0	A	T	7.9%	175.1 / 5.5	63	intergenic (+29/+123)	pvdS/pvdY	alternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGAT  >  NC_002947/4831919‑4832196                                                                                                                                         |                                                                                                                                               cGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGctct                                                                                                                                   >  5:254635/1‑149 (MQ=255)  gCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCtctc                                                                                                                                  >  7:55475/1‑149 (MQ=255)    aGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGc                                                                                                                                <  6:303875/149‑1 (MQ=255)          aCATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgc                                                                                                                          >  2:238858/1‑149 (MQ=255)            aTTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgctt                                                                                                                        >  8:225412/1‑149 (MQ=255)             ttGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                                       <  3:27592/149‑1 (MQ=255)             ttGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                                       <  1:292175/149‑1 (MQ=255)                ccAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAg                                                                                                                    >  7:100435/1‑149 (MQ=255)                ccAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAg                                                                                                                    <  2:283632/149‑1 (MQ=37)                 cAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                                       >  4:6972/1‑145 (MQ=255)                 cAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                                       <  3:6972/145‑1 (MQ=255)                 cAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgc                                                                                                                   >  1:52434/1‑149 (MQ=255)                  aaGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcg                                                                                                                  >  1:17370/1‑149 (MQ=255)                   aGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcgc                                                                                                                 >  5:59093/1‑149 (MQ=255)                   aGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcgc                                                                                                                 >  5:46430/1‑149 (MQ=255)                    ggAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCt                                                                                                                <  4:10865/149‑1 (MQ=255)                     gAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATGCACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCt                                                                                                                <  1:238858/149‑1 (MQ=255)                      aaCTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTttgcttg                                                                                                                       >  7:229519/1‑140 (MQ=37)                      aaCTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTttgcttg                                                                                                                       <  8:229519/140‑1 (MQ=37)                              gTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGgtgt                                                                                                      >  4:129789/1‑149 (MQ=255)                                 tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                                   <  7:225412/149‑1 (MQ=255)                                 tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                                   <  6:46430/149‑1 (MQ=255)                                 tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGttt                                                                                                    <  1:168161/148‑1 (MQ=255)                                  cGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTg                                                                                                  >  6:191187/1‑149 (MQ=255)                                             gTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGc                                                                                        <  5:127473/148‑1 (MQ=255)                                             gTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGc                                                                                        <  3:129789/148‑1 (MQ=255)                                                   ttCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAAc                                                                                 >  8:85682/1‑149 (MQ=255)                                                     cATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTc                                                                               >  6:52711/1‑149 (MQ=255)                                                             gcgATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTg                                                                                                  >  6:223467/1‑122 (MQ=255)                                                             gcgATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTg                                                                                                  <  5:223467/122‑1 (MQ=255)                                                              cgATGCGCTGGTGCATTGCCGCAAGACAGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAg                                                                           <  2:289647/144‑1 (MQ=255)                                                                  gcgcTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTACGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCaaaaaa                                                                  >  4:210512/1‑149 (MQ=255)                                                                             ttGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGtt                                                       <  2:17370/149‑1 (MQ=255)                                                                                ccGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATg                                                    >  1:200122/1‑149 (MQ=255)                                                                                    aaGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCg                                                >  5:93286/1‑149 (MQ=255)                                                                                     aGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGt                                               >  4:247431/1‑149 (MQ=255)                                                                                     aGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGt                                               <  5:52711/149‑1 (MQ=255)                                                                                          gCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATACCTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTAttt                                          <  6:254635/149‑1 (MQ=255)                                                                                          gCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTAttt                                          <  2:52434/149‑1 (MQ=255)                                                                                          gCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTAttt                                          <  8:156936/149‑1 (MQ=255)                                                                                          gCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTACGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTAttt                                          <  3:210512/149‑1 (MQ=255)                                                                                            cAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAAt                                                                                    >  6:179418/1‑105 (MQ=255)                                                                                            cAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAAt                                                                                    <  5:179418/105‑1 (MQ=255)                                                                                            cAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCg                                                >  5:248706/1‑141 (MQ=255)                                                                                            cAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCg                                                <  6:248706/141‑1 (MQ=255)                                                                                            cAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                                        >  6:86230/1‑149 (MQ=255)                                                                                             aGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                                        >  5:173482/1‑148 (MQ=255)                                                                                              gTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGa                                      <  8:100435/149‑1 (MQ=255)                                                                                                 gCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTAt                                                     <  6:227993/131‑1 (MQ=255)                                                                                                 gCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTAt                                                     >  5:227993/1‑131 (MQ=255)                                                                                                           aaGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTa                          <  4:261498/148‑1 (MQ=255)                                                                                                               tG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc                     >  2:176439/1‑149 (MQ=255)                                                                                                                   cTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTAc                  <  1:176439/149‑1 (MQ=255)                                                                                                                   cTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTAc                  <  4:186168/149‑1 (MQ=255)                                                                                                                      caagcaaAAGCGAGAGCGATCCAGGTTCGCTCTGGGGTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAAGGCCAATACCTCGCAGAGCAAAAAACTCCGCAAGTTATAGTCGTTATTTGTGAAACCCCACACTAGCCCCTCCGt                >  4:212593/1‑148 (MQ=25)                                                                                                                       aagcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCa              <  3:247431/149‑1 (MQ=255)                                                                                                                           aaaaGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCtt          >  3:155910/1‑149 (MQ=38)                                                                                                                            aaCGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTg         <  5:191187/149‑1 (MQ=255)                                                                                                                             aCGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTg         >  2:124482/1‑148 (MQ=255)                                                                                                                               gCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTAc                  >  4:163264/1‑137 (MQ=255)                                                                                                                               gCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTAc                  <  3:163264/137‑1 (MQ=255)                                                                                                                                ccAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAAc     >  2:43746/1‑149 (MQ=255)                                                                                                                                   gagCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGAt  <  6:59093/149‑1 (MQ=255)                                                                                                                                         |                                                                                                                                               CGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATG‑ACTGCAAGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGAT  >  NC_002947/4831919‑4832196

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