Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,639,021	A→G	6.1%	A130A (GCT→GCC)	czcR‑III ←	response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,639,021	0	A	G	6.1%	189.9 / 2.3	65	A130A (GCT→GCC)	czcR‑III	response regulator
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/32)
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.30e-01

CCACACCAGCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGAC  >  NC_002947/1638874‑1639169                                                                                                                                                    |                                                                                                                                                     ccACACCAGCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAg                                                                                                                                                     >  2:173927/1‑149 (MQ=255)   acacCAGCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGc                                                                                                                                                    >  7:104900/1‑148 (MQ=255)     acCAGCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGGGGCCGGGGTCCAGTTCCAGGTCAGCGAc                                                                                                                                                 <  1:173927/149‑1 (MQ=255)        aGCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGAcgcg                                                                                                                                              >  1:47863/1‑149 (MQ=255)         gCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCTGGGTCCAGTTCCAGGTCAGCGACgcgc                                                                                                                                             <  1:247319/149‑1 (MQ=255)           gagaTGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGc                                                                                                                                                    <  3:110405/140‑1 (MQ=255)           gagaTGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGc                                                                                                                                                    >  4:110405/1‑140 (MQ=255)           gagaTGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAg                                                                                                                                           >  1:306137/1‑149 (MQ=255)                  tCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTg                                                                                                                                    >  6:208969/1‑149 (MQ=255)                     gcgtgcgCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGc                                                                                                                                 >  7:189997/1‑149 (MQ=255)                         gcgcGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTg                                                                                                                               >  7:86869/1‑147 (MQ=255)                          cgcgACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCa                                                                                                                             <  8:104900/148‑1 (MQ=255)                           gcgACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGt                                                                                                                           >  5:100956/1‑149 (MQ=255)                             gACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCGGCGACGCGCAGGGTGCTGGGCTGCAGTTg                                                                                                                         >  6:214564/1‑149 (MQ=255)                               caacaCTTCGCCAGTCCGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCt                                                                                                                       <  7:227358/149‑1 (MQ=255)                                     ttCGCCAGTCTGGCGCATTAGCAGGTGCTGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGt                                                                                                                 >  1:173791/1‑149 (MQ=255)                                     ttCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGt                                                                                                                 >  1:237059/1‑149 (MQ=255)                                           aGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTgcggc                                                                                                           >  7:2774/1‑149 (MQ=255)                                             tCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTCGCGGCCGGGGTCGAGTTCCAGGTCAGCGACGCGCAGGGTGCTCGGCTGCAGTTGCTGGTCGTAgcggcg                                                                                                          >  4:71487/1‑148 (MQ=255)                                                     aTTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCt                                                                                                 >  7:126181/1‑149 (MQ=255)                                                      ttcgCAGGTGCAGCAGGGCGAACTCCTTGGCGGGCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGGGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTg                                                                                                <  8:61887/146‑1 (MQ=255)                                                        aGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCGGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGc                                                                                               <  5:214564/148‑1 (MQ=255)                                                        aGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGc                                                                                               <  7:190247/148‑1 (MQ=255)                                                         gCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCgg                                                                                             <  3:186892/149‑1 (MQ=255)                                                           aGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGAt                                                                                           <  7:69366/149‑1 (MQ=255)                                                             gTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATgc                                                                                         <  2:306137/149‑1 (MQ=255)                                                              tGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATgcg                                                                                        <  2:47863/149‑1 (MQ=255)                                                                 agcagGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATgcgcgc                                                                                     <  8:2774/149‑1 (MQ=255)                                                                  gcagGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATgcgcgc                                                                                     >  5:158954/1‑148 (MQ=255)                                                                   cagGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCa                                                                                   <  2:237059/149‑1 (MQ=255)                                                                      ggCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCa                                                                                                     >  1:102055/1‑128 (MQ=255)                                                                      ggCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCa                                                                                                     <  2:102055/128‑1 (MQ=255)                                                                          atcTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTc                                                                            <  3:111178/147‑1 (MQ=255)                                                                             tCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAgg                                                                         >  8:176758/1‑149 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCGGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAggg                                                                        >  2:105171/1‑148 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCa                                                                                                     >  8:65343/1‑119 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCa                                                                                                     <  7:65343/119‑1 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATgcgc                                                                                       >  7:135264/1‑133 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATgcgc                                                                                       <  8:135264/133‑1 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTc                                                                            >  4:111178/1‑144 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAggg                                                                        >  3:129114/1‑148 (MQ=255)                                                                               cTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAggg                                                                        >  3:269537/1‑148 (MQ=255)                                                                                  ggCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGaa                                                                      <  5:262344/147‑1 (MQ=255)                                                                                  ggCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAAc                                                                     <  7:269080/148‑1 (MQ=255)                                                                                  ggCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAAc                                                                     <  7:176758/148‑1 (MQ=255)                                                                                       tCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAg                                                               <  5:208969/149‑1 (MQ=255)                                                                                         aGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAgg                                                              <  8:189997/148‑1 (MQ=255)                                                                                            tCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCtt                                                          <  8:80556/149‑1 (MQ=255)                                                                                                     tGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCaaa                                                                <  4:231599/134‑1 (MQ=255)                                                                                                     tGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCaaa                                                                >  3:231599/1‑134 (MQ=255)                                                                                                       cccGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTa                                                 <  2:173791/147‑1 (MQ=255)                                                                                                          ggCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCGGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCAt                                            <  1:105171/149‑1 (MQ=255)                                                                                                                gTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCgcgc                                      >  6:301421/1‑149 (MQ=255)                                                                                                                       cgGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCg                                                                                              >  2:225429/1‑86 (MQ=255)                                                                                                                       cgGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCg                                                                                              <  1:225429/86‑1 (MQ=255)                                                                                                                            gcggcCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGc                          >  3:240931/1‑149 (MQ=255)                                                                                                                               gcCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCtt                       >  8:97037/1‑149 (MQ=255)                                                                                                                                    ggTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTAcc                  <  5:301421/149‑1 (MQ=255)                                                                                                                                    ggTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTAcc                  <  7:97037/149‑1 (MQ=255)                                                                                                                                      tCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCg                >  4:292469/1‑149 (MQ=255)                                                                                                                                       ccAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCgg               <  7:289821/149‑1 (MQ=255)                                                                                                                                        cAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGt              >  4:293908/1‑149 (MQ=255)                                                                                                                                              cAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTAcc                  >  2:255727/1‑139 (MQ=255)                                                                                                                                              cAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTAcc                  <  1:255727/139‑1 (MQ=255)                                                                                                                                                 gTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGc     >  1:112778/1‑149 (MQ=255)                                                                                                                                                    aGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGAc  <  4:269537/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCACACCAGCGAGATGATCTGCGTGCGCGACAACACTTCGCCAGTCTGGCGCATTAGCAGGTGCAGCAGGGCGAACTCCTTGGCGGTCAGGTCGATGCGCTGCCCGGCGCGGTAGGCGCGGTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGAC  >  NC_002947/1638874‑1639169

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