Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,335,001	Δ1 bp	5.9%	coding (604/1704 nt)	PP_2932 ←	amidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,335,001	0	A	.	5.9%	219.5 / 9.9	85	coding (604/1704 nt)	PP_2932	amidase family protein
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base . (3/2);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

GCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAG  >  NC_002947/3334854‑3335149                                                                                                                                                     |                                                                                                                                                     gcCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACg                                                                                                                                                    <  2:1870005/150‑1 (MQ=255)  ccAAAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGGTTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                                   <  1:379733/150‑1 (MQ=255)   cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CTCGACggg                                                                                                                                                  <  1:1523935/150‑1 (MQ=255)   cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGCCggg                                                                                                                                                  >  2:804331/1‑150 (MQ=255)   cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                                  <  1:2095008/150‑1 (MQ=255)   cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                                  <  1:1523852/150‑1 (MQ=255)   cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTAACACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                                  <  1:1946536/150‑1 (MQ=255)       ggTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGAACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgtgt                                                                                                                                              >  1:1071502/1‑150 (MQ=255)          cgccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAgg                                                                                                                                           <  2:14779/150‑1 (MQ=255)             cgcgCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgcgc                                                                                                                                        >  1:750538/1‑150 (MQ=255)              gcgcTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGGGTAGGCGca                                                                                                                                       >  1:327041/1‑150 (MQ=255)                 cTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGAGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAg                                                                                                                                    <  2:1998215/150‑1 (MQ=255)                          gTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGATATCACCCCG‑CGCGTCGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                           <  2:1465791/150‑1 (MQ=255)                           tCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCTGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                          <  2:396628/150‑1 (MQ=255)                           tCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                          <  1:600846/150‑1 (MQ=255)                            cATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                                   <  2:2323358/140‑1 (MQ=255)                            cATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                                   >  1:2323358/1‑140 (MQ=255)                             aTCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGa                                                                                                                        <  1:1619811/150‑1 (MQ=255)                              tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGcc                                                                                                                                 >  1:1941887/1‑140 (MQ=255)                              tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGcc                                                                                                                                 <  2:1941887/140‑1 (MQ=255)                              tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGAc                                                                                                                       >  1:2131898/1‑150 (MQ=255)                              tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGAc                                                                                                                       >  2:2147222/1‑150 (MQ=255)                                 gCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgggg                                                                                                                                                 >  1:470154/1‑121 (MQ=255)                                 gCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgggg                                                                                                                                                 <  2:470154/121‑1 (MQ=255)                                    accaccACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACAACCACGTCCATGGTCGACAAGGACGGCCAGTTGCCAAGCAAAGAAATAATTGCG‑TGCGAGGGGGTGTAGGCGCAGAGGCCGTTGTTGGACGCCggg                                                                                                                 >  1:1117171/1‑150 (MQ=255)                                     ccaccaCGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGc                                                                                                                <  1:433928/150‑1 (MQ=255)                                      caccacGTCGAGGATTTCAAGGAGGTCGGCTATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCCGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGcc                                                                                                               >  2:644723/1‑150 (MQ=255)                                      caccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGcc                                                                                                               >  2:1420951/1‑150 (MQ=255)                                              cGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGc                                                                                                       <  2:766268/150‑1 (MQ=255)                                               gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACg                                                                                                            >  2:774518/1‑144 (MQ=255)                                               gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACg                                                                                                            <  1:774518/144‑1 (MQ=255)                                                 ggATTGCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTgg                                                                                                    >  1:1365258/1‑150 (MQ=255)                                                  gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAg                                                                                                                                            <  2:1523491/109‑1 (MQ=255)                                                  gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAg                                                                                                                                            >  1:1523491/1‑109 (MQ=255)                                                  gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGAGAAACGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGa                                                                                                   <  1:883508/150‑1 (MQ=255)                                                    tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgatttc                                                                                                                                          <  1:679675/109‑6 (MQ=255)                                                    tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgatttc                                                                                                                                          >  2:679675/1‑104 (MQ=255)                                                    tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAcc                                                                                                 <  1:947405/150‑1 (MQ=255)                                                      tCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                               <  2:1538518/150‑1 (MQ=255)                                                       cAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTAGGCGTTAGCGGCCCGTTGCCACGCAGCGCAATCACACCG‑CGCGACGGGGTGGGGGCTCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAgg                                                                                              >  2:2101074/1‑150 (MQ=255)                                                        aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCA‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                                     <  1:108083/142‑1 (MQ=255)                                                        aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCA‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                                     >  2:108083/1‑142 (MQ=255)                                                               tCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                                      <  1:476552/150‑1 (MQ=255)                                                                cGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCgg                                                                                     <  2:1433169/150‑1 (MQ=255)                                                                 ggCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                             >  1:1153840/1‑141 (MQ=255)                                                                 ggCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                             <  2:1153840/141‑1 (MQ=255)                                                                    cATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccag                                                                                 <  1:1420951/150‑1 (MQ=255)                                                                     aTGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccagg                                                                                <  1:419442/150‑1 (MQ=255)                                                                              gCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAg                                                                                                                                    >  2:513123/1‑89 (MQ=255)                                                                              gCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAg                                                                                                                                    <  1:513123/89‑1 (MQ=255)                                                                                aTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCt                                                                     >  1:2186430/1‑150 (MQ=255)                                                                                   gggCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGaa                                                                  <  1:2147222/150‑1 (MQ=255)                                                                                     gaaccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGc                                                                <  2:1071502/148‑1 (MQ=255)                                                                                     gCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGc                                                                >  2:545688/1‑150 (MQ=255)                                                                                        ccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTgg                                                             >  1:2390159/1‑150 (MQ=255)                                                                                         caccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAACAGGTCTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggc                                                            >  2:302533/1‑150 (MQ=255)                                                                                           ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggcgg                                                          <  1:1825316/150‑1 (MQ=255)                                                                                                  cATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggc                                                   >  1:1748400/1‑150 (MQ=255)                                                                                                   aTGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcg                                                  <  2:613321/150‑1 (MQ=255)                                                                                                     ggTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcggt                                                <  1:230605/150‑1 (MQ=255)                                                                                                       tCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTAc                                              <  1:2361676/150‑1 (MQ=255)                                                                                                         ggCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTCTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCg                                            >  2:171646/1‑150 (MQ=255)                                                                                                                  ggCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTg                                   <  2:1748400/150‑1 (MQ=255)                                                                                                                   gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                                  >  2:922472/1‑150 (MQ=255)                                                                                                                   gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                                  <  2:2186430/150‑1 (MQ=255)                                                                                                                   gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                                  >  1:1358625/1‑150 (MQ=255)                                                                                                                     cAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGAc                                <  1:20051/150‑1 (MQ=255)                                                                                                                       gTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGa                              >  2:1743698/1‑150 (MQ=255)                                                                                                                            cACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCAGGAACCGGGTTCTTCGGCCAGGCCGGCGGCGAAAAAGCCGGCATTGCTTTCTGTACCTGCACCGTTGGACGAGCCCg                         >  1:532393/1‑150 (MQ=255)                                                                                                                             aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGa                        >  2:486390/1‑150 (MQ=255)                                                                                                                             aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGCCCGTTGTTGGACGCCGGGCCCCGGCCGCTGGACCCGGTCTCTTCGGCCCGGATCTCGGCCCTCGACCAGTCTGCGGGGGCCGTCCCGGCACCGTTGGACGAGCCCGa                        >  2:2367419/1‑150 (MQ=255)                                                                                                                                   cGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGaa                  >  2:423849/1‑150 (MQ=255)                                                                                                                                   cGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGaa                  >  1:915425/1‑150 (MQ=255)                                                                                                                                    gAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGANGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAg                 <  1:495849/150‑1 (MQ=255)                                                                                                                                    gAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAg                 >  1:1089192/1‑150 (MQ=255)                                                                                                                                     aaaTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAgg                <  2:1358625/150‑1 (MQ=255)                                                                                                                                      ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                               <  2:204091/67‑1 (MQ=37)                                                                                                                                      ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                               >  1:204091/4‑70 (MQ=37)                                                                                                                                       tacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCTCGGCCACAGCCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgc              >  1:1123668/3‑150 (MQ=255)                                                                                                                                         cACCCCG‑CGCGACGGGGTGGAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgcgc            <  2:1119646/150‑1 (MQ=255)                                                                                                                                            cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                                   <  2:1229788/92‑1 (MQ=255)                                                                                                                                            cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                                   >  1:1229788/1‑92 (MQ=255)                                                                                                                                            cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt         >  1:1168993/1‑150 (MQ=255)                                                                                                                                            cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt         >  2:1637018/1‑150 (MQ=255)                                                                                                                                                    gACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGATTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGACAGGTAg  >  1:701327/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAG  >  NC_002947/3334854‑3335149

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