Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,902,672	T→C	7.4%	V196A (GTG→GCG) ‡	PP_4310 →	hydantoin racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,902,672	0	T	C	7.4%	190.6 / 5.0	67	V196A (GTG→GCG) ‡	PP_4310	hydantoin racemase
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base C (3/2);  total (34/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 = 7.86e-01

GAACTGGAAGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTA  >  NC_002947/4902527‑4902813                                                                                                                                                  |                                                                                                                                               gaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑Gcc                                                                                                                                            >  6:72445/1‑149 (MQ=255) gaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑Gc                                                                                                                                             <  3:61482/148‑1 (MQ=255)  aactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCg                                                                                                                                           >  1:279332/1‑149 (MQ=255)      ggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCggtgg                                                                                                                                       >  3:214007/1‑149 (MQ=255)       taaGCCGACCCGCTGCGGTCCGTGGAGGCGATTGTCGAGAAGGCCGAACGATCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCTGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCggtggt                                                                                                                                      <  7:159755/148‑1 (MQ=255)        aaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGgtggtg                                                                                                                                     >  3:122104/1‑149 (MQ=255)          gCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGa                                                                                                                                   >  4:40018/1‑149 (MQ=255)          gCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCAC‑GCCggtgcgttg                                                                                                                                   <  5:220299/149‑6 (MQ=255)                ccGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATTgtgt                                                                                                                             <  2:260576/149‑1 (MQ=255)                  gCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGgtgt                                                                                                                             >  7:17028/1‑147 (MQ=255)                  gCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGa                                                                                                                           <  4:139249/149‑1 (MQ=255)                        ggCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAg                                                                                                                     >  8:129445/1‑149 (MQ=255)                          ccGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCg                                                                                                                   >  1:36796/1‑149 (MQ=255)                                   cGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGGTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCAtt                                                                                                          <  6:268/149‑1 (MQ=255)                                    gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTg                                                                                                         >  1:225477/1‑149 (MQ=255)                                    gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTg                                                                                                         >  2:267253/1‑149 (MQ=255)                                    gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTg                                                                                                         >  1:220145/1‑149 (MQ=255)                                    gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTg                                                                                                         <  4:214007/149‑1 (MQ=255)                                               aGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTg                                                                                              <  2:225477/149‑1 (MQ=255)                                               aGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTg                                                                                              <  2:220145/149‑1 (MQ=255)                                               aGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTg                                                                                              >  5:155097/1‑149 (MQ=255)                                               aGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTg                                                                                              <  5:72445/149‑1 (MQ=255)                                                ggCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTCGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGTCGAATCGCTg                                                                                              <  1:267253/149‑1 (MQ=255)                                                   cGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgc                                                                                          <  1:46887/149‑1 (MQ=255)                                                    gAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgcg                                                                                         >  4:94019/1‑149 (MQ=255)                                                      aCGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCa                                                                                       <  2:279332/149‑1 (MQ=255)                                                      aCGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCa                                                                                       <  3:40018/149‑1 (MQ=255)                                                       cgcgCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCAt                                                                                      <  2:295926/149‑1 (MQ=255)                                                       cgcgCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCAt                                                                                      <  1:171351/149‑1 (MQ=255)                                                           ccGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGt                                                                                  <  4:122104/149‑1 (MQ=255)                                                                  cgacgaCAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTg                                                                               >  1:67365/1‑145 (MQ=255)                                                                  cgacgaCAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTg                                                                               <  2:67365/145‑1 (MQ=255)                                                                      gacAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTc                                                                       >  3:298807/1‑149 (MQ=255)                                                                        cAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTcc                                                                      >  1:202747/1‑148 (MQ=255)                                                                         aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGt                                                                                            <  8:246002/125‑1 (MQ=255)                                                                         aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGt                                                                                            >  7:246002/1‑125 (MQ=255)                                                                         aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCa                                                                     >  6:101107/1‑148 (MQ=255)                                                                          aGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                   <  2:202747/149‑1 (MQ=255)                                                                          aGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaa                                                                    <  1:136173/148‑1 (MQ=255)                                                                          aGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaa                                                                    >  2:136173/1‑148 (MQ=255)                                                                            gCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                   <  1:169976/147‑1 (MQ=255)                                                                            gCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                   >  2:169976/1‑147 (MQ=255)                                                                            gCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGt                                                                 <  6:155097/149‑1 (MQ=255)                                                                                       cTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATg                                                      >  8:211364/1‑149 (MQ=255)                                                                                        tGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGc                                                     >  5:266456/1‑149 (MQ=255)                                                                                        tGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGc                                                     >  1:154183/1‑149 (MQ=255)                                                                                          ccTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGc                                                     <  7:129445/147‑1 (MQ=255)                                                                                          ccTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTACCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGc                                                     <  8:17028/147‑1 (MQ=255)                                                                                                          aTGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCgaa                                      <  8:228648/146‑1 (MQ=255)                                                                                                          aTGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCgaa                                      >  7:228648/1‑146 (MQ=255)                                                                                                             gCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAgt                                <  7:211364/149‑1 (MQ=255)                                                                                                               cGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAgtgg                              >  2:131089/1‑149 (MQ=255)                                                                                                                      gACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATTGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTg                       >  7:204863/1‑149 (MQ=255)                                                                                                                        cGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATTGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGc                     <  8:204863/149‑1 (MQ=255)                                                                                                                         gAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGcc                    >  4:94218/1‑149 (MQ=255)                                                                                                                         gAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGcc                    >  1:16808/1‑149 (MQ=255)                                                                                                                               aTTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATgcg              >  7:22255/1‑149 (MQ=255)                                                                                                                               aTTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATgcg              >  2:231704/1‑149 (MQ=255)                                                                                                                                   catccaccACCGGCAC‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTc          >  4:231775/8‑149 (MQ=255)                                                                                                                                     tccaccACCGGCAC‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCgg        >  5:336983/6‑149 (MQ=255)                                                                                                                                     cAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCgg        <  5:101107/149‑1 (MQ=255)                                                                                                                                        accACCGGCAC‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATgc               >  1:195762/3‑139 (MQ=255)                                                                                                                                        accACCGGCAC‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATgc               <  2:195762/137‑1 (MQ=255)                                                                                                                                         gCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc    <  2:36796/149‑1 (MQ=255)                                                                                                                                          cACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCt   <  4:298807/149‑1 (MQ=255)                                                                                                                                           aCCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc    >  1:196669/1‑147 (MQ=255)                                                                                                                                           aCCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTa  <  1:131089/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               GAACTGGAAGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGT‑GCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTA  >  NC_002947/4902527‑4902813

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