Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,723,275	G→C	6.2%	L97V (CTG→GTG)	htpG ←	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,723,275	0	G	C	6.2%	206.2 / 4.0	65	L97V (CTG→GTG)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base C (2/2);  total (29/36)
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.45e-01

GCGGCTGTACACGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGC  >  NC_002947/4723129‑4723417                                                                                                                                                    |                                                                                                                                               gcggcTGTACACGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGACAGTACCGGACTTGGCGATGGTGCCCAGGt                                                                                                                                              <  6:206332/149‑1 (MQ=255)   ggcTGTACACGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCAACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTATCGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGCATTGTCCATGGTGCACaggagg                                                                                                                                            >  8:128179/1‑149 (MQ=255)    tcTGTACACGTCGACCTTGTCGGCAACGATGTACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGTGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTggg                                                                                                                                           <  1:262353/148‑1 (MQ=255)    gcTGTACACGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTggg                                                                                                                                           <  6:268802/149‑1 (MQ=255)        tACACGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGCACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAAt                                                                                                                                       >  2:271352/1‑149 (MQ=255)         acacGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAAt                                                                                                                                       >  1:83315/1‑148 (MQ=255)          cacGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGa                                                                                                                                     <  3:299323/149‑1 (MQ=255)           acGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGAc                                                                                                                                    >  1:45534/1‑149 (MQ=255)           acGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCACCTGGGCaccatc                                                                                                                                    <  8:279249/149‑6 (MQ=255)             gTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCCCGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGt                                                                                                                                  >  4:157789/1‑149 (MQ=255)                       tCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCt                                                                                                                        <  2:45534/149‑1 (MQ=255)                          gCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCat                                                                                                                     >  4:265821/1‑149 (MQ=255)                           cAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCatg                                                                                                                    <  8:137893/149‑1 (MQ=255)                              cGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCatgccg                                                                                                                 <  8:197521/149‑1 (MQ=255)                                      gcgGATTAGAAGCCCACGCCTAACTGACCTATCAGGTGCTAGTGCTTCTTGTGGTGACGGGTGAGGTTG‑TTCATTAAGTCGGCAGTACCGGACTTGTCGATGTTGCCCATGTGTGCAATGACGTCTTCGCGGCTCATGCCGATGCCGtt                                                                                                         <  7:128179/149‑1 (MQ=255)                                      gcgGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGtt                                                                                                         >  8:37803/1‑149 (MQ=255)                                       cgGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTg                                                                                                        <  3:265821/149‑1 (MQ=255)                                        ggAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGCCAATGACGTCTTCGCGGCTCATGCCGATGCCGTAGt                                                                                                       >  5:99365/1‑149 (MQ=255)                                           gTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCt                                                                                                    <  7:37803/149‑1 (MQ=255)                                            tAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTc                                                                                                   >  3:53388/1‑149 (MQ=255)                                              gAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGa                                                                                                 <  2:83315/149‑1 (MQ=255)                                               aaGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAg                                                                                                >  1:169597/1‑149 (MQ=255)                                                aGCCCACGCCTAACTGACCGATCAGGTTCTAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTTGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGatg                                                                                               <  5:84503/149‑3 (MQ=255)                                                         cGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTa                                                                                      <  2:169597/149‑1 (MQ=255)                                                                 ccGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCGGTTGTCCTCGAGGGTCACGGTACCGGCGTc                                                                              <  1:26559/149‑1 (MQ=255)                                                                  cGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTCTTTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTc                                                                              >  1:224554/1‑149 (MQ=255)                                                                  cGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCCCGAGGGTCACGGTACCGGCGTcc                                                                             <  3:126688/149‑1 (MQ=255)                                                                       aGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGt                                                                         <  6:246814/148‑1 (MQ=255)                                                                        ggTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTc                                                                        >  7:37418/1‑148 (MQ=255)                                                                         gTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCatgccgatgccg                                                                                                           >  8:14034/1‑112 (MQ=255)                                                                         gTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCatgccgatgccg                                                                                                           <  7:14034/112‑1 (MQ=255)                                                                         gTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCg                                                                       >  8:75403/1‑148 (MQ=255)                                                                          tGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGaa                                                                     <  8:37418/149‑1 (MQ=255)                                                                                  cttcttCTGGTCACAGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGACCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGc                                                             >  5:117519/1‑149 (MQ=255)                                                                                          agTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGt                                                                         <  2:202449/128‑1 (MQ=255)                                                                                              aCCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCAAGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAgg                                                  >  5:135503/1‑148 (MQ=255)                                                                                               ccgcagAGCTTC‑TTCATGAAGTCGTCAGTACCGCACTTGGCGATGGTGCCCAGGTTGGCAAAAAAGTCTTCGCGGCTCATACCCAAGCAGTTCTCCTCCACTCACACGGTACCACCGTACTTCTCGAAGCAAAGGCGAAACTTCAgaa                                                 >  8:76645/6‑146 (MQ=255)                                                                                               ccGGTGAGGTTA‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCa                                                    <  7:76645/145‑1 (MQ=255)                                                                                                cGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCa                                               <  7:75403/149‑1 (MQ=255)                                                                                                  gTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAg                                              >  7:4913/1‑148 (MQ=255)                                                                                                      ggTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCg                                         >  5:162239/1‑149 (MQ=255)                                                                                                       gTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCa                                               >  6:251175/1‑142 (MQ=255)                                                                                                       gTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCa                                               <  5:251175/142‑1 (MQ=255)                                                                                                       gTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGc                                        >  6:213602/1‑149 (MQ=255)                                                                                                              cATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTc                                    <  6:135503/147‑1 (MQ=255)                                                                                                               aTGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGt                                                                                       >  7:166540/1‑95 (MQ=255)                                                                                                               aTGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGt                                                                                       <  8:166540/95‑1 (MQ=255)                                                                                                               aTGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCATTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGcc                                       >  5:247106/1‑143 (MQ=255)                                                                                                               aTGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCATTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGcc                                       <  6:247106/143‑1 (MQ=255)                                                                                                                     tCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCa                                                                                           <  7:159345/85‑1 (MQ=255)                                                                                                                     tCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCa                                                                                           >  8:159345/1‑85 (MQ=255)                                                                                                                     tCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGctc                           <  6:162239/149‑1 (MQ=255)                                                                                                                       ggCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAgg                                                              >  8:301743/1‑112 (MQ=255)                                                                                                                       ggCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAgg                                                              <  7:301743/112‑1 (MQ=255)                                                                                                                           gTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCtt                     <  5:213602/149‑1 (MQ=255)                                                                                                                            tACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTg                    >  1:309161/1‑149 (MQ=255)                                                                                                                            tACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGGCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTg                    <  4:53388/149‑1 (MQ=255)                                                                                                                            tACCGGACTTGGCGAGGGTGCCCAGGGGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTg                    <  8:267049/149‑1 (MQ=255)                                                                                                                              ccGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTgg                   >  2:249298/1‑148 (MQ=255)                                                                                                                              ccGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTgg                   <  1:249298/148‑1 (MQ=255)                                                                                                                               cGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGc                  >  2:10807/1‑148 (MQ=255)                                                                                                                                 gACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAg               <  2:224554/149‑1 (MQ=255)                                                                                                                                         gtcattGCCCACCTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCg        >  8:185991/6‑148 (MQ=255)                                                                                                                                         gtcattGCCCACCTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCg        >  8:6418/6‑148 (MQ=255)                                                                                                                                             gTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACg   <  1:271352/149‑1 (MQ=255)                                                                                                                                               gCCCACCTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGc  <  7:185991/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCGGCTGTACACGTCGACCTTGTCGGCAACGATGAACGCGGAGTAGAAGCCCACGCCGAACTGACCGATCAGGTGCGAGTCCTTCTTCTGGTCACCGGTGAGGTTC‑TTCATGAAGTCGGCAGTACCGGACTTGGCGATGGTGCCCAGGTGGGCAATGACGTCTTCGCGGCTCATGCCGATGCCGTTGTCCTCGAGGGTCACGGTACCGGCGTCCTTGTCGAAGCTCAGGCGAATCTTCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGC  >  NC_002947/4723129‑4723417

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑