Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	494,451	A→T	6.7%	D71V (GAT→GTT)	djlA →	DnaJ‑like protein DjlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	494,451	0	A	T	6.7%	178.2 / 2.5	59	D71V (GAT→GTT)	djlA	DnaJ‑like protein DjlA
Reads supporting (aligned to +/- strand):  ref base A (33/22);  new base T (2/2);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTGCTCGGTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTGCGGCGCAT  >  NC_002947/494306‑494598                                                                                                                                                   |                                                                                                                                                    ttGCTCGGTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGTTgg                                                                                                                                                  >  4:53353/1‑149 (MQ=255) ttGCTCGGTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGTTgg                                                                                                                                                  >  4:48435/1‑149 (MQ=255)   gCTCGGTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGcc                                                                                                                                                >  6:117532/1‑149 (MQ=255)   gCTCGGTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCAAAGTGCGATGGcc                                                                                                                                                >  8:9016/1‑149 (MQ=255)       ggTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGTTGGCCgg                                                                                                                                              <  3:53353/147‑1 (MQ=255)       ggTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGTTGGCCgg                                                                                                                                              <  3:48435/147‑1 (MQ=255)       ggTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCgg                                                                                                                                              <  5:358776/147‑1 (MQ=255)         tGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGAt                                                                                                                                                    <  7:78359/139‑1 (MQ=255)         tGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGAt                                                                                                                                                    >  8:78359/1‑139 (MQ=255)                cTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGagca                                                                                                                                   <  2:319662/149‑1 (MQ=255)                     gCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCata                                                                                                                              <  1:122837/149‑1 (MQ=255)                     gCAGGCCATCGACCGCAGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCata                                                                                                                              >  5:294946/1‑149 (MQ=255)                           cATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCcag                                                                                                                         >  6:141544/1‑148 (MQ=255)                           cATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCcag                                                                                                                         >  4:37332/1‑148 (MQ=255)                           cATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCca                                                                                                                          <  2:178061/147‑1 (MQ=255)                           cATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCca                                                                                                                          >  1:178061/1‑147 (MQ=255)                           cATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCca                                                                                                                          >  5:352876/1‑147 (MQ=255)                           cATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCca                                                                                                                          <  6:352876/147‑1 (MQ=255)                            aTCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCcag                                                                                                                         >  2:172666/1‑147 (MQ=255)                                 ccgccgGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGgcgc                                                                                                                  <  2:78257/149‑1 (MQ=255)                                      ggCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAg                                                                                                             <  5:141544/149‑1 (MQ=255)                                      ggCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAg                                                                                                             <  1:84733/149‑1 (MQ=255)                                       gCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAg                                                                                                             >  8:131939/1‑148 (MQ=255)                                           cgcTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGAt                                                                                                        >  7:352957/1‑149 (MQ=255)                                            gcTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATg                                                                                                       <  7:311764/149‑1 (MQ=255)                                                 cAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATc                                                                                                                            >  4:286386/1‑123 (MQ=255)                                                 cAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATc                                                                                                                            <  3:286386/123‑1 (MQ=255)                                                       tGGGAAGACATGCGGGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGa                                                                                            >  7:320029/1‑149 (MQ=255)                                                                    gcgAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGTAGATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAgg                                                                                 >  1:92027/1‑148 (MQ=255)                                                                    gcgAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAgg                                                                                 >  5:337797/1‑147 (MQ=255)                                                                      gAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAg                                                                                  <  5:97269/144‑1 (MQ=255)                                                                      gAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAg                                                                                  >  6:97269/1‑144 (MQ=255)                                                                          gCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCgg                                                                         >  7:34170/1‑149 (MQ=255)                                                                            cTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGtt                                                                       >  7:48430/1‑149 (MQ=255)                                                                               ggtgggCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTgcg                                                                    >  4:132178/1‑149 (MQ=255)                                                                                gtgggCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTgcgc                                                                   <  7:9016/149‑1 (MQ=255)                                                                                        ccTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCaa                                                           >  1:335603/1‑149 (MQ=255)                                                                                         cTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAAt                                                          <  7:131939/149‑1 (MQ=255)                                                                                            cATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCg                                                       >  2:308663/1‑149 (MQ=255)                                                                                                 cgcgACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCaa                                                  <  8:48430/149‑1 (MQ=255)                                                                                                      cgacgaGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTg                                             <  8:352957/149‑1 (MQ=255)                                                                                                       gacgaGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTgg                                            <  8:320029/149‑1 (MQ=255)                                                                                                       gacgaGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTgg                                            >  7:167516/1‑149 (MQ=255)                                                                                                          gaGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGc                                           >  3:198780/1‑147 (MQ=255)                                                                                                            gttgttgttCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAgg                                       >  3:102734/1‑149 (MQ=255)                                                                                                                ttgttCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCgg                                    <  1:172666/148‑1 (MQ=255)                                                                                                                  gttCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGc                                  <  3:37332/148‑1 (MQ=255)                                                                                                                              tCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCg                                                                          <  6:3043/97‑1 (MQ=255)                                                                                                                              tCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCg                                                                          >  5:3043/1‑97 (MQ=255)                                                                                                                                    gcCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCggg                >  1:301728/1‑149 (MQ=255)                                                                                                                                     ccTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATTAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCggg                >  4:120247/1‑148 (MQ=255)                                                                                                                                           cAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTgc         >  5:373000/1‑149 (MQ=255)                                                                                                                                           cAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTgc         >  3:222787/1‑149 (MQ=255)                                                                                                                                            aaGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCAc             <  6:287719/144‑1 (MQ=255)                                                                                                                                            aaGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCAc             >  5:287719/1‑144 (MQ=255)                                                                                                                                             aGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTgcgg       <  2:335603/149‑1 (MQ=255)                                                                                                                                                gCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTgcggcg     >  5:107469/1‑148 (MQ=255)                                                                                                                                                  gATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTGCGGCGCAt  >  7:18189/1‑149 (MQ=255)                                                                                                                                                  gATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTGCGGCGCAt  >  1:336076/1‑149 (MQ=255)                                                                                                                                                   aTGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTGCGGCGCAt  >  4:170343/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TTGCTCGGTGCGTTGCTTGGGCAGGCCATCGACCGCCGGCTGCGCTTGCAAGGCTGGGAAGACATGCGCGAACGCCTGGGTGGGCGCCCTGCATTACGCGACGACGAGTTGTTGTTCGT‑GATGCTCGGGCGCCTGGCCAAGTGCGATGGCCGGGTAGCCGAGCAGCATATCCAGCAGGCGCGTCAGGAGATGGTGCGCCTGGACCTGGCCGAGGCAGCCCGGTTGCGCGCCATCAATGCGTTCAACCGTGGCAAGGCGGGCAAGGACCGGCTTGGCGGGCACCTGCGGCGCAT  >  NC_002947/494306‑494598

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