Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,281,073	T→A	6.1%	intergenic (+88/‑149)	fadH → / → PP_2009	2,4‑dienoyl‑CoA reductase/1‑aminocyclopropane‑1‑carboxylate deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,281,073	0	T	A	6.1%	204.5 / 4.0	66	intergenic (+88/‑149)	fadH/PP_2009	2,4‑dienoyl‑CoA reductase/1‑aminocyclopropane‑1‑carboxylate deaminase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base A (2/2);  total (34/32)
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.19e-01

CGAGCTGGATGCCAAGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCTGTAAACTGCGAT  >  NC_002947/2280928‑2281219                                                                                                                                                  |                                                                                                                                                    cGAGCTGGATGCCAAGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGCACTCCCGACACCTGCACGCGAGGGCGGAGAGATGATCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGc‑c                                                                                                                                                >  1:624570/1‑150 (MQ=255) cGAGCTGGATGCCAAGCGGGCGATCAACCAAGGCGCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGc‑c                                                                                                                                                <  1:401921/150‑1 (MQ=255)      tGGATGCCAAGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGgtgt                                                                                                                                           >  2:209161/1‑150 (MQ=255)       ggATGCCAAGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGtt                                                                                                                                          >  1:721924/1‑150 (MQ=255)            ccAAGCGGGCGATCAACCAAGGCTCGCGGTTGGGGGCTGGGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCgg                                                                                                                                     <  2:68447/150‑1 (MQ=255)              aaGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGt                                                                                                                                   <  1:993504/150‑1 (MQ=255)              aaGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGt                                                                                                                                   <  1:989731/150‑1 (MQ=255)              aaGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGt                                                                                                                                   <  1:848201/150‑1 (MQ=255)                gCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGa                                                                                                                                 >  1:706252/1‑150 (MQ=255)                  gggCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAAc                                                                                                                               >  2:828097/1‑150 (MQ=255)                            ccAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTac                                                                                                                     >  2:1005267/1‑150 (MQ=255)                            ccAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTac                                                                                                                     >  2:662753/1‑150 (MQ=255)                             cAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTacg                                                                                                                    >  1:22865/1‑150 (MQ=255)                             cAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTacg                                                                                                                    >  2:299986/1‑150 (MQ=255)                             cAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGAGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTacg                                                                                                                    >  1:376356/1‑150 (MQ=255)                                   tCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGc                                                                                                              <  1:65660/150‑1 (MQ=255)                                        gTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGc                                                                                                         >  2:912464/1‑150 (MQ=255)                                        gTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGc                                                                                                         <  1:912464/150‑1 (MQ=255)                                           ggcggcTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCgg                                                                                                      <  1:815581/150‑1 (MQ=255)                                            gcggcTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCggg                                                                                                     <  2:721924/150‑1 (MQ=255)                                                 tGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCtgt                                                                                                           >  1:813848/1‑139 (MQ=255)                                                 tGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCtgt                                                                                                           <  2:813848/139‑1 (MQ=255)                                                  gAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCGACTACGGGCTGTGCCGGTGCATCg                                                                                               >  2:140366/1‑150 (MQ=255)                                                      tctGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGc                                                                                           >  2:668276/1‑150 (MQ=255)                                                      tctGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGc                                                                                           <  1:828097/150‑1 (MQ=255)                                                        tGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCt                                                                                         <  2:22865/150‑1 (MQ=255)                                                         gAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGAt                                                                                             >  2:463775/1‑145 (MQ=255)                                                         gAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGAt                                                                                             <  1:463775/145‑1 (MQ=255)                                                          aGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCtgt                                                                                       <  1:662753/150‑1 (MQ=255)                                                           ggTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCtgtg                                                                                      >  2:738964/1‑150 (MQ=255)                                                             ttGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTc                                                                                    >  2:1130214/1‑150 (MQ=255)                                                                tCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAg                                                                                 >  2:950944/1‑150 (MQ=255)                                                                       tGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGTATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCt                                                                          <  1:12396/150‑1 (MQ=255)                                                                         cGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTTGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTgg                                                                        <  1:738964/150‑1 (MQ=255)                                                                          gTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGc                                                                       >  1:1025722/1‑150 (MQ=255)                                                                               cccGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACTCTGGCACtgt                                                                  >  1:725577/1‑150 (MQ=255)                                                                               cccGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACtgt                                                                  <  2:439328/150‑1 (MQ=255)                                                                                ccGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACtgtg                                                                 >  2:1121345/1‑150 (MQ=255)                                                                                 cGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGa                                                                >  2:192034/1‑150 (MQ=255)                                                                                  gACACCTGCACGCGTGGGCGCGCGCCTGCCCGTGAAACAGGTGACGCGGTGGATGGCCCTGTTTTTTCTCGCTG‑TCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGAt                                                               <  1:299986/150‑1 (MQ=255)                                                                                                  ttcgCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCACCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCggcgggggcgg                                               >  2:1061109/4‑150 (MQ=255)                                                                                                  ggAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCggcgggggcgg                                               <  1:154310/150‑1 (MQ=255)                                                                                                  ggAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCAACGGCCTTGA‑CGGTGTTCAAGCGTAATTGCGCCCCTACTACGGTCTTTTCCTGTGCATCTATTCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCggcgggggcgg                                               <  2:970897/150‑1 (MQ=25)                                                                                                   gAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGc                                                             <  1:1024921/135‑1 (MQ=255)                                                                                                   gAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGc                                                             >  2:1024921/1‑135 (MQ=255)                                                                                                   gAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCggcgggggcggg                                              <  1:1005267/150‑1 (MQ=255)                                                                                                     gCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTc                                            >  1:140633/1‑150 (MQ=255)                                                                                                           aTGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGAt                                      <  1:552807/150‑1 (MQ=255)                                                                                                               ccGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACa                                  >  2:1072838/1‑150 (MQ=255)                                                                                                               ccGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACa                                  >  2:497538/1‑150 (MQ=255)                                                                                                               ccGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCAAGCCACCCTGGCACTGTGATGCCGCGGCGTGCGCGGGTCATCGATTACa                                  >  2:57325/1‑150 (MQ=255)                                                                                                                cgcgAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACaa                                 <  1:209161/150‑1 (MQ=255)                                                                                                                 gcgAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAg                                >  1:638570/1‑150 (MQ=255)                                                                                                                  cgAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAgc                               <  1:497538/150‑1 (MQ=255)                                                                                                                  cgAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAgc                               <  1:1072838/150‑1 (MQ=255)                                                                                                                            gACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGa                     <  1:1121345/150‑1 (MQ=255)                                                                                                                                 ggTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTAcc                <  2:1025722/150‑1 (MQ=255)                                                                                                                                   tGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCt              <  1:668276/150‑1 (MQ=255)                                                                                                                                   ccgcgaacACCGGCAATGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCt              <  1:147365/143‑1 (MQ=255)                                                                                                                                     gATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCTGt            >  1:657379/1‑150 (MQ=255)                                                                                                                                        aacACCGGCAATGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGGCGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCTGTaaa         >  2:430904/3‑150 (MQ=255)                                                                                                                                          cACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCggcgggggcgg                                               >  1:91107/1‑110 (MQ=255)                                                                                                                                          cACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCggcgggggcgg                                               <  2:91107/110‑1 (MQ=255)                                                                                                                                          cACCGGCAATGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGa                                                                                  >  1:1106900/1‑75 (MQ=255)                                                                                                                                          cACCGGCAATGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGa                                                                                  <  2:1106900/75‑1 (MQ=255)                                                                                                                                               gCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCTGTAAACTGCGAt  >  2:1037699/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CGAGCTGGATGCCAAGCGGGCGATCAACCAAGGCTCGCGGTTGGCGGCTGAGCTCTGAGGTTGTCGGGGCTGCGTGAGCCCGACACCTGCACGCGTGGGAGCGGGCATGCCCGCGAAACAGGTGACGCGGTGGATGGCACCGGCATTGC‑CGGTGTTCGCGGGTGAACCCGCTCCTACGACGGGCTGTGCCGGGGCATCGATGCCTGTGTCGAGCCACCCTGGCACTGTGATGCCTCGGCGGGGGCGGGTCATCGATTACAAGCGCTTGCGGGATTACCCTGTAAACTGCGAT  >  NC_002947/2280928‑2281219

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