Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,376,248	G→A	6.3%	H439H (CAC→CAT)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,376,248	0	G	A	6.3%	189.5 / 2.6	63	H439H (CAC→CAT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/33)
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.71e-01

TCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCA  >  NC_002947/5376100‑5376396                                                                                                                                                      |                                                                                                                                                     tCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACg                                                                                                                                                      >  2:323335/1‑149 (MQ=255)     aCCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGaa                                                                                                                                                  >  8:38951/1‑149 (MQ=255)     aCCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGaa                                                                                                                                                  >  4:98984/1‑149 (MQ=255)       cTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATg                                                                                                                                                <  7:38951/149‑1 (MQ=255)          cGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTc                                                                                                                                             <  2:114883/149‑1 (MQ=255)          cGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTc                                                                                                                                             <  4:53711/149‑1 (MQ=255)              cAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACgg                                                                                                                                         >  2:320563/1‑149 (MQ=255)                   tGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGctc                                                                                                                                    >  1:31623/1‑149 (MQ=255)                     ccGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTg                                                                                                                                  <  5:24111/149‑1 (MQ=255)                             gTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTgtcgt                                                                                                                          >  7:86152/1‑149 (MQ=255)                              tCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTgtcgtc                                                                                                                         <  4:237492/149‑1 (MQ=255)                                gATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGtcgtcg                                                                                                                        <  1:320563/148‑1 (MQ=255)                                 aTGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGTCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACATGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                       >  5:37157/1‑148 (MQ=255)                                 aTGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACATGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                       >  5:40154/1‑148 (MQ=255)                                        aGGTCACTTCGATGTGCGGCACGCCGCGCGGTTCTGGCGGAATGTCCGCCAGGGCGAACTTGCCCAGCGACTTGTTCTTCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCgag                                                                                                               <  8:284608/149‑1 (MQ=255)                                               ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                                                                                                        >  2:280241/1‑149 (MQ=255)                                               ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                                                                                                        >  4:194012/1‑149 (MQ=255)                                                 cGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                                                                                                        <  1:323335/147‑1 (MQ=255)                                                 cGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                                                                                                        <  5:267619/147‑1 (MQ=255)                                                       gCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                       <  5:284966/126‑1 (MQ=255)                                                       gCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                       >  6:284966/1‑126 (MQ=255)                                                        cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGaa                                                                                                                                                  <  2:42675/98‑1 (MQ=255)                                                        cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGaa                                                                                                                                                  >  1:42675/1‑98 (MQ=255)                                                        cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                       >  7:270520/1‑125 (MQ=255)                                                        cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCgg                                                                                                                       <  8:270520/125‑1 (MQ=255)                                                        cGGCACGCCGCGCGGTGCTGGCGGAATGGCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc                                                                                                <  7:159372/148‑1 (MQ=255)                                                            aCGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGt                                                                                           <  1:280241/149‑1 (MQ=255)                                                              gccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCg                                                                                         <  2:31623/149‑1 (MQ=255)                                                                   gcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCGTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtg                                                                                    <  3:98984/149‑1 (MQ=255)                                                                            gCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGATCACTGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCTACTTCTTGGTCGGGATGGTGGTGtt                                                                            <  1:64950/148‑1 (MQ=255)                                                                            gCGGAATGGCAGCCTGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACATTGCA‑GCACGTGAATGATCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCTACTTCTTGGTCGGGATGGTGGTGtt                                                                            <  1:60374/148‑1 (MQ=255)                                                                                 aTGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTtctc                                                                      <  8:23798/149‑1 (MQ=255)                                                                                    tCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGa                                                                    >  8:277904/1‑148 (MQ=255)                                                                                        ccAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTtcttct                                                                       <  8:225829/141‑1 (MQ=255)                                                                                        ccAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTtcttct                                                                       >  7:225829/1‑141 (MQ=255)                                                                                        ccAGGTCGAACTTGCCCAGCGACTTGTTCTGCACAGACTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGACGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAgc                                                               >  7:157892/1‑149 (MQ=255)                                                                                          aGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGACGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAgcgc                                                             >  8:122104/1‑149 (MQ=255)                                                                                           ggTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGcgcg                                                            <  5:320125/149‑1 (MQ=255)                                                                                            gTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTGTTCTCGATCAGCGCgg                                                           >  8:286950/1‑149 (MQ=255)                                                                                            gTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGCTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGCCTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCgg                                                           >  8:54475/1‑149 (MQ=255)                                                                                                aaCTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtca                                                        >  1:80032/1‑149 (MQ=255)                                                                                                 aCTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                                                                                 <  8:210169/122‑1 (MQ=255)                                                                                                 aCTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                                                                                 >  7:210169/1‑122 (MQ=255)                                                                                                   ttGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATAAc                                                    >  4:292062/1‑149 (MQ=255)                                                                                                       ccAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGcgcg                                                            >  1:40226/1‑137 (MQ=255)                                                                                                       ccAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGcgcg                                                            <  2:40226/137‑1 (MQ=255)                                                                                                          gCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                                                                                                      <  5:109155/92‑1 (MQ=255)                                                                                                          gCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                                                                                                      >  6:109155/1‑92 (MQ=255)                                                                                                           cGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCa                                            >  2:229546/1‑149 (MQ=255)                                                                                                               ttGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg                                                                                   >  4:171866/1‑106 (MQ=255)                                                                                                               ttGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg                                                                                   <  3:171866/106‑1 (MQ=255)                                                                                                                  ttCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTc                                     >  6:315975/1‑149 (MQ=255)                                                                                                                  ttCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACATGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGttt                                      <  6:40154/148‑1 (MQ=255)                                                                                                                  ttCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACATGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGttt                                      <  6:37157/148‑1 (MQ=255)                                                                                                                       cgcAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGa                                                                                                    >  5:122914/1‑81 (MQ=255)                                                                                                                       cgcAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGa                                                                                                    <  6:122914/81‑1 (MQ=255)                                                                                                                       cgcAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATAACGCCGCCCATGGTTTCGa                                   <  3:292062/146‑1 (MQ=255)                                                                                                                          aGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtcat                                                       <  5:197276/123‑1 (MQ=255)                                                                                                                          aGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtcat                                                       >  6:197276/1‑123 (MQ=255)                                                                                                                            ccTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                           >  7:257560/1‑149 (MQ=255)                                                                                                                                 ttACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATg                                          >  1:125166/1‑129 (MQ=255)                                                                                                                                 ttACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATg                                          <  2:125166/129‑1 (MQ=255)                                                                                                                                   aCGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                     <  8:86152/148‑1 (MQ=255)                                                                                                                                              gca‑gcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc          <  3:194012/148‑1 (MQ=255)                                                                                                                                                      gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACTTCcagcagca  <  2:80032/149‑1 (MQ=255)                                                                                                                                                      gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcag    >  7:148809/1‑147 (MQ=255)                                                                                                                                                      gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCCCGCTGAccc           >  4:285969/1‑138 (MQ=255)                                                                                                                                                      gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCCCGCTGAccc           <  3:285969/140‑3 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCA‑GCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCA  >  NC_002947/5376100‑5376396

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