Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,547,118	A→G	6.0%	E138G (GAG→GGG)	PP_1357 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,547,118	0	A	G	6.0%	198.8 / 2.4	66	E138G (GAG→GGG)	PP_1357	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTACGCGCTGGACCGCT  >  NC_002947/1546972‑1547267                                                                                                                                                   |                                                                                                                                                      ggcaggcaAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGaa                                                                                                                                                    <  2:475178/150‑1 (MQ=255)    aggcaAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGa                                                                                                                                                 <  2:1095459/150‑1 (MQ=255)       caAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATTAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGacac                                                                                                                                              >  1:419304/1‑150 (MQ=255)       caAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGATAACGacac                                                                                                                                              >  1:1035757/1‑150 (MQ=255)          ggggACCTTGGCCAATGGCGATAAACCGCAGACGGTGGCCCTGGCCATCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACcga                                                                                                                                           <  2:233945/147‑1 (MQ=255)           gTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACcgac                                                                                                                                          <  1:616810/150‑1 (MQ=255)                cTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTa                                                                                                                                     >  2:720254/1‑150 (MQ=255)                 ttGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTAt                                                                                                                                    <  1:720254/150‑1 (MQ=255)                       aaTGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTg                                                                                                                              <  2:804388/150‑1 (MQ=255)                             gATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGAGCTGAACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTc                                                                                                                        >  1:385479/1‑150 (MQ=255)                              aTTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGTTGTGGCGCTGGACTACAAGAGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCa                                                                                                                       <  2:650136/150‑1 (MQ=255)                                 cACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATTAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGa                                                                                                                    <  2:419304/150‑1 (MQ=255)                                 cACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGa                                                                                                                    <  1:515539/150‑1 (MQ=255)                                  aCCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGAc                                                                                                                   >  1:1267074/1‑150 (MQ=255)                                    cGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCa                                                                                                                 >  1:1191847/1‑150 (MQ=255)                                    cGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCa                                                                                                                 <  2:101153/150‑1 (MQ=255)                                        gACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGc                                                                                                             >  2:183819/1‑150 (MQ=255)                                                       cAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGCAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTg                                                                                              >  1:344036/1‑150 (MQ=255)                                                           aTCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTggcgg                                                                                          >  1:816393/1‑150 (MQ=255)                                                             cGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGca                                                                                        <  2:783247/150‑1 (MQ=255)                                                              gAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTgg                                                                                                                             <  2:654155/112‑1 (MQ=255)                                                              gAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTgg                                                                                                                             >  1:654155/1‑112 (MQ=255)                                                                gCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACaa                                                                                     >  1:205417/1‑150 (MQ=255)                                                                   gATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGATGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGc                                                                                  >  2:905971/1‑150 (MQ=255)                                                                       aTGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAg                                                                              <  1:188390/150‑1 (MQ=255)                                                                       aTGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAg                                                                              >  1:973071/1‑150 (MQ=255)                                                                       aTGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAg                                                                              >  1:1222681/1‑150 (MQ=255)                                                                       aTGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTAGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAg                                                                              <  1:1379296/150‑1 (MQ=255)                                                                          aaGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGGGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGt                                                                           >  2:1332827/1‑150 (MQ=255)                                                                          aaGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGGGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGt                                                                           >  2:1333000/1‑150 (MQ=255)                                                                          aaGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGt                                                                           <  1:183819/150‑1 (MQ=255)                                                                           aGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAg                                                                              <  2:346584/146‑1 (MQ=255)                                                                           aGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAg                                                                              >  1:346584/1‑146 (MQ=255)                                                                            gCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGa                                                                         >  2:110646/1‑150 (MQ=255)                                                                             cccAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAg                                                                        >  2:1425309/1‑150 (MQ=255)                                                                              ccAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGt                                                                       <  2:205417/150‑1 (MQ=255)                                                                                 aGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACa                                                                    >  2:1337095/1‑150 (MQ=255)                                                                                        ggTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGGGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCga                                                             <  1:1333000/150‑1 (MQ=255)                                                                                        ggTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGGGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCga                                                             <  1:1332827/150‑1 (MQ=255)                                                                                        ggTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCga                                                             <  2:109044/150‑1 (MQ=255)                                                                                          tCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCgaga                                                           <  1:1425309/150‑1 (MQ=255)                                                                                          tCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGAAACCGACGACCCACACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCgaga                                                           <  1:883761/150‑1 (MQ=255)                                                                                            gAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGAcc                                                         >  1:1302207/1‑150 (MQ=255)                                                                                               gATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAg                                                      >  2:803742/1‑150 (MQ=255)                                                                                               gATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAg                                                      <  1:803742/150‑1 (MQ=255)                                                                                                 tttCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGga                                                    >  1:1368315/1‑150 (MQ=255)                                                                                                    cGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGgacga                                                 <  1:1337095/150‑1 (MQ=255)                                                                                                     gtgtGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGgacgac                                                <  1:834491/150‑1 (MQ=255)                                                                                                      tgtgGACAGGTAATGTCGCCGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACACCGCCGAGGGTGAGTACAACCGCGAGACCAAGGacgacg                                               >  2:205944/1‑150 (MQ=255)                                                                                                      tgtgGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGacgacg                                               >  2:268147/1‑150 (MQ=255)                                                                                                       gtgGAAAGGTAATGTCGACGTGGCGCTGGACTCAACGCGCGCCGAGACCGCCACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGt                                              <  2:1366703/150‑1 (MQ=255)                                                                                                          gAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTcac                                           >  1:1188781/1‑150 (MQ=255)                                                                                                                tAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGACTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTcaccaccac                                     >  1:360096/1‑150 (MQ=255)                                                                                                                    gTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGa                                                                         <  1:926917/110‑1 (MQ=255)                                                                                                                    gTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGa                                                                         >  2:926917/1‑110 (MQ=255)                                                                                                                           tGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGa                          <  2:1191847/150‑1 (MQ=255)                                                                                                                           tGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGACTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGa                          <  2:360096/150‑1 (MQ=255)                                                                                                                               gcTGGACTACAAGCGCGCCGAGAACGACACCGACGACTATGAAATGGGCTTTAACACCACTTCACGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAgcgc                      >  1:739982/1‑150 (MQ=255)                                                                                                                                  ggACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGa                   <  1:110646/150‑1 (MQ=255)                                                                                                                                     cTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACCCCGAGCGGGAGTACAACCGCGAGACCAAGGACGACGGCACCACCACCAACAACTGGAGCGCCGAGTa                >  2:332973/1‑150 (MQ=255)                                                                                                                                       aCAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTAcg              <  2:973071/150‑1 (MQ=255)                                                                                                                                            ttttacGAGGACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGga                                                    >  2:1150680/6‑107 (MQ=255)                                                                                                                                            ttttacGAGGACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGga                                                    <  1:1150680/102‑1 (MQ=255)                                                                                                                                            cgcgcCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTACGCGCTg         <  2:385479/150‑1 (MQ=255)                                                                                                                                              cgcCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTACGCGCTGGa       <  1:268147/150‑1 (MQ=255)                                                                                                                                               gcCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTACGCGCTGGAc      <  2:739982/150‑1 (MQ=255)                                                                                                                                                   agaACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTACGCGCTGGACCGCt  <  1:905971/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACGTCACCACCACCAACAACTGGAGCGCCGAGTACGCGCTGGACCGCT  >  NC_002947/1546972‑1547267

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