Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,094	A→G	12.0%	intergenic (+62/‑68)	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,094	0	A	G	12.0%	162.7 / 12.7	67	intergenic (+62/‑68)	tdcG‑I/PP_5441	L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base G (6/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-02

ACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357946‑358242                                                                                                                                                     |                                                                                                                                                         aCCATGCGCGATACCGGCGCCGATATGCACGACAAGTACACAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCaa                                                                                                                                                          <  4:95122/149‑1 (MQ=255)      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                                                                                                                                                  <  6:10187/149‑1 (MQ=255)         cgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                                  <  1:268555/149‑1 (MQ=255)         cgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                                  <  1:119085/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                                                             >  2:58133/1‑129 (MQ=255)                                                                                                                  ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATggg                                                             <  1:58133/129‑1 (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                      >  5:142955/1‑149 (MQ=37)                                                                                                                                 gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAg                      >  2:169206/1‑149 (MQ=37)                                                                                                                                 gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAg                      >  3:119780/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      >  3:289436/1‑145 (MQ=255)                                                                                                                                                     aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      <  4:289436/145‑1 (MQ=255)                                                                                                                                                     aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGaca  <  8:284317/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                         ACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGC‑‑‑C‑CCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357946‑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: ‑