Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	664,889	C→T	6.3%	A584A (GCC→GCT)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	664,889	0	C	T	6.3%	261.7 / 4.8	90	A584A (GCC→GCT)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base T (3/3);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCTACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGA  >  NC_002947/664741‑665038                                                                                                                                                     |                                                                                                                                                      gCTACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCa                                                                                                                                                      <  2:1224818/150‑1 (MQ=255)     ccAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCt                                                                                                                                                  >  1:618076/1‑150 (MQ=255)      cAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTg                                                                                                                                                 >  1:76815/1‑150 (MQ=255)       aaTTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGc                                                                                                                                                >  1:1171303/1‑150 (MQ=255)        aTTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCa                                                                                                                                               >  2:1116725/1‑150 (MQ=255)          tCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAgg                                                                                                                                             >  2:444504/1‑150 (MQ=255)                      ccGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggcc                                                                                                                                 <  2:76815/150‑1 (MQ=255)                        gAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCCTGCAGGCTGGCCTGCacggtc                                                                                                                                          <  2:930183/139‑6 (MQ=255)                        gAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCCTGCAGGCTGGCCTGCacggtc                                                                                                                                          >  1:930183/1‑134 (MQ=255)                                 cTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCAcc                                                                                                                      <  2:1668899/150‑1 (MQ=255)                                  tGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCAccc                                                                                                                     >  2:1586116/1‑150 (MQ=255)                                    aaCCAGGACCTGGAACTGGACGAAGTCACCCTTGACCCGGCCGGCAAGGCGGTGGGCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACGGACCGTGCAGGCCAGCCTGCAGGAGTCGGGGGGCCGGCCGATCACCCgg                                                                                                                   >  1:603934/1‑150 (MQ=255)                                      ccAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGgcg                                                                                                                 >  1:1463137/1‑150 (MQ=255)                                             cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                          >  1:1191108/1‑150 (MQ=255)                                             cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                          >  2:1233637/1‑150 (MQ=255)                                              tGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagc                                                                                                         >  2:749577/1‑150 (MQ=255)                                                gAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcag                                                                                                       <  2:685174/150‑1 (MQ=255)                                                 aaCTTGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcagc                                                                                                      <  1:1023455/150‑1 (MQ=255)                                                 aaCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcagc                                                                                                      <  2:1730099/150‑1 (MQ=255)                                                 aaCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTCCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGTCCGATCACCCGGCGCCTGGagcagc                                                                                                      <  1:1023882/150‑1 (MQ=255)                                                     ggACGAAGTCACCCTTGACCAGGTCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGATCGGCCGGCCACCCgact                                                                                                                 <  1:195177/135‑4 (MQ=255)                                                     ggACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTAt                                                                                                  <  2:1843346/150‑1 (MQ=255)                                                     ggACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGATCGGCCGGCCACCCgact                                                                                                                 >  2:195177/1‑132 (MQ=255)                                                      gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                                 >  2:372583/1‑150 (MQ=255)                                                                 ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGc                                                                                                                                                    <  1:1883076/88‑1 (MQ=255)                                                                 ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGc                                                                                                                                                    >  2:1883076/1‑88 (MQ=255)                                                                  cTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAg                                                                                     <  2:1725459/150‑1 (MQ=255)                                                                         aGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGc                                                                              >  1:342432/1‑150 (MQ=255)                                                                          ggCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGGCGGGTGGCCGGCCTATCACCCGGCGCCTGGCGCAGCCTATCTGGCCGGCTGAGCGCCTGcc                                                                             >  1:1888659/1‑150 (MQ=255)                                                                            ccGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  1:372583/150‑1 (MQ=255)                                                                                 aaGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTg                                                                      >  1:637389/1‑150 (MQ=255)                                                                                  aGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGc                                                                     <  1:1607498/150‑1 (MQ=255)                                                                                  aGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGc                                                                     <  1:1502263/150‑1 (MQ=255)                                                                                    gCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGa                                                                   <  2:618076/150‑1 (MQ=255)                                                                                       gtggtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                <  2:784201/150‑1 (MQ=255)                                                                                         ggtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCt                                                              <  1:749577/150‑1 (MQ=255)                                                                                          gtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTTAGCGCCTGCCGGGCCTGCGAGGCCTg                                                             >  1:360272/1‑150 (MQ=255)                                                                                          gtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           <  2:1073837/136‑1 (MQ=255)                                                                                          gtCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                           >  1:1073837/1‑136 (MQ=255)                                                                                           tCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGt                                                            <  1:1233637/150‑1 (MQ=255)                                                                                            cGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCGGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGtt                                                           <  1:787396/150‑1 (MQ=255)                                                                                                   gAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGt                                                    >  1:1596389/1‑150 (MQ=255)                                                                                                     aaGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                                  <  1:444504/150‑1 (MQ=255)                                                                                                      aGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                <  1:887529/135‑1 (MQ=255)                                                                                                      aGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                >  2:887529/1‑135 (MQ=255)                                                                                                         cGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgaggag                                              >  1:1875151/1‑150 (MQ=255)                                                                                                               gCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGTGAGGCCTGTTCGAAGGTGAGGAGACCGAt                                        >  1:1548347/1‑150 (MQ=255)                                                                                                               gCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGAt                                        >  1:586152/1‑150 (MQ=255)                                                                                                                ccGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                       >  2:634921/1‑150 (MQ=255)                                                                                                                 cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAgagct                                                              >  1:972762/1‑122 (MQ=255)                                                                                                                 cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                >  2:196608/1‑124 (MQ=255)                                                                                                                 cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                <  1:196608/124‑1 (MQ=255)                                                                                                                 cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCt                                                              <  2:972762/126‑1 (MQ=255)                                                                                                                 cGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAg                                      >  1:1899675/1‑150 (MQ=255)                                                                                                                  gAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGc                                     >  2:362849/1‑150 (MQ=255)                                                                                                                  gAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGAAGCCTATCTGGCCGGCTGAGAGCCTGCCGGTCCTGCTAGGCCTGTTCGAAGGTGAGGAGACCGATAGc                                     <  2:1191108/150‑1 (MQ=255)                                                                                                                    aGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGACGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGa                                   >  2:1440907/1‑150 (MQ=255)                                                                                                                      cccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACg                                 >  2:1055823/1‑150 (MQ=255)                                                                                                                        cGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                                  >  1:1896211/1‑131 (MQ=255)                                                                                                                        cGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                                  <  2:1896211/131‑1 (MQ=255)                                                                                                                         gCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGACGGTGAGGAGACCGATAGCGACGGGc                              >  2:621953/1‑150 (MQ=255)                                                                                                                           tctcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGccc                            >  2:1917508/1‑150 (MQ=255)                                                                                                                             tcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                                  <  1:955511/126‑1 (MQ=255)                                                                                                                             tcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                                  >  2:955511/1‑126 (MQ=255)                                                                                                                              ccGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTg                         <  2:1171303/150‑1 (MQ=255)                                                                                                                              ccGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTg                         <  1:1586116/150‑1 (MQ=255)                                                                                                                               cGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACgg                                >  1:635030/1‑142 (MQ=255)                                                                                                                               cGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACgg                                <  2:635030/142‑1 (MQ=255)                                                                                                                                 cTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGaa                      <  2:1596389/150‑1 (MQ=255)                                                                                                                                  tGCAACTGACCGTGCAGGCCAGCCTGCAGGAGGCGGGTGTCCGGCCGATCACCCGGCTCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAAt                     >  2:423701/1‑150 (MQ=255)                                                                                                                                   gCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCGGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCGGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAAtt                    <  1:1116725/150‑1 (MQ=255)                                                                                                                                    cAACTGACCGTGCAGGCCCGCCTGCAGGAGGAGGGTGGCCGGCCGATCACCGGGCGCCTGGAGGAGCCGATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCTACGGGCCCGTGGAATTc                   <  2:1888659/150‑1 (MQ=255)                                                                                                                                     aaCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCg                  <  2:1434385/150‑1 (MQ=255)                                                                                                                                     aaCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCg                  >  2:381976/1‑150 (MQ=255)                                                                                                                                      aCTGCCCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  2:1527694/150‑1 (MQ=255)                                                                                                                                      aCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  2:748073/150‑1 (MQ=255)                                                                                                                                      aCTGACCGGGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  1:1165488/150‑1 (MQ=255)                                                                                                                                       cTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAg                >  1:257519/1‑150 (MQ=255)                                                                                                                                        tGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGt               >  2:1022351/1‑150 (MQ=255)                                                                                                                                        tGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGt               >  2:1524889/1‑150 (MQ=255)                                                                                                                                          aCCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTc             <  2:586152/150‑1 (MQ=255)                                                                                                                                             ctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGcc                                                                                            >  1:876575/2‑68 (MQ=255)                                                                                                                                             ctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGcc                                                                                            <  2:876575/67‑1 (MQ=255)                                                                                                                                              tGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 >  1:969555/1‑142 (MQ=255)                                                                                                                                              tGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa                 <  2:969555/142‑1 (MQ=255)                                                                                                                                              tGCAGGCCAGCCTGGAGTAGTCGGGTGGCCTGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtgg         <  2:603934/150‑1 (MQ=255)                                                                                                                                              tGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtgg         >  1:1664951/1‑150 (MQ=255)                                                                                                                                              tGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCATGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtgg         >  1:33854/1‑150 (MQ=255)                                                                                                                                                cAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggtg       >  2:1697032/1‑150 (MQ=255)                                                                                                                                                   gCCAGCCTGCAGGAGTCGGGTTGCCGGCCGGTCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGCGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGTcc    >  2:348523/1‑150 (MQ=255)                                                                                                                                                     cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGa  >  1:652419/1‑150 (MQ=255)                                                                                                                                                     cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGa  >  2:202381/1‑150 (MQ=255)                                                                                                                                                     cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTGCCTGGTGCCCGa  >  2:517468/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCTACCAATTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGA  >  NC_002947/664741‑665038

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