Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,812,702	A→G	11.1%	D328G (GAC→GGC)	PP_3368 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,812,702	0	A	G	11.1%	183.3 / 12.4	74	D328G (GAC→GGC)	PP_3368	MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (36/29);  new base G (3/5);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

TGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTGC  >  NC_002947/3812561‑3812849                                                                                                                                               |                                                                                                                                                    tGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGccc                                                                                                                                             <  2:77491/150‑1 (MQ=255)   ggCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCAt                                                                                                                                           >  2:179245/1‑150 (MQ=255)     cgcggaaCTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGa                                                                                                                                         <  1:1043133/144‑1 (MQ=255)        ggTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑GCAATGTCCATGAGgc                                                                                                                                      <  2:372773/150‑1 (MQ=255)         gTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑GCAATGTCCATGAGgcg                                                                                                                                     <  1:949793/150‑1 (MQ=255)           aCTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGcc                                                                                                                                   <  1:1018695/150‑1 (MQ=255)             tGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCTGCCTCATGG‑GCAATGTCCATGAGgcggatc                                                                                                                                 <  2:653199/150‑5 (MQ=255)                   ggtcggtcGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTgg                                                                                                                           <  2:16511/150‑1 (MQ=255)                    gtcggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGc                                                                                                                          >  2:904786/1‑150 (MQ=255)                    gtcggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGc                                                                                                                          >  1:41441/1‑150 (MQ=255)                     tcggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTGGGCGGCACCGCTACAGATCCGCCTCATGG‑ACAGTTCCCATGAGGCGCCTAGCCTGGCa                                                                                                                         <  1:426164/150‑1 (MQ=255)                      cggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAg                                                                                                                        <  2:1006038/150‑1 (MQ=255)                         tcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACa                                                                                                                                                   >  2:961695/1‑120 (MQ=255)                         tcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACa                                                                                                                                                   <  1:961695/120‑1 (MQ=255)                         tcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAgcgg                                                                                                                     >  1:1099865/1‑150 (MQ=255)                            ccGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGt                                                                                                                  <  1:707636/150‑1 (MQ=255)                              gTGCTGCTGTTCTTTACCTTCACCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTcc                                                                                                                <  2:838874/150‑1 (MQ=255)                                     tgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATg                                                                                                         >  2:910013/1‑150 (MQ=255)                                     tgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATg                                                                                                         >  2:908314/1‑150 (MQ=255)                                     tgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATg                                                                                                         <  1:465354/150‑1 (MQ=255)                                      gTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgc                                                                                                        <  2:358807/150‑1 (MQ=255)                                      gTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgc                                                                                                        <  2:1098488/150‑1 (MQ=255)                                       ttCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGCCATCGGCCTGGCGCTCACTATGCTCGCCTTGGCGGCACCGCTACAGATCCCCCCCATGG‑ACATTGCCCCTGCGGCGCCTAGCCTGGCAGCCGCGTCCAAccccgcg                                                                                                       >  2:427526/1‑145 (MQ=255)                                           ttACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGt                                                                                                   >  2:147728/1‑150 (MQ=255)                                               cTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCaa                                                                                               >  2:603562/1‑150 (MQ=255)                                                ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                              >  2:624315/1‑150 (MQ=255)                                                ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                              >  1:813839/1‑150 (MQ=255)                                                ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATAGACTTGGCGGCACCGCTACAGATCGGAAGAGCAC‑ACGTTTGCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                              >  1:252464/1‑150 (MQ=255)                                                   gccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTg                                                                                           >  2:568663/1‑150 (MQ=255)                                                    ccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑GCAATGTCCATGAGgc                                                                                                                                      <  1:1006661/106‑1 (MQ=255)                                                    ccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑GCAATGTCCATGAGgc                                                                                                                                      >  2:1006661/1‑106 (MQ=255)                                                     cgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCAGCACCACTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGc                                                                                         <  1:7283/150‑1 (MQ=255)                                                       cccAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCa                                                                                       >  2:482217/1‑150 (MQ=255)                                                               cTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTg                                                                               >  1:140397/1‑150 (MQ=255)                                                                tgtgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTgg                                                                              <  2:497844/150‑1 (MQ=255)                                                                  tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTgggg                                                                            >  2:152590/1‑150 (MQ=255)                                                                    gAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGA‑ACATTGCCCATGAGGCGAGTAGGCTCGGTCTGACCTCCAACCAGGCGGCGTTCAACCTGGCCAATGCCCTGGGGcc                                                                          >  1:1097925/1‑150 (MQ=255)                                                                      gCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGt                                                                        <  1:482217/150‑1 (MQ=255)                                                                              ctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGa                                                                <  2:1099865/150‑1 (MQ=255)                                                                                gggCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCACATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAgg                                                              >  2:401045/1‑150 (MQ=255)                                                                                   cATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCAt                                                           <  2:813839/150‑1 (MQ=255)                                                                                       ggCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCg                                                       >  1:817445/1‑150 (MQ=255)                                                                                        gCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGa                                                      >  1:693364/1‑150 (MQ=255)                                                                                        gCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGa                                                      >  1:495908/1‑150 (MQ=255)                                                                                          cTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCgctg                                                    <  2:211439/150‑4 (MQ=255)                                                                                           tggtggGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCa                                                   >  1:887366/1‑150 (MQ=255)                                                                                               gggCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGCGCCTTGGCTGGGAGGCATGGCGATCACTGc                                               <  1:261148/150‑1 (MQ=255)                                                                                                 gCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCg                                             >  1:281540/1‑150 (MQ=255)                                                                                                   aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGc                                           >  2:483531/1‑150 (MQ=255)                                                                                                    ccATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCa                                          <  1:568663/150‑1 (MQ=255)                                                                                                     cATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCAt                                         >  2:51863/1‑150 (MQ=255)                                                                                                      aTGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATg                                        <  2:495908/150‑1 (MQ=255)                                                                                                        gATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATggg                                      <  1:904786/150‑1 (MQ=255)                                                                                                        gATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATggg                                      >  1:883878/1‑150 (MQ=255)                                                                                                          tCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTgg                                                                                          >  2:71209/1‑96 (MQ=255)                                                                                                          tCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTgg                                                                                          <  1:71209/96‑1 (MQ=255)                                                                                                              cTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAc                                <  1:1006054/150‑1 (MQ=255)                                                                                                              cTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAc                                >  1:183929/1‑150 (MQ=255)                                                                                                              cTTGGCGGCACCGCTACAGA‑CCCCCTTGTGGTTCATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAc                                <  1:840678/150‑1 (MQ=255)                                                                                                               ttGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAcc                               >  1:1100139/1‑150 (MQ=255)                                                                                                               ttGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAcc                               >  2:756464/1‑150 (MQ=255)                                                                                                                  gcggcACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGc                            <  2:1023253/150‑1 (MQ=255)                                                                                                                       aCCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCgg                       >  2:381364/1‑150 (MQ=255)                                                                                                                       aCCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGACGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCgg                       >  1:277491/1‑150 (MQ=255)                                                                                                                        ccGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                      <  2:252464/150‑1 (MQ=255)                                                                                                                        ccGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                      <  1:147728/150‑1 (MQ=255)                                                                                                                             aCAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACAt                 >  1:306652/1‑150 (MQ=255)                                                                                                                               aGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCg               <  2:41441/150‑1 (MQ=255)                                                                                                                                    cGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgcc          >  2:580051/1‑150 (MQ=255)                                                                                                                                      ccTCATGG‑GCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTgg                                  <  2:724352/124‑1 (MQ=255)                                                                                                                                      ccTCATGG‑GCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTgg                                  >  1:724352/1‑124 (MQ=255)                                                                                                                                        tCATGG‑GCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCACTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCGCCGGATACATCGGCGCCGCGa      >  1:862812/1‑150 (MQ=255)                                                                                                                                           tGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTg   >  1:420348/1‑150 (MQ=255)                                                                                                                                            gg‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTgc  <  1:483531/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    TGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGG‑ACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTGC  >  NC_002947/3812561‑3812849

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