Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,883,081	T→G	L331R (CTG→CGG)	tnaA →	tryptophanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,883,081	0	T	G	92.7%	329.0 / 14.1	97	L331R (CTG→CGG)	tnaA	tryptophanase
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base G (49/40);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-03

CTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCA  >  NZ_CP009273/3882947‑3883213                                                                                                                                       |                                                                                                                                     cttcttTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGAATGAATCTCGACTggcgggct                                                                                                                                  >  1:138270/1‑139 (MQ=255)    cttTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTta                                                                                                                                <  1:750499/138‑1 (MQ=255)    cttTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTta                                                                                                                                <  2:1056667/138‑1 (MQ=255)         aTGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTtatcgtat                                                                                                                          <  2:1138421/139‑1 (MQ=255)         aTGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTtatcgtat                                                                                                                          <  2:514689/139‑1 (MQ=255)           gtgtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTtatcgtatcg                                                                                                                        >  2:880056/1‑139 (MQ=255)            tgtACACCGAGTGCAGAGCCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTTTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgc                                                                                                                       <  1:315897/139‑1 (MQ=255)             gtACACCGAGTGCAGAACCCTTTGCGCCGTTCATGCACGCAACCCGACATATGGCGGCCCGGATCGCCGCGCAATGGAGCGTCTCCCGGTAGGTCTGTATGACCGCATGCATCCCGACCGCCCGGCTTATCGTATcgcg                                                                                                                      <  1:1100970/139‑1 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      <  1:747505/138‑1 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      >  1:874812/1‑138 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      >  2:960938/1‑138 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      >  1:28197/1‑138 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      <  1:907436/138‑1 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      >  1:295441/1‑138 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCgcgc                                                                                                                     <  2:349162/139‑1 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCgcgc                                                                                                                     <  2:622109/139‑1 (MQ=255)              tACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCgcgc                                                                                                                     <  1:1106423/139‑1 (MQ=255)                   cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  1:314172/139‑1 (MQ=255)                   cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  1:610260/139‑1 (MQ=255)                   cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  1:460169/139‑1 (MQ=255)                   cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  1:1127623/139‑1 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCa                                                                                                                    <  1:663174/133‑1 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCa                                                                                                                    >  2:663174/1‑133 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCCggga                                                                                                                >  2:1127623/1‑135 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAgggg                                                                                                                >  2:314172/1‑135 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  2:812478/1‑137 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  2:990857/1‑137 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  2:528078/1‑137 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  1:521879/1‑137 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  1:470564/1‑137 (MQ=255)                     aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  1:1238466/1‑137 (MQ=255)                      gTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAg                                                                                                             >  1:144963/1‑139 (MQ=255)                       tGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCCGGGACAg                                                                                                             >  2:472726/1‑138 (MQ=255)                        gCAGAACCCTTTGGGTGTTGCAGGAAGGCTTTCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGATGGCTTGGTAGTAGACTGGCGGGCTTGTCGTATCGCGCGGGTACAGTa                                                                                                           >  2:1100970/1‑139 (MQ=255)                                 tttGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  >  1:7368/1‑139 (MQ=255)                                 tttGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  >  1:462338/1‑139 (MQ=255)                                   tGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  <  1:4795/137‑1 (MQ=255)                                      gtggtgCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTc                                                                                             <  1:401177/139‑1 (MQ=255)                                         gtgCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTtatcgtat                                                                                                                          >  2:183333/1‑107 (MQ=255)                                         gtgCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTtatcgtat                                                                                                                          <  1:183333/107‑1 (MQ=255)                                          tgCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGa                                                                                         <  1:528078/139‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATcgcg                                                                                                                      <  2:902178/100‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATcgcg                                                                                                                      >  1:902178/1‑100 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCgcgc                                                                                                                     >  1:159583/1‑101 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCgcgc                                                                                                                     <  2:159583/101‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATTGCGCCGGGTCAGGATCTGGGCGGTGGTCTGGGAGAGATTGGc                                                                               >  2:618510/1‑139 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCTCGCAGGTa                                                                                                                >  2:982328/1‑106 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  1:1080148/106‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  2:435192/106‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  2:1080148/1‑106 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                <  1:982328/106‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTa                                                                                                                >  1:435192/1‑106 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  <  1:718167/120‑1 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  >  2:718167/1‑120 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGc                                                                               >  1:729368/1‑139 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGc                                                                               >  1:28686/1‑139 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGGACAGTATCTGGTCGATGGTCTGGAAGAGATTGGc                                                                               >  2:648158/1‑139 (MQ=255)                                                    cTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGGATCGCGCAGGGACAGTATCTGGTCGATGGTCTGGGAGAGATTGGc                                                                               >  2:1102035/1‑139 (MQ=255)                                                        ccGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGtt                                                                            <  1:880056/138‑1 (MQ=255)                                                             atatGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGc                                                                      <  2:934677/139‑1 (MQ=255)                                                               atGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTc                                                                                             >  2:893885/1‑114 (MQ=255)                                                               atGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTc                                                                                             <  1:893885/114‑1 (MQ=255)                                                                      ccTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTTGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcg                                                             >  1:86481/1‑139 (MQ=255)                                                                      ccTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGaa                                                                                        <  1:1216324/112‑1 (MQ=255)                                                                      ccTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGaa                                                                                        >  2:1216324/1‑112 (MQ=255)                                                                      ccTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcg                                                             <  2:144963/139‑1 (MQ=255)                                                                      ccTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGGATCGCGCAGGGACAGTATCTGGTCGATGGTCTGGAAGAGAATGGCGTTGTCTGCCAGCAggcg                                                             >  2:959692/1‑139 (MQ=255)                                                                          gAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcgggcg                                                         <  1:960938/139‑1 (MQ=255)                                                                                ggcgCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCAcg                                                   >  1:691178/1‑139 (MQ=255)                                                                                 gcgcGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCAcgc                                                  >  2:487118/1‑139 (MQ=255)                                                                                  cgcgATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                                                                              >  2:1258635/1‑110 (MQ=255)                                                                                  cgcgATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                                                                              <  1:1258635/110‑1 (MQ=255)                                                                                    cgATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  <  2:170411/88‑1 (MQ=255)                                                                                    cgATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                                                                                  >  1:170411/1‑88 (MQ=255)                                                                                     gATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCa                                              <  1:487118/139‑1 (MQ=255)                                                                                              tCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      >  1:52759/1‑58 (MQ=255)                                                                                              tCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATcgcg                                                                                                                      <  2:52759/58‑1 (MQ=255)                                                                                              tCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGGACAGTATCTGGTCGATGGTCTGGCAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGAt                                     >  2:714931/1‑139 (MQ=255)                                                                                                        ggTCTGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAAc                           <  2:691178/139‑1 (MQ=255)                                                                                                            tGTATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGtt                       >  2:521978/1‑139 (MQ=255)                                                                                                              tATGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTg                      <  2:1116410/138‑1 (MQ=255)                                                                                                                tGACGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCg                   <  2:1020501/139‑1 (MQ=255)                                                                                                                   cGGCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATTCCGGTAAACTGTTGCCGCa                 >  1:794719/1‑138 (MQ=255)                                                                                                                    ggCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATg                                                                                                >  2:430020/1‑58 (MQ=255)                                                                                                                    ggCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATg                                                                                                <  1:430020/58‑1 (MQ=255)                                                                                                                    ggCATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCata               >  1:1133419/1‑139 (MQ=255)                                                                                                                      cATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATc             >  2:471179/1‑139 (MQ=255)                                                                                                                      cATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATc             >  1:172804/1‑139 (MQ=255)                                                                                                                      cATGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGGCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGGGGGCACGCGGCATTCGTTGGTGCCGGGAAAATGTTGCCGCATATc             >  2:287739/1‑139 (MQ=255)                                                                                                                        tGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATccc           >  2:529934/1‑139 (MQ=255)                                                                                                                        tGAATCTCGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATccc           >  1:710951/1‑139 (MQ=255)                                                                                                                              tcGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCg                   <  1:131933/125‑1 (MQ=255)                                                                                                                              tcGACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCg                   >  2:131933/1‑125 (MQ=255)                                                                                                                                gACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCgg                                                <  2:665602/94‑1 (MQ=255)                                                                                                                                gACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCgg                                                >  1:665602/1‑94 (MQ=255)                                                                                                                                gACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGAc    >  2:696041/1‑138 (MQ=255)                                                                                                                                gACTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGAc    <  1:521978/138‑1 (MQ=255)                                                                                                                                 aCTGGCGGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCa  <  2:7368/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     CTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCA  >  NZ_CP009273/3882947‑3883213

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