Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,440	G→A	16.9%	intergenic (‑107/+34)	hemO ← / ← PP_1006	heme oxygenase/heme receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,145,440	0	G	A	16.9%	132.8 / 17.4	59	intergenic (‑107/+34)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base G (27/22);  new base A (7/3);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAA  >  NC_002947/1145295‑1145588                                                                                                                                                   |                                                                                                                                                     gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTg                                                                                                                                                    <  8:96832/148‑1 (MQ=255)      ggCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGt                                                                                                                                                     >  4:131043/1‑142 (MQ=255)      ggCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGt                                                                                                                                                     <  3:131043/142‑1 (MQ=255)        cgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCaaa                                                                                                                                            <  5:77661/149‑1 (MQ=255)        cgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCaaa                                                                                                                                            <  2:165012/149‑1 (MQ=255)             ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAAGAGGc                                                                                                                                       <  5:120509/149‑1 (MQ=255)             ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAAGAGGc                                                                                                                                       >  6:120509/1‑149 (MQ=255)             ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGc                                                                                                                                       >  5:11557/1‑149 (MQ=255)               aGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAAGAGGCCg                                                                                                                                     <  1:232476/149‑1 (MQ=255)                    aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGtc                                                                                                                                  >  7:249351/1‑146 (MQ=255)                        cTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTCGCGTTGCCAGTACCGGCCTCTTCGCG‑AGCGTGCCCGCAAAGAGGCCGGTACAGACg                                                                                                                            >  6:172030/1‑149 (MQ=255)                           ggCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTaa                                                                                                                         <  6:11557/149‑1 (MQ=255)                           ggCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTaa                                                                                                                         >  8:220795/1‑149 (MQ=255)                                  ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAagag                                                                                                                                         <  3:248329/126‑1 (MQ=255)                                  ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAagag                                                                                                                                         >  4:248329/1‑126 (MQ=255)                                    tcATGGGGGAGTACTTGGAAAGAGGTCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCACGCAAAGAGGCCTGTACACACGTAACGCCt                                                                                                                    >  6:50747/1‑145 (MQ=255)                                             aGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                       <  7:220795/149‑1 (MQ=255)                                                 cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                    >  3:293660/1‑148 (MQ=255)                                                    ggAAAGTGGGCTTCTTTAAACAACACTAAACAGCGGGACCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcaca                                                                                                  <  2:173731/147‑1 (MQ=255)                                                     gAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                                >  6:4986/1‑148 (MQ=255)                                                        agagGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                              >  4:110096/1‑147 (MQ=255)                                                        agagGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                              <  3:110096/147‑1 (MQ=255)                                                           gggCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTCGCGTTGCCAGTACCGGCCTCTTCGCG‑AGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCg                                                                                         <  5:172030/149‑1 (MQ=255)                                                             gCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                       <  3:269796/149‑1 (MQ=255)                                                                       aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                             >  6:77674/1‑149 (MQ=255)                                                                       aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTAGCGTTCCCAGTACCGCCCTCTACGCG‑GGCGTGCCCGCAAAGAGGCCCGAACAGACGTAACGCCTCAGAAATCACAAACCACAATGAAACCGAAGTAGCGaccc                                                                             >  2:264143/1‑149 (MQ=255)                                                                        cAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                       >  4:7474/1‑138 (MQ=255)                                                                        cAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                       <  3:7474/138‑1 (MQ=255)                                                                                aaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                     >  8:223544/1‑148 (MQ=255)                                                                                       ggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                                       <  7:15569/139‑1 (MQ=255)                                                                                       ggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                                       >  8:15569/1‑139 (MQ=255)                                                                                        taGCGAGAAATTAAAAAAAGAAGTTATTTTCGATGCCATTACCGGCCTATAAGCA‑GGCGTGCACGCAAAGAGGCCGGTACAGACGTAACTACTAAGATATCAAAAACAACAATGATACCGAAGTTGCGCACCGGCATGGTCAGGCGGTc                                                            >  6:68286/2‑149 (MQ=255)                                                                                        gAGCGAGACAGTAAAAAAAGCCGGTATTTTCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGTCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc                                                            <  5:50747/149‑1 (MQ=255)                                                                                          gCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGa                                                          >  8:246543/1‑149 (MQ=255)                                                                                           cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                         <  1:145985/149‑1 (MQ=255)                                                                                                   tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                                       <  1:152161/127‑1 (MQ=255)                                                                                                   tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                                       >  2:152161/1‑127 (MQ=255)                                                                                                        aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                            <  5:77674/149‑1 (MQ=255)                                                                                                        aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                            <  5:4986/149‑1 (MQ=255)                                                                                                         aaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTggg                                           >  1:34574/1‑149 (MQ=255)                                                                                                         aaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTggg                                           <  4:293660/149‑1 (MQ=255)                                                                                                          aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                                          <  7:223544/149‑1 (MQ=255)                                                                                                            ccGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGTGGGGTCCCCGCAAAGGGGCCGGTACAGCGGTACCGCCGCAGATTTCCCAACCCACTTTGTTACCGAGGTTGCGCCCCGGCATGGTCAGGCGGGCGATTTTGCCTGGCTGGGTc                                         >  1:54805/1‑149 (MQ=255)                                                                                                              ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                                      <  2:34574/149‑1 (MQ=255)                                                                                                                          gCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAccc                          >  7:62161/1‑149 (MQ=255)                                                                                                                                    gCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                 >  2:169261/1‑148 (MQ=255)                                                                                                                                      ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg              >  4:222268/1‑149 (MQ=255)                                                                                                                                      ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg              >  5:159207/1‑149 (MQ=255)                                                                                                                                      ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg              >  8:195131/1‑149 (MQ=255)                                                                                                                                      ctctTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg              >  5:101450/1‑149 (MQ=255)                                                                                                                                       tctTCGCG‑GGCGTGTCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  6:310689/1‑149 (MQ=255)                                                                                                                                       tctTCGCG‑GGCGTGCCCGCAAAGAGGCCGCTACAGACGTAACGCCTCAGATATCCCAACCCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  3:299686/1‑149 (MQ=255)                                                                                                                                       tctTCGCG‑GGCGTGCCAGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGACACAGCCGCCTCACACTCGCCGTCGTAGc             >  7:194241/1‑149 (MQ=255)                                                                                                                                         tttgcg‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTg           >  5:90118/4‑149 (MQ=255)                                                                                                                                         ttCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                                      <  4:252193/106‑1 (MQ=255)                                                                                                                                         ttCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                                      >  3:252193/1‑106 (MQ=255)                                                                                                                                            gcg‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgcgc        >  5:28/1‑149 (MQ=255)                                                                                                                                              g‑ggCGTGCCCGCAAAGAGGGCGGTACAGACGTAACGGCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCGCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGGCACAGCCGCCTCAGCCTGGCCGTGGTAGCTGCGCAc      <  4:299686/149‑1 (MQ=255)                                                                                                                                                 gCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCAcc     >  3:2974/1‑148 (MQ=255)                                                                                                                                                   gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa  <  7:246543/149‑1 (MQ=255)                                                                                                                                                   gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa  <  6:101450/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAA  >  NC_002947/1145295‑1145588

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