Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,134,023	C→G	5.0%	A351P (GCA→CCA)	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,134,023	0	C	G	5.0%	252.4 / 4.2	80	A351P (GCA→CCA)	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base G (2/2);  total (40/40)
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.38e-01

CCACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGC  >  NC_002947/4133874‑4134170                                                                                                                                                      |                                                                                                                                                    ccACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcctgc                                                                                                                                                     <  2:335140/150‑1 (MQ=255)   acacCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcctgcct                                                                                                                                                   <  2:671078/150‑1 (MQ=255)   acacCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcctgcct                                                                                                                                                   <  2:1011039/150‑1 (MQ=255)    caccacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGcctgcctg                                                                                                                                                  <  2:753915/150‑1 (MQ=255)     agcggtGGGGCCAGGGATGCTCACGGCGTGGCCTGGTGCCAGCTGCAGAGCGAACTGCGTCAGCACCTGCTGCGCACGACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTgg                                                                                                                                                 <  2:10414/145‑1 (MQ=255)       cacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCg                                                                                                                                               >  1:353919/1‑150 (MQ=255)        acTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACAACGCGCTCATGCAACGCCTGCCTGGCGt                                                                                                                                              <  2:425729/150‑1 (MQ=255)         cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGtt                                                                                                                                             <  1:1134138/150‑1 (MQ=255)            gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCa                                                                                                                                          <  2:606695/150‑1 (MQ=255)              cccAAGGATTCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCCCGCGCGTACCCCAAGCGGATGTCACGGGCCTCCAGCAACTCCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                                        <  2:544652/149‑1 (MQ=255)                   ggATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGgc                                                                                                                                   >  1:794684/1‑150 (MQ=255)                       gCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGgcgcgc                                                                                                                               >  2:1192032/1‑150 (MQ=255)                         tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCaa                                                                                                                             >  1:732181/1‑150 (MQ=255)                             gACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGc                                                                                                                         >  1:975269/1‑150 (MQ=255)                              aCCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCa                                                                                                                        <  2:254724/150‑1 (MQ=255)                                  cACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTAGTGTTGCAGTTCGGCGCGCAATTGCACCAg                                                                                                                    >  1:1071668/1‑150 (MQ=255)                                     cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCt                                                                                                                 >  1:10012/1‑150 (MQ=255)                                     cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCACAATTGCACCAGGCt                                                                                                                 <  2:503491/150‑1 (MQ=255)                                       ggTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAGCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgt                                                                                                               <  2:975269/150‑2 (MQ=255)                                        gTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTggg                                                                                                              >  2:1135395/1‑150 (MQ=255)                                          gCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                             >  2:1414635/1‑149 (MQ=255)                                          gCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                             <  1:1414635/149‑1 (MQ=255)                                             aGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         <  2:353919/150‑1 (MQ=255)                                             aGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         <  2:732181/150‑1 (MQ=255)                                                   aGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCg                                                                                                   <  1:105468/150‑1 (MQ=255)                                                   aGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCg                                                                                                   >  1:622298/1‑150 (MQ=255)                                                     gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTGACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                 >  1:473534/1‑150 (MQ=255)                                                     gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCCCGAGGCGAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                 <  1:374432/150‑1 (MQ=255)                                                     gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                 >  1:928171/1‑150 (MQ=255)                                                     gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                 >  2:93590/1‑150 (MQ=255)                                                         aaCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGa                                                                                             >  2:186188/1‑150 (MQ=255)                                                           cTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAg                                                                                           >  2:1262215/1‑150 (MQ=255)                                                              cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                        <  1:1135395/150‑1 (MQ=255)                                                              cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                        <  1:263812/150‑1 (MQ=255)                                                               gcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGCGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCCCAATTGCACCAGGCTGCGCGTGACGCTCGGCACGAAAGccc                                                                                       >  2:908300/1‑149 (MQ=255)                                                                    cACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTcc                                                                                  >  2:1214469/1‑150 (MQ=255)                                                                       ctgctgCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGc                                                                               >  2:781991/1‑150 (MQ=255)                                                                             ccAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTAGGCGCGCAATTGCACAAGCCCGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcgg                                                                         <  2:49597/150‑1 (MQ=255)                                                                              cAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                        >  2:1292068/1‑150 (MQ=255)                                                                                accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                      <  2:928171/150‑1 (MQ=255)                                                                                accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                      <  2:401121/150‑1 (MQ=255)                                                                                  caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                                    >  1:1237852/1‑150 (MQ=255)                                                                                  caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                                    >  1:1237727/1‑150 (MQ=255)                                                                                     ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGTGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                 <  2:1071668/150‑1 (MQ=255)                                                                                     ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         <  1:1148671/110‑1 (MQ=255)                                                                                     ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         >  2:1148671/1‑110 (MQ=255)                                                                                     ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                 <  2:1344712/150‑1 (MQ=255)                                                                                     ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                 <  1:1192032/150‑1 (MQ=255)                                                                                      cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                <  2:794684/150‑1 (MQ=255)                                                                                        cACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCGATGGCgggagtcggc                                                              <  1:1262215/150‑8 (MQ=255)                                                                                         aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGGAGGCCGACTCCCGCCATCGGCGGGCAAAcccc                                                             >  2:315348/1‑149 (MQ=255)                                                                                           ccgGGTAGCCCAAGCGGATTTCACGGGACTCCAGCAACACCGCGCTCATGCAACGCCTGCTTGGCGTTGCAGTTCTGCGCGCAATTGCCCCAGGCTGGGCGTGACGATCGCCACGAAGGCCGTCTCCCGCCACCGGCGGGCAAACCCTTc                                                           <  2:359037/149‑1 (MQ=255)                                                                                                   cccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                                   >  1:1146168/1‑150 (MQ=255)                                                                                                    ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTc                                                  >  2:742522/1‑150 (MQ=255)                                                                                                        gCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCag                                              >  2:830803/1‑150 (MQ=255)                                                                                                         cGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTCCGGCGCGCAATTGCACCAGGCTGGGCATGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                             >  2:212449/1‑150 (MQ=255)                                                                                                          ggATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                            <  1:186188/150‑1 (MQ=255)                                                                                                          ggATGTCACGGGCATACAGCAACACCACTCTCGACCGGCCGTTGACTTGCGTTGCGTTTCTTCGCGTAGTTGCACCATGCTGTGCGTGACGATCGGCACGAAGGCCGCCTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                            <  2:666347/150‑1 (MQ=255)                                                                                                           gATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCGTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCGCGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAgg                                           >  1:1360509/1‑150 (MQ=255)                                                                                                                     gCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgc                                 <  1:1292068/150‑1 (MQ=255)                                                                                                                        tCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                             >  2:680857/1‑71 (MQ=255)                                                                                                                        tCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                             <  1:680857/71‑1 (MQ=255)                                                                                                                        tCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                              >  1:666112/1‑150 (MQ=255)                                                                                                                              aaCACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcgg                                                                         <  2:853039/101‑1 (MQ=255)                                                                                                                              aaCACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcgg                                                                         >  1:853039/1‑101 (MQ=255)                                                                                                                               acacCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGc                       <  1:830803/150‑1 (MQ=255)                                                                                                                                  ccGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    <  1:212449/150‑1 (MQ=255)                                                                                                                                   cgcgCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc                   >  1:470004/1‑150 (MQ=255)                                                                                                                                        tCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc              >  2:282676/1‑150 (MQ=255)                                                                                                                                          aTGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                             <  2:546255/117‑1 (MQ=255)                                                                                                                                          aTGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                             >  1:546255/1‑117 (MQ=255)                                                                                                                                          aTGCAACGCCTGCCTGCCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGAGCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa            <  1:781991/150‑1 (MQ=255)                                                                                                                                           tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                       >  1:632887/1‑90 (MQ=255)                                                                                                                                           tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                       <  2:632887/90‑1 (MQ=255)                                                                                                                                            gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                     >  1:296898/1‑139 (MQ=255)                                                                                                                                            gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                     <  2:296898/139‑1 (MQ=255)                                                                                                                                            gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt          >  1:592693/1‑150 (MQ=255)                                                                                                                                            gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt          <  2:1237852/150‑1 (MQ=255)                                                                                                                                            gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt          <  2:1237727/150‑1 (MQ=255)                                                                                                                                            gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt          >  2:702401/1‑150 (MQ=255)                                                                                                                                            gCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg           <  2:813947/149‑1 (MQ=255)                                                                                                                                            gCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg           >  1:813947/1‑149 (MQ=255)                                                                                                                                              aaCGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCag                                              >  2:1185073/1‑112 (MQ=255)                                                                                                                                              aaCGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCag                                              <  1:1185073/112‑1 (MQ=255)                                                                                                                                                    tgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACcgc  <  1:908300/150‑1 (MQ=255)                                                                                                                                                    tgcctgGCGTTGCAGTTCGCCGCGATATTCCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGCCCCCCGCCATCGGCGGCCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTGGCCGCCACTGGCCTGCAACTCCAGTTGTACcgc  <  1:315348/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGC  >  NC_002947/4133874‑4134170

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