Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	666,931	Δ1 bp	5.8%	coding (3794/4902 nt)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	666,931	0	A	.	5.8%	254.4 / 10.0	86	coding (3794/4902 nt)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base . (2/3);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GGCAGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGAACGTCGCGC  >  NC_002947/666786‑667078                                                                                                                                                  |                                                                                                                                                     ggCAGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTg                                                                                                                                                  <  2:1257005/150‑1 (MQ=255)   cAGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCa                                                                                                                                                <  2:175321/150‑1 (MQ=255)   cAGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCa                                                                                                                                                <  2:859921/150‑1 (MQ=255)    aGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAg                                                                                                                                               <  1:1143418/150‑1 (MQ=255)    aGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAg                                                                                                                                               <  1:1321759/150‑1 (MQ=255)    aGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTT‑CCTGCAgg                                                                                                                                              <  1:955969/150‑1 (MQ=255)       tttGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑aa                                                                                                                                           >  2:1394337/1‑150 (MQ=255)       tttGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑aa                                                                                                                                           >  2:711174/1‑150 (MQ=255)       tttGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACACAGGGTCAGGCGTGCACGACGCGGCCCGGATTAAGTGCAGTGTGGGGCTGGTTGGTTGCCTTCAGG‑aa                                                                                                                                           <  2:1454979/150‑1 (MQ=255)        ttGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AAc                                                                                                                                          <  2:1547072/150‑1 (MQ=255)          ggCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGc                                                                                                                                        <  2:1693527/150‑1 (MQ=255)           gCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCa                                                                                                                                       >  1:731080/1‑150 (MQ=255)            ccTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTGCCTGCTGCTGGCCCGTGAACCTGGGTACGGCGTGCCAGCCGGGGGACTGGGGATGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCaa                                                                                                                                      <  1:976250/150‑1 (MQ=255)            ccTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCATATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCaa                                                                                                                                      <  1:635076/150‑1 (MQ=255)             cTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAAc                                                                                                                                     <  2:281884/150‑1 (MQ=255)               gtgGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTT‑CCTGCAGGTAACGCagcag                                                                                                                                   <  2:89461/150‑5 (MQ=255)                  gAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGAt                                                                                                                                >  1:247607/1‑150 (MQ=255)                   aGTTCGGACAGAGAGGAAGAGTACTGGGTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATc                                                                                                                               <  1:371630/150‑1 (MQ=255)                    gTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑aa                                                                                                                                           <  2:1738931/137‑1 (MQ=255)                    gTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑aa                                                                                                                                           >  1:1738931/1‑137 (MQ=255)                     ttCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGa                                                                                                                             >  1:485701/1‑150 (MQ=255)                         gACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTc                                                                                                                         <  2:808709/150‑1 (MQ=255)                         gACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTc                                                                                                                         <  2:1769233/150‑1 (MQ=255)                           cAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTAAGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGa                                                                                                                       <  2:731080/150‑1 (MQ=255)                            aGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGAc                                                                                                                      >  2:1579006/1‑150 (MQ=255)                            aGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGAc                                                                                                                      <  2:1264797/150‑1 (MQ=255)                                 ggAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAg                                                                                                                 >  2:1624189/1‑150 (MQ=255)                                    agagTACTGGCTGACCGCCTATGTGAACAACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCga                                                                                                              >  1:1789928/1‑150 (MQ=255)                                           tGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGcc                                                                                                       <  1:1394337/150‑1 (MQ=255)                                            ggCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACa                                                                                                                  >  2:519962/1‑138 (MQ=255)                                            ggCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACa                                                                                                                  <  1:519962/138‑1 (MQ=255)                                             gCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGa                                                                                                     <  2:1902637/150‑1 (MQ=255)                                             gCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGa                                                                                                     >  1:1625630/1‑150 (MQ=255)                                             gCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGa                                                                                                     >  1:861345/1‑150 (MQ=255)                                              cTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGAc                                                                                                    >  1:1878924/1‑150 (MQ=255)                                                    gCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACAcg                                                                                              >  1:1217620/1‑150 (MQ=255)                                                       tATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCTACCACACgcgc                                                                                           >  2:578929/1‑150 (MQ=255)                                                            gactgactTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGc                                                                                      >  2:354611/1‑150 (MQ=255)                                                             actgactTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGAt                                                                                                                                <  2:126022/107‑1 (MQ=255)                                                             actgactTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGAt                                                                                                                                >  1:126022/1‑107 (MQ=255)                                                             actgactTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCg                                                                                     <  2:665100/150‑1 (MQ=255)                                                               tgactTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGt                                                                                   >  1:1015185/1‑150 (MQ=255)                                                                 actTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGc                                                                                 >  2:1523473/1‑150 (MQ=255)                                                                  ctTCCTGCTACGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGACACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCa                                                                                >  1:210914/1‑150 (MQ=255)                                                                   ttCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAg                                                                               <  1:677228/150‑1 (MQ=255)                                                                      ctgctgCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGcc                                                                            >  1:1143409/1‑150 (MQ=255)                                                                           gcGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTAcgc                                                                       <  2:1347016/150‑1 (MQ=255)                                                                                 ccGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCt                                                                 >  2:868649/1‑150 (MQ=255)                                                                                  cGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTg                                                                <  1:1579006/150‑1 (MQ=255)                                                                                  cGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTg                                                                <  1:1259397/150‑1 (MQ=255)                                                                                  cGTGACCAAGGGGACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTCCGCGGCGCTg                                                                <  2:177821/150‑1 (MQ=255)                                                                                    tGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGt                                                              <  1:1496081/150‑1 (MQ=255)                                                                                       ccAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCt                                                           >  1:764997/1‑150 (MQ=255)                                                                                         aaGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGt                                                         >  1:176607/1‑150 (MQ=255)                                                                                         aaGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGt                                                         >  1:1458426/1‑150 (MQ=255)                                                                                             gTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGAc                                                                                                    >  1:1133365/1‑103 (MQ=255)                                                                                             gTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGAc                                                                                                    <  2:1133365/103‑1 (MQ=255)                                                                                                  gCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTc                                                                                                                         <  1:1920297/77‑1 (MQ=255)                                                                                                  gCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTc                                                                                                                         >  2:1920297/1‑77 (MQ=255)                                                                                                   cGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGcagcca                                               <  1:193942/150‑1 (MQ=255)                                                                                                      gCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCcagca                                            <  1:525085/150‑1 (MQ=255)                                                                                                       cccGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCcagcag                                           >  1:1761956/1‑150 (MQ=255)                                                                                                       cccGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCcagcag                                           <  2:1625630/150‑1 (MQ=255)                                                                                                       cccGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGGGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCcagcag                                           >  1:64566/1‑150 (MQ=255)                                                                                                        ccgccgAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAgg                                          >  2:1566644/1‑150 (MQ=255)                                                                                                                 cACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTgg                                 <  2:614209/150‑1 (MQ=255)                                                                                                                  aCTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTggg                                >  1:643299/1‑150 (MQ=255)                                                                                                                  aCTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTggg                                >  1:785700/1‑150 (MQ=255)                                                                                                                   cTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTgg                                 <  1:483558/148‑1 (MQ=255)                                                                                                                   cTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTgg                                 >  2:483558/1‑148 (MQ=255)                                                                                                                    tGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTg                              >  2:390480/1‑150 (MQ=255)                                                                                                                            gCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTgcgc                      >  1:1304470/1‑150 (MQ=255)                                                                                                                            gCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTgcgc                      >  1:1890935/1‑150 (MQ=255)                                                                                                                             ccAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGcgcg                     >  2:721206/1‑150 (MQ=255)                                                                                                                              cAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCgg                    >  1:1806462/1‑150 (MQ=255)                                                                                                                              cAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCgg                    >  1:1389904/1‑150 (MQ=255)                                                                                                                                 tgggcgcCTGCTGCGTTACTTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCt                 >  2:27939/4‑150 (MQ=255)                                                                                                                                 tGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCt                 <  2:1480105/150‑1 (MQ=255)                                                                                                                                     cgcCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTc             <  1:721206/150‑1 (MQ=255)                                                                                                                                        ctgctgCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCTACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGaa          <  1:578929/150‑1 (MQ=255)                                                                                                                                          gttgcGTT‑CCTGCAGGTAACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGAACg        >  1:1769681/3‑150 (MQ=255)                                                                                                                                           ttgcGTT‑CCTGCAGGTAACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGcagcca                                               <  1:1918822/109‑1 (MQ=255)                                                                                                                                           ttgcGTT‑CCTGCAGGTAACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGcagcca                                               >  2:1918822/2‑110 (MQ=255)                                                                                                                                              cGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGAACGTcgc    <  2:210914/150‑1 (MQ=255)                                                                                                                                              cGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTACTTTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGAACGTcgc    >  1:1454862/1‑150 (MQ=255)                                                                                                                                                ttACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCATACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGAACGTCgcgc  >  1:1559590/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GGCAGCTTTGGCCTGTGGAGTTCGGACAGCGAGGAAGAGTACTGGCTGACCGCCTATGTGACTGACTTCCTGCTGCGTGCCCGTGACCAAGGGTACGGCGTGCCCGCCGAGGCACTGAAGAAGGCCAGTGAGCGCCTGCTGCGTTACCTGCAGG‑AACGCAACCTGATCGAAGTCGACTACAGCGAGAATGCCGACCACACGCGCTTTGCGGTGCAGGCCTACGCGGCGCTGGTGCTGTCGCGCAGCCAGCAGGCGCCGTTGGGTGCCCTGCGCGGGCTGTTCGAACGTCGCGC  >  NC_002947/666786‑667078

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