Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,096	C→T	11.4%	intergenic (+64/‑66)	tdcG‑I → / → PP_5441	L‑serine dehydratase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	358,096	0	C	T	11.4%	175.3 / 12.2	70	intergenic (+64/‑66)	tdcG‑I/PP_5441	L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (6/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-02

GCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357951‑358242                                                                                                                                                  |                                                                                                                                                       gcgcGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                     >  7:284317/1‑149 (MQ=255)  cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                     <  3:317030/148‑1 (MQ=255)  cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                     <  2:78701/148‑1 (MQ=255)  cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                     <  5:270747/148‑1 (MQ=255)   gcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGt                                                                                                                                                   <  8:171683/149‑1 (MQ=255)    cgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                                  <  1:119085/149‑1 (MQ=255)    cgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                                  <  6:10187/149‑1 (MQ=255)    cgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                                  <  1:268555/149‑1 (MQ=255)       tACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGa                                                                                                                                               <  8:64105/149‑1 (MQ=37)         ccGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGc                                                                                                                                             >  7:271987/1‑149 (MQ=255)            gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCggg                                                                                                                                          <  5:163029/149‑1 (MQ=37)            gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                                          >  6:377432/1‑149 (MQ=255)            gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                                          >  4:248524/1‑149 (MQ=255)               ccGATATGCACGACTAGTACTAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCtt                                                                                                                                       <  4:205228/149‑1 (MQ=255)               ccGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCtt                                                                                                                                       <  5:87090/149‑1 (MQ=255)                  atatGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTAGCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCGCTTCGCGGGCACGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGcc                                                                                                                                    <  1:277579/149‑1 (MQ=255)                      gCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACCGGGACGGTGCAATCCCTGTAGGAGCGGGCTTTCCcgc                                                                                                                                 >  4:54250/1‑148 (MQ=255)                        aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCGGTAGGCGCGGGCTTGCCcgcg                                                                                                                                >  7:46739/1‑147 (MQ=255)                                   aaaGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCCCTCCGAAAGGAACGGTGCAATCCCGGTAGGGGCGGGCTTGCCCCCGAAGAGGCCCTTgt                                                                                                                   >  4:104630/1‑148 (MQ=38)                                    aaGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAg                                                                                                                  >  4:350719/1‑149 (MQ=255)                                     agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGGCATCCCTGTAGGAGCGGGCGTGCGCGCGGAGAGGCCCTTGAgg                                                                                                                 >  5:281663/1‑149 (MQ=38)                                     agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCATCAGGGACAGTGGCATCCATGTAGGAGCGGGGTTGCCCGCGAAGACGCCCTTGAgg                                                                                                                 >  7:161185/1‑149 (MQ=38)                                     agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTCCTGATTTTCTGAGCTGACCTTCTGGGCCGCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAgg                                                                                                                 <  3:75607/149‑1 (MQ=255)                                      gagaCCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGACCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                                                <  1:31600/149‑1 (MQ=255)                                      gagaCCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGCAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                                                <  8:68358/149‑1 (MQ=255)                                         aCCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGGAAGCCCCCTCCGACAGGGACTGTGCGATCGCGGGAGGGGCGGGGTGGCCCGCGCAGAGGGCCGTGAGGCTgg                                                                                                             >  1:353463/1‑148 (MQ=255)                                         aCCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGa                                                                                                             <  5:285606/149‑1 (MQ=255)                                           cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATg                                                                                                           >  7:96926/1‑149 (MQ=255)                                            tCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCGACAGGAACTGTGCAATCCCGGTAGGAGCGGGCTTGGCCGCGAAGAGGCCCTTGAGGCTGATg                                                                                                           >  7:79000/1‑148 (MQ=255)                                              gcccggcggCCTGGGGGTCCGCGCTACCGAGGCCTGATTTTCGGCGCTCCCCTTCGGGGCCGTTTCGCGGGCAAGCCCGCTCCTACCGCACCTGTGCAACCCCTGTCGGTGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGc                                                                                                          <  6:281663/144‑1 (MQ=255)                                                gcggcgGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGTAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGt                                                                                                      <  2:225794/149‑1 (MQ=255)                                                  ggcggcCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCGAGCGGGACTGTGCCATCCCTGTTGGCGCGGGGTTGGCCGCGCAGAGGGCCGTGAGGCGGATTGAGGTg                                                                                                     >  5:202092/1‑148 (MQ=255)                                                     ggcCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGGTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAac                                                                                                 <  5:377432/149‑1 (MQ=255)                                                     ggcCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAac                                                                                                 <  3:248524/149‑1 (MQ=255)                                                                             tGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATg                                                                          <  8:271987/148‑1 (MQ=255)                                                                              gCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGt                                                                        >  6:86751/1‑149 (MQ=255)                                                                               cTGATTTGCTGAGCGCCCCTTCTGGCCCTCTTCGCGGGAACGCCCGCTCCTACCGGAACTGTCCAACCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTg                                                                       <  8:79000/149‑1 (MQ=37)                                                                                          cgCTCACCTTCGGGGCATCTTCCCGGGCACGCCCGCTCCTACCGCACCTGGGCACCCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATg                                                               <  6:202092/145‑1 (MQ=37)                                                                                                        gCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGAAGCAGGTGCAACACAACCGCCGTTGCCGCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCCCCTACCGa                                          >  6:97379/1‑149 (MQ=37)                                                                                                            ctTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGAGGGGCGC‑‑‑C‑CCGCCTACCGAtcgt                                      >  3:228653/1‑149 (MQ=255)                                                                                                             ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATggg                                                             <  1:58133/129‑1 (MQ=255)                                                                                                             ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATggg                                                             >  2:58133/1‑129 (MQ=255)                                                                                                             ttCGCGGGCAAGCCCGCTCCAACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCCCCTATGCCTCTTTATGGTGCGGTGATGGGCCC‑‑‑C‑CCGCCTACCCAacgtc                                     >  8:76949/1‑149 (MQ=38)                                                                                                                 cgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGAt                                 >  4:229765/1‑149 (MQ=255)                                                                                                                      aaGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTcgc                            <  3:229765/149‑1 (MQ=255)                                                                                                                       aGCCCGCTCCTCCAGGAACTGTCCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTcgcg                           <  3:350719/149‑1 (MQ=255)                                                                                                                        gCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGa                          >  8:182942/1‑149 (MQ=37)                                                                                                                        gCCAGCTCCTACATTAAGTGTGCAATCCTGTTAGTATGGGGATTGCACGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATTCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCTCCTACCGATCGTCTGATTTCTCGa                          <  4:228653/149‑1 (MQ=255)                                                                                                                          ccGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                                                >  8:218588/1‑65 (MQ=255)                                                                                                                          ccGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                                                <  7:218588/65‑1 (MQ=255)                                                                                                                           cGCTCCTACAGGAACTGTGCAATCCCTGTTGGAGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCa                       <  2:70175/149‑1 (MQ=255)                                                                                                                           cGCTCCTACAGGAACGGTGCAATCCCTGTACGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCa                       <  3:54250/149‑1 (MQ=255)                                                                                                                           cGCTCCTACAGAAACTGGGAAATCCCTGTAGGCGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCa                       <  4:137723/149‑1 (MQ=255)                                                                                                                            gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACGACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAg                      >  6:289186/1‑149 (MQ=37)                                                                                                                            gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAg                      >  2:169206/1‑149 (MQ=37)                                                                                                                            gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAg                      >  3:119780/1‑149 (MQ=37)                                                                                                                            gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAg                      >  5:142955/1‑149 (MQ=37)                                                                                                                              tCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAgtg                    >  1:301565/1‑149 (MQ=37)                                                                                                                                cTACAGGACCTGTGCAAGCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAgtgt                   <  6:301577/148‑1 (MQ=255)                                                                                                                                 tACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg                 >  3:262715/1‑149 (MQ=255)                                                                                                                                 tACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCACTAGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg                 >  2:57336/1‑149 (MQ=255)                                                                                                                                   cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg                 <  4:262715/147‑1 (MQ=255)                                                                                                                                    aGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGAACGTCGGATGTCGCGATCAGTGTGGGCg              >  2:109771/1‑149 (MQ=255)                                                                                                                                       aaCTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTg           >  1:11782/1‑149 (MQ=255)                                                                                                                                             gCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTTTGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      >  6:333970/1‑148 (MQ=255)                                                                                                                                              cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCAGGCGGAGCAAACCTAGGGTCGGATGTCGCGATCAGTTGGGGCTTTGTCg        >  1:125358/1‑149 (MQ=255)                                                                                                                                              cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGa    >  6:357057/1‑149 (MQ=255)                                                                                                                                                aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      <  4:289436/145‑1 (MQ=255)                                                                                                                                                aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      >  3:289436/1‑145 (MQ=255)                                                                                                                                                aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGaca  <  8:284317/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       GCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357951‑358242

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