Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,448	A→G	10.0%	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	10.0%	159.5 / 7.0	61	intergenic (‑115/+26)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base A (31/23);  new base G (3/3);  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.42e-01

CTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACT  >  NC_002947/1145304‑1145595                                                                                                                                                  |                                                                                                                                                     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                                                                                                            >  8:93717/1‑148 (MQ=255)                                        cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                            >  2:205287/1‑148 (MQ=255)                                        cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                            >  5:279097/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                                                                                             >  6:281821/1‑149 (MQ=255)                                                     cGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAg                                                                                             <  3:320559/149‑1 (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                                                               >  4:239793/1‑149 (MQ=255)                                                                                   gagaCAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                               >  7:225425/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                                                   >  7:55871/1‑149 (MQ=255)                                                                                                aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                                   <  6:151551/149‑1 (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                    >  1:243920/1‑149 (MQ=255)                                                                                                                               tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc                    >  6:289062/1‑149 (MQ=255)                                                                                                                               tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc                    >  6:106543/1‑149 (MQ=255)                                                                                                                               tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc                    >  5:198117/1‑149 (MQ=255)                                                                                                                                          gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa         >  6:330547/1‑149 (MQ=255)                                                                                                                                             cccGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATcc      <  2:243920/149‑1 (MQ=255)                                                                                                                                                 cAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGTGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACt  <  3:159327/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACT  >  NC_002947/1145304‑1145595

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