Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,951,341	A→C	5.8%	intergenic (‑120/+41)	oguA ← / ← PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,951,341	0	A	C	5.8%	218.0 / 4.6	69	intergenic (‑120/+41)	oguA/PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base C (2/2);  total (36/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 = 5.61e-01

CGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAG  >  NC_002947/2951195‑2951487                                                                                                                                                    |                                                                                                                                                   cGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTaga                                                                                                                                                  <  5:274283/149‑1 (MQ=255)  gTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTagag                                                                                                                                                 <  8:307799/149‑1 (MQ=255)       tGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGTATCGGCGATTACCGATGTCATAGACGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCcgcg                                                                                                                                            >  1:325141/1‑149 (MQ=255)         cAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGATCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGaa                                                                                                                                          <  6:237661/149‑1 (MQ=255)           gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGG‑CTCTAGCCGCCAaga                                                                                                                                        <  4:268794/149‑1 (MQ=255)              tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGc                                                                                                                                     <  8:230254/149‑1 (MQ=255)              tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGc                                                                                                                                     >  3:49758/1‑149 (MQ=255)                  tAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCgg                                                                                                                                  >  7:283008/1‑148 (MQ=255)                  tAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCgg                                                                                                                                  >  6:179180/1‑148 (MQ=255)                  tAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCgg                                                                                                                                  <  1:316983/148‑1 (MQ=255)                  tAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCgg                                                                                                                                  <  6:201311/148‑1 (MQ=255)                   aGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCa                                                                                                                                >  1:50326/1‑149 (MQ=255)                         cATGGTGCAAGCCTGCGCAATACGTATCGGCGATTACCGATGTCATAGACGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGtt                                                                                                                          <  2:325141/149‑1 (MQ=255)                                   gCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCATCCATg                                                                                                                 >  4:296945/1‑148 (MQ=255)                                     cTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGaaa                                                                                                              <  5:179180/149‑1 (MQ=255)                                     cTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGATTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCACTTGGCTG‑CTAGAGCAGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                >  8:220375/1‑147 (MQ=255)                                      tGCGCAATACACATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTAAACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                >  2:89303/1‑146 (MQ=255)                                           aaTACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATca                                                                                                         >  8:242769/1‑148 (MQ=255)                                             tACACATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTAAACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATcaca                                                                                                       <  1:89303/149‑1 (MQ=255)                                                 cATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGcc                                                                                                  <  7:220375/149‑1 (MQ=255)                                                   tCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCcgcg                                                                                                                                            <  1:87281/105‑1 (MQ=255)                                                   tCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCcgcg                                                                                                                                            >  2:87281/1‑105 (MQ=255)                                                     ggCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTcgc                                                                                              <  7:263172/149‑1 (MQ=255)                                                      gCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACggg                                                                                                                            <  6:71085/118‑1 (MQ=255)                                                      gCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACggg                                                                                                                            >  5:71085/1‑118 (MQ=255)                                                       cGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCgg                                                                                            <  8:283008/149‑1 (MQ=255)                                                           tACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgctg                                                                                         >  5:223429/1‑148 (MQ=255)                                                             ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgctgct                                                                                       >  1:25738/1‑148 (MQ=255)                                                             ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTc                                                                                      >  8:10290/1‑149 (MQ=255)                                                              cGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTca                                                                                     >  1:215086/1‑149 (MQ=255)                                                               gATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGCGTTGCAGCCATGAAAGATCACACCGCCTCGCCGCTGCTca                                                                                     >  8:158167/1‑148 (MQ=255)                                                                   tCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACAcc                                                                                 <  2:50326/148‑1 (MQ=255)                                                                     gtAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGG‑CTCTAGCCGCCAag                                                                                                                                         >  4:72706/2‑90 (MQ=37)                                                                     aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgc                                                                                           >  8:77685/1‑136 (MQ=255)                                                                     aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgc                                                                                           <  7:77685/136‑1 (MQ=255)                                                                       aGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCTTGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                              >  3:258422/1‑147 (MQ=255)                                                                        ggCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                            <  7:314826/148‑1 (MQ=255)                                                                         gcgcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAg                                                                           >  6:160823/1‑148 (MQ=255)                                                                            cTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAAc                                                                       <  7:242769/149‑1 (MQ=255)                                                                             tCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAAcc                                                                      >  8:74915/1‑149 (MQ=255)                                                                                gTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCggg                                                                   <  2:215086/149‑1 (MQ=255)                                                                                    tGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAAt                                                               >  1:201606/1‑149 (MQ=255)                                                                                      tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATc                                                              <  2:25738/148‑1 (MQ=255)                                                                                           aagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                        >  5:280071/1‑149 (MQ=255)                                                                                              cTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCggg                                                     >  2:10922/1‑149 (MQ=255)                                                                                                ttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                                    >  1:53494/1‑148 (MQ=255)                                                                                                   tACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCTTGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAg                                                 <  4:258422/148‑1 (MQ=255)                                                                                                    aCTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGAGAGCg                                               >  1:229679/1‑149 (MQ=255)                                                                                                    aCTGATTGTATACAGATCGAAGTGCCCGTGCCGGCATCTTGGAGG‑CTAGATCGGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCTACTCGCGGCTGCTCACACCCTCCAGCAACAGGGCAATCCCCAGCCGGTTAGCg                                               >  1:235227/1‑149 (MQ=255)                                                                                                         ttGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAg                                                         >  2:300779/1‑134 (MQ=255)                                                                                                         ttGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAg                                                         <  1:300779/134‑1 (MQ=255)                                                                                                         ttGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttttt                                          <  2:53494/149‑1 (MQ=255)                                                                                                         ttGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttttt                                          <  4:49758/149‑1 (MQ=255)                                                                                                            tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                        <  5:256090/132‑1 (MQ=255)                                                                                                             ataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCTTTTTTCaa                                       <  7:10290/148‑1 (MQ=255)                                                                                                              taCAGCTCGATGTTACCGTGCCGTCCTCTTGTCGG‑CTAGAGCCTCGAAGAGGCCGGGACGGGTTGCAGCCATGAAAGATCACACCTCCTCGCGGCTGCTAACACCCTCCAGGAACCGGTCAATCCCCAGCTGTTTAGCTTTTTTCAAc                                      <  2:235227/148‑1 (MQ=255)                                                                                                              taCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                      <  7:74915/148‑1 (MQ=255)                                                                                                              taCAGCGCGATGTGCCCGTGCCGGGCGCTTGGCGG‑CTAGAGCCTCGAAGAGGGCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                      <  2:229679/148‑1 (MQ=255)                                                                                                                   cTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATcc                                >  5:239119/1‑149 (MQ=255)                                                                                                                      gAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAAc                                                                       >  3:166269/1‑107 (MQ=255)                                                                                                                      gAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAAc                                                                       <  4:166269/107‑1 (MQ=255)                                                                                                                      gAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                              <  6:223429/148‑1 (MQ=255)                                                                                                                           gCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGcagcag                        >  6:181550/1‑149 (MQ=255)                                                                                                                            cccGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                       >  5:91476/1‑149 (MQ=255)                                                                                                                            cccGTGCCGGCCTCTTCGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGTTCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                       <  7:158167/149‑1 (MQ=255)                                                                                                                                tGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGt                   >  6:58858/1‑149 (MQ=255)                                                                                                                                  ccGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCACACACCCTCCCGCACCCCGGCAATCCCCAGAGCGTTAGCGTTTTTCAACGCCTCCGGCAGCAGTGcctaa                 >  3:235823/1‑145 (MQ=255)                                                                                                                                   cGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAg                >  3:271070/1‑149 (MQ=255)                                                                                                                                            cgcgg‑cTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGc                     <  4:220489/134‑1 (MQ=255)                                                                                                                                            cgcgg‑cTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGc                     >  3:220489/2‑135 (MQ=255)                                                                                                                                                  cTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAg  <  3:296945/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGG‑CTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAG  >  NC_002947/2951195‑2951487

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