Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,448	A→G	18.3%	intergenic (‑115/+26)	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,448	0	A	G	18.3%	124.1 / 22.8	60	intergenic (‑115/+26)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base G (8/3);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

CGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACTG  >  NC_002947/1145302‑1145596                                                                                                                                                    |                                                                                                                                                     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                                                                                                                                 >  8:220795/1‑149 (MQ=255)                    ggCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTaa                                                                                                                                 <  6:11557/149‑1 (MQ=255)                           ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAagag                                                                                                                                                 >  4:248329/1‑126 (MQ=255)                           ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAagag                                                                                                                                                 <  3:248329/126‑1 (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                                                                                               <  3:7474/138‑1 (MQ=255)                                                                 cAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                               >  4:7474/1‑138 (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                      >  5:159207/1‑149 (MQ=255)                                                                                                                               ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg                      >  8:195131/1‑149 (MQ=255)                                                                                                                               ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg                      >  4:222268/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                                                              >  3:252193/1‑106 (MQ=255)                                                                                                                                  ttCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                                              <  4:252193/106‑1 (MQ=255)                                                                                                                                     gcg‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgcgc                >  5:28/1‑149 (MQ=255)                                                                                                                                       g‑ggCGTGCCCGCAAAGAGGGCGGTACAGACGTAACGGCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCGCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGGCACAGCCGCCTCAGCCTGGCCGTGGTAGCTGCGCAc              <  4:299686/149‑1 (MQ=255)                                                                                                                                          gcacgCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCAcc             >  3:89569/5‑148 (MQ=255)                                                                                                                                          gCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCAcc             >  3:2974/1‑148 (MQ=255)                                                                                                                                            gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa          <  7:246543/149‑1 (MQ=255)                                                                                                                                            gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa          <  6:101450/149‑1 (MQ=255)                                                                                                                                                ccGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATcc       >  6:182704/1‑148 (MQ=255)                                                                                                                                                ccGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTAGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATcc       >  7:107438/1‑148 (MQ=255)                                                                                                                                                  gCAAAGAGGCCTGTACAGACGTAACGCATCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCGGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAc    <  8:249351/149‑1 (MQ=255)                                                                                                                                                    aaaGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACTg  >  8:56668/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACTG  >  NC_002947/1145302‑1145596

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