Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,090	T→C	12.7%	intergenic (+58/‑72)	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,090	0	T	C	12.7%	144.2 / 13.5	65	intergenic (+58/‑72)	tdcG‑I/PP_5441	L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (6/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-02

CGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTG  >  NC_002947/357943‑358233                                                                                                                                                    |                                                                                                                                                cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                                <  6:76349/149‑1 (MQ=255) cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgt                                                                                                                                                 <  2:107113/148‑1 (MQ=255)  gCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGAGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                                >  5:200065/1‑148 (MQ=255)    aCCATGCGCGATACCGGCGCCGATATGCACGACAAGTACACAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCaa                                                                                                                                             <  4:95122/149‑1 (MQ=255)         gcgcGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                        >  7:284317/1‑149 (MQ=255)          cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                        <  5:270747/148‑1 (MQ=255)          cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                        <  3:317030/148‑1 (MQ=255)          cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                        <  2:78701/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                                                                                                                             >  4:248524/1‑149 (MQ=255)                    gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                             >  6:377432/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)                                                                                                                gCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGAAGCAGGTGCAACACAACCGCCGTTGCCGCTTTATGGTGCGGTGATGGGCGCCCCCCCTACCGa                                 >  6:97379/1‑149 (MQ=37)                                                                                                                    ctTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGAGGGGCGCCCCGCCTACCGAtcgt                             >  3:228653/1‑149 (MQ=255)                                                                                                                     ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATggg                                                <  1:58133/129‑1 (MQ=255)                                                                                                                     ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATggg                                                >  2:58133/1‑129 (MQ=255)                                                                                                                     ttCGCGGGCAAGCCCGCTCCAACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCCCCTATGCCTCTTTATGGTGCGGTGATGGGCCCCCCGCCTACCCAacgtc                            >  8:76949/1‑149 (MQ=38)                                                                                                                         cgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                        >  4:229765/1‑149 (MQ=255)                                                                                                                              aaGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTcgc                   <  3:229765/149‑1 (MQ=255)                                                                                                                               aGCCCGCTCCTCCAGGAACTGTCCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTcgcg                  <  3:350719/149‑1 (MQ=255)                                                                                                                                gCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGa                 >  8:182942/1‑149 (MQ=37)                                                                                                                                gCCAGCTCCTACATTAAGTGTGCAATCCTGTTAGTATGGGGATTGCACGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATTCCTCTTTATGGTGCGGTGATGGGCGCCCCTCCTACCGATCGTCTGATTTCTCGa                 <  4:228653/149‑1 (MQ=255)                                                                                                                                  ccGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                                   <  7:218588/65‑1 (MQ=255)                                                                                                                                  ccGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGc                                                                                                   >  8:218588/1‑65 (MQ=255)                                                                                                                                   cGCTCCTACAGGAACTGTGCAATCCCTGTTGGAGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa              <  2:70175/149‑1 (MQ=255)                                                                                                                                   cGCTCCTACAGGAACGGTGCAATCCCTGTACGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa              <  3:54250/149‑1 (MQ=255)                                                                                                                                   cGCTCCTACAGAAACTGGGAAATCCCTGTAGGCGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa              <  4:137723/149‑1 (MQ=255)                                                                                                                                    gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACGACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAg             >  6:289186/1‑149 (MQ=37)                                                                                                                                    gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAg             >  5:142955/1‑149 (MQ=37)                                                                                                                                    gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAg             >  3:119780/1‑149 (MQ=37)                                                                                                                                    gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAg             >  2:169206/1‑149 (MQ=37)                                                                                                                                      tCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg           >  1:301565/1‑149 (MQ=37)                                                                                                                                        cTACAGGACCTGTGCAAGCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtgt          <  6:301577/148‑1 (MQ=255)                                                                                                                                         tACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg        >  3:262715/1‑149 (MQ=255)                                                                                                                                         tACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCACTAGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg        >  2:57336/1‑149 (MQ=255)                                                                                                                                           cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg        <  4:262715/147‑1 (MQ=255)                                                                                                                                            aGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAACGTCGGATGTCGCGATCAGTGTGGGCg     >  2:109771/1‑149 (MQ=255)                                                                                                                                               aaCTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTg  >  1:11782/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                CGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTG  >  NC_002947/357943‑358233

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