Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,440	G→A	8.0%	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	8.0%	171.1 / 2.9	61	intergenic (‑107/+34)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (3/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAA  >  NC_002947/1145295‑1145588                                                                                                                                                   |                                                                                                                                                      gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                                    >  8:34684/1‑149 (MQ=255)  ttCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGcc                                                                                                                                                   <  2:104734/149‑1 (MQ=255)     agtcgcTGGGAGAGCAGGGCTGGGCGGTCTGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGc                                                                                                                                                <  7:152613/146‑1 (MQ=255)          cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                           <  4:315160/149‑1 (MQ=255)          cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                           <  1:53602/149‑1 (MQ=255)             ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAagag                                                                                                                                          <  3:231434/147‑1 (MQ=255)             ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                        >  2:272039/1‑149 (MQ=255)                  gcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCGGGTa                                                                                                                                   >  8:104453/1‑149 (MQ=255)                     gggCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTActg                                                                                                                                <  1:216376/149‑3 (MQ=255)                        cTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACg                                                                                                                             <  5:65622/149‑1 (MQ=255)                                tcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGcc                                                                                                                      >  6:307847/1‑148 (MQ=255)                                    tcATGGGTGAGTCCTTGGAAAGAGAGCGTATAGAAACAAGACGAAACATCGGGAGCTAGACAGTAAAA‑AAAGCCGCTATTTGCGTTGCCAGTACCGGCCTCTTCGCCGCCGAGCCCGCAAAGAGGCAGATACACACGTAACgaatcac                                                                                                                  >  4:296635/1‑142 (MQ=255)                                         ggggAGTCCTTGGAAAGATGGCGTCTAGATACAAGACGACACAGCGGGAGCGAGACAGTTAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATcc                                                                                                            <  4:46218/149‑1 (MQ=255)                                               tCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAAcc                                                                                                      <  6:170460/149‑1 (MQ=255)                                                 cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCTAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                     <  3:304703/148‑1 (MQ=255)                                                 cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                     >  5:279097/1‑148 (MQ=255)                                                 cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                     >  8:93717/1‑148 (MQ=255)                                                 cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                     >  2:205287/1‑148 (MQ=255)                                                     gAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                                 >  6:121304/1‑148 (MQ=255)                                                      aaaGAGGGCGTCTAGATACAAGACGAACCAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                               >  1:254903/1‑149 (MQ=255)                                                         gagGGCGTCTAGATACAAGACGAAACAGCGGTAACGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TAc                                                                                           <  1:112448/149‑1 (MQ=255)                                                             gCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCg                                                                                         >  4:320559/1‑147 (MQ=255)                                                              cGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAg                                                                                      <  3:320559/149‑1 (MQ=255)                                                              cGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAg                                                                                      >  6:281821/1‑149 (MQ=255)                                                                    gATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTgcgc                                                                                <  5:112337/149‑1 (MQ=255)                                                                       aCAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGcccc                                                                             >  4:105027/1‑149 (MQ=255)                                                                        cAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCg                                                                            >  4:201394/1‑149 (MQ=255)                                                                           gACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCa                                                                         >  8:295779/1‑149 (MQ=255)                                                                              gAAACAGCGGGAGCGAGACAGTAAAAGAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATg                                                                       >  1:3665/1‑149 (MQ=255)                                                                              gAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATgg                                                                      >  4:190381/1‑149 (MQ=255)                                                                              gAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATgg                                                                      >  8:164279/1‑149 (MQ=255)                                                                                  cAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAg                                                                  <  3:201394/149‑1 (MQ=255)                                                                                        gAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc                                                            <  5:307847/149‑1 (MQ=255)                                                                                            gagaCAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                        >  7:225425/1‑149 (MQ=255)                                                                                            gagaCAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                        >  4:239793/1‑149 (MQ=255)                                                                                            gagaCAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCACAAACCACATTGA‑TACGGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                        >  6:172935/1‑149 (MQ=255)                                                                                              gaCAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                                      >  4:271487/1‑149 (MQ=255)                                                                                                cAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtgg                                                    <  5:121304/149‑1 (MQ=255)                                                                                                cAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                     <  6:279097/148‑1 (MQ=255)                                                                                                cAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                     <  2:254903/148‑1 (MQ=255)                                                                                                 aGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAAGGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCGGGCATGGTCAGGCGGTCGATGTtggc                                                   <  3:296635/149‑1 (MQ=255)                                                                                                  gTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACAGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggct                                                  >  5:161404/1‑149 (MQ=255)                                                                                                   tAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctg                                                 >  5:151551/1‑149 (MQ=255)                                                                                                    aaaa‑aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctgg                                                >  2:60716/1‑149 (MQ=255)                                                                                                         aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                            <  6:151551/149‑1 (MQ=255)                                                                                                         aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                            >  7:55871/1‑149 (MQ=255)                                                                                                            gCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGgctggctgtgac                                         <  7:164279/149‑5 (MQ=255)                                                                                                              cGGTATTTGCGTTGCCAGTACCGGCCTGTTCGCGGGCGGGCCCGCTAAGAGGCCGGTACAGACGCAAGGCCGCAGATATCCGATTCCACATTGA‑TACCGACGTTGCGCCCCGGCATGGTCAGGGGGTCGATGTTGGCTGGCTGGGTcac                                       <  4:200434/149‑1 (MQ=255)                                                                                                               ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGTGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                                      >  4:159327/1‑149 (MQ=255)                                                                                                                   tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGccc                                                                              >  3:217742/1‑105 (MQ=255)                                                                                                                   tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGccc                                                                              <  4:217742/105‑1 (MQ=255)                                                                                                                   tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgccg                                  <  1:60716/149‑1 (MQ=255)                                                                                                                   tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgcc                                   <  7:295779/148‑1 (MQ=255)                                                                                                                       cGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGACTACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGTGACCCAgccagcc                               <  1:272039/149‑5 (MQ=255)                                                                                                                         ttGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAc                            <  8:225425/149‑1 (MQ=255)                                                                                                                            ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCt                         >  2:47896/1‑149 (MQ=255)                                                                                                                                        tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  1:312917/1‑149 (MQ=255)                                                                                                                                        tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  6:289062/1‑149 (MQ=255)                                                                                                                                        tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  1:243920/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  5:198117/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  6:106543/1‑149 (MQ=255)                                                                                                                                                   gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa  >  6:330547/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAA  >  NC_002947/1145295‑1145588

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