Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,448	A→G	7.4%	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	7.4%	188.5 / 4.8	67	intergenic (‑115/+26)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base G (2/3);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAC  >  NC_002947/1145301‑1145594                                                                                                                                                    |                                                                                                                                                   gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTGGCGGGCGTGCCCGCa                                                                                                                                                    >  6:28035/1‑148 (MQ=255) gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCa                                                                                                                                                    >  1:36242/1‑148 (MQ=255)  cgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCaaa                                                                                                                                                  <  7:330461/149‑1 (MQ=255)    cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAaga                                                                                                                                                <  7:75481/149‑1 (MQ=255)             cAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTa                                                                                                                                        >  7:270243/1‑148 (MQ=255)                ggCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGTGCGTGCCCGCAAAGAGGCCGGTACAGa                                                                                                                                    >  5:92724/1‑149 (MQ=255)                 gCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGAc                                                                                                                                   <  7:249024/149‑1 (MQ=255)                   tGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAagag                                                                                                                                               <  6:113934/135‑1 (MQ=255)                   tGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAagag                                                                                                                                               >  5:113934/1‑135 (MQ=255)                          tcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCt                                                                                                                          >  7:27450/1‑149 (MQ=255)                           cggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACg                                                                                                                                  <  7:280295/140‑1 (MQ=255)                           cggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACg                                                                                                                                  >  8:280295/1‑140 (MQ=255)                              tcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAagag                                                                                                                                               <  4:95107/124‑1 (MQ=255)                              tcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAagag                                                                                                                                               >  3:95107/1‑124 (MQ=255)                              tcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGAGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCt                                                                                                                          >  7:304924/1‑145 (MQ=255)                                   ggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATcc                                                                                                                 <  6:310491/149‑1 (MQ=255)                                          ccTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                          >  1:356488/1‑149 (MQ=255)                                           cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATc                                                                                                                  >  4:210672/1‑140 (MQ=255)                                           cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATc                                                                                                                  <  3:210672/140‑1 (MQ=255)                                           cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                          >  1:177134/1‑148 (MQ=255)                                             tGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                       >  1:329911/1‑149 (MQ=255)                                              ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTGGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                       <  5:28035/148‑1 (MQ=255)                                              ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                                <  8:331975/107‑1 (MQ=255)                                              ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                                >  7:331975/1‑107 (MQ=255)                                                     gggCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGGCCGCAAAGAGGCCGGTACAGACGTAACGCCGCAGCTATCCCAAACCACATTgcgacc                                                                                                >  6:59555/1‑143 (MQ=255)                                                            tAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                        >  3:2292/1‑149 (MQ=255)                                                              gATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAAc                                                                                                            >  7:235804/1‑127 (MQ=255)                                                                    aGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCgg                                                                                 >  7:47042/1‑148 (MQ=255)                                                                     gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCa                                                                               >  5:21996/1‑149 (MQ=255)                                                                       cGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                                             <  2:356488/149‑1 (MQ=255)                                                                         aaaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                           <  4:224194/149‑1 (MQ=255)                                                                         aaaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                           >  6:89319/1‑149 (MQ=255)                                                                          aaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTc                                                                          <  2:27415/149‑1 (MQ=255)                                                                              gCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTTACAGTACCGGCATCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGc                                                                      >  4:325399/1‑149 (MQ=255)                                                                                gggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCgg                                                                    >  5:288740/1‑149 (MQ=255)                                                                                 ggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGt                                                                   <  6:21996/149‑1 (MQ=255)                                                                                 ggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGt                                                                   <  3:259271/149‑1 (MQ=255)                                                                                     cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCTGGCGGTCGAt                                                               <  6:92724/149‑1 (MQ=255)                                                                                     cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                               <  2:177134/149‑1 (MQ=255)                                                                                        gaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                                            <  8:270243/149‑1 (MQ=255)                                                                                         aCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                           <  8:47042/149‑1 (MQ=255)                                                                                           aGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGCTATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggc                                                         >  8:256742/1‑149 (MQ=255)                                                                                             tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGGCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctg                                                       >  1:314937/1‑149 (MQ=255)                                                                                             tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctg                                                       >  2:280621/1‑149 (MQ=255)                                                                                             tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctg                                                       >  7:109873/1‑149 (MQ=255)                                                                                              aaaaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg                                                                     >  4:21618/1‑134 (MQ=255)                                                                                              aaaaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg                                                                     <  3:21618/134‑1 (MQ=255)                                                                                               aaaacaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                             <  6:356500/93‑1 (MQ=255)                                                                                               aaaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGCGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                             >  5:356500/1‑93 (MQ=255)                                                                                                  aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                                  <  8:27450/149‑1 (MQ=255)                                                                                                    aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAg                                                                        <  7:112449/125‑1 (MQ=255)                                                                                                    aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAg                                                                        >  8:112449/1‑125 (MQ=255)                                                                                                         gTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAg                                           >  4:24926/1‑149 (MQ=255)                                                                                                                         tACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                       <  6:142693/73‑1 (MQ=255)                                                                                                                           ccGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcg                          <  5:89319/148‑1 (MQ=255)                                                                                                                              gCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTggg                                                 >  7:162362/1‑122 (MQ=255)                                                                                                                              gCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTggg                                                 <  8:162362/122‑1 (MQ=255)                                                                                                                              gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAGTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGGCGATGTTGGCTGGCTGGGTCACCGCCGCCTCACCCTGGCcgtcg                       >  6:80773/1‑148 (MQ=255)                                                                                                                                     cgcgGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCt                  <  2:36242/146‑1 (MQ=255)                                                                                                                                      gcgGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgcgc              >  4:70943/1‑149 (MQ=255)                                                                                                                                          gCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggct                                                        <  3:336444/103‑1 (MQ=255)                                                                                                                                          gCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggct                                                        >  4:336444/1‑103 (MQ=255)                                                                                                                                            acgCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGa         >  8:152696/3‑148 (MQ=255)                                                                                                                                             tGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAAt       >  2:10496/1‑149 (MQ=255)                                                                                                                                               cccGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATcc     >  4:288944/1‑149 (MQ=255)                                                                                                                                                ccGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATcc     >  3:172013/1‑148 (MQ=255)                                                                                                                                                 cGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCa   <  3:288944/149‑1 (MQ=255)                                                                                                                                                 cGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCa   <  1:280621/149‑1 (MQ=255)                                                                                                                                                  gCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAc  >  1:346394/1‑149 (MQ=255)                                                                                                                                                  gCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGCTCACAGCCGCCTCACCCTGGACGTCGTAGCTGCGCACCGAATCCCAc  >  4:287299/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAC  >  NC_002947/1145301‑1145594

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