Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,433	C→T	10.9%	intergenic (‑100/+41)	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,433	0	C	T	10.9%	170.4 / 9.1	64	intergenic (‑100/+41)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (4/3);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

GATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGC  >  NC_002947/1145289‑1145577                                                                                                                                                  |                                                                                                                                                  gATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCgg                                                                                                                                               >  2:216376/1‑149 (MQ=255)  aTCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCgg                                                                                                                                               <  5:145730/148‑1 (MQ=255)       gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGcac                                                                                                                                           <  7:240202/147‑3 (MQ=255)       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                                                                                                                                <  1:53602/149‑1 (MQ=255)                cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                <  4:315160/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                                                                           >  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                                                           >  8:164279/1‑149 (MQ=255)                                                                                    gAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATgg                                                           >  4:190381/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                                 <  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  >  6:289062/1‑149 (MQ=255)                                                                                                                                              tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc  >  1:312917/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)                                                                                                                                                  |                                                                                                                                                  GATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAA‑AAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGA‑TACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGC  >  NC_002947/1145289‑1145577

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