Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	667,326	C→G	5.3%	R1397G (CGT→GGT)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	667,326	0	C	G	5.3%	240.9 / 4.4	76	R1397G (CGT→GGT)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (32/40);  new base G (2/2);  total (34/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 = 7.20e-01

ATTAGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGAAGCTTGAGGGC  >  NC_002947/667179‑667475                                                                                                                                                    |                                                                                                                                                       aTTAGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGt                                                                                                                                                      <  2:84278/150‑1 (MQ=255)  ttAGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGGTCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGtt                                                                                                                                                     <  1:1388890/150‑1 (MQ=255)  ttAGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGtt                                                                                                                                                     >  2:1351233/1‑150 (MQ=255)   tAGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTg                                                                                                                                                    <  1:46011/150‑1 (MQ=255)    aGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTgg                                                                                                                                                   <  2:521688/150‑1 (MQ=255)     gCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGc                                                                                                                                                  >  1:1544708/1‑150 (MQ=255)      ccGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCt                                                                                                                                                 <  2:1299155/150‑1 (MQ=255)       cGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTg                                                                                                                                                >  1:510407/1‑150 (MQ=255)        gTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGt                                                                                                                                               <  1:746129/150‑1 (MQ=255)                ggctggctGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCaa                                                                                                                                       <  1:627414/150‑1 (MQ=255)                  ctggctggCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAaga                                                                                                                                     <  1:1398122/150‑1 (MQ=255)                  ctggctggCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGCTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAaga                                                                                                                                     <  1:791037/150‑1 (MQ=255)                         gcCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACGGTTGGCTgc                                                                                                                                               >  2:1476359/1‑132 (MQ=255)                         gcCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACGGTTGGCTgc                                                                                                                                               <  1:1476359/133‑2 (MQ=255)                             aTTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTAGAAGAGAGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCAc                                                                                                                          >  1:328448/1‑150 (MQ=255)                              ttATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACt                                                                                                                         <  2:135196/150‑1 (MQ=255)                                  ggCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTc                                                                                                                     <  1:749933/150‑1 (MQ=255)                                    cAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCt                                                                                                                   <  2:1299943/150‑1 (MQ=255)                                     aGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACGGTTGGCTgc                                                                                                                                               >  1:1449538/1‑120 (MQ=255)                                     aGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACGGTTGGCTgc                                                                                                                                               <  2:1449538/121‑2 (MQ=255)                                      gcgcCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCg                                                                                                                 >  1:387073/1‑150 (MQ=255)                                        gcCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑Ct                                                                                                              >  1:1463784/1‑150 (MQ=255)                                         cccTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTg                                                                                                             >  1:1411560/1‑150 (MQ=255)                                         cccTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTg                                                                                                             >  1:763435/1‑150 (MQ=255)                                          ccTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGGACTGGCAGCCAACCGTTGGCTGTCGACCCCAACGCGCAATGCACTGTTCCTCG‑CTgg                                                                                                            <  1:521868/150‑1 (MQ=255)                                                tGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTgg                                                                                                      <  2:1769076/150‑1 (MQ=255)                                                  aCCAGGCACTGATCCTGGCGCTGCTGCAGGACCGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGGGCGCAAGGCATTGGTCCTCG‑GGGGCCGTGGcc                                                                                                    <  2:1454862/150‑1 (MQ=255)                                                       gCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCctgctg                                                                                               >  2:1386120/1‑150 (MQ=255)                                                            gATCCTGGCGCTGCGGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGTCCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCaa                                                                                          <  2:458239/150‑1 (MQ=255)                                                              tCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAaga                                                                                                                                     <  2:183849/106‑1 (MQ=255)                                                              tCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAaga                                                                                                                                     >  1:183849/1‑106 (MQ=255)                                                               ccTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGcc                                                                                       <  1:691197/150‑1 (MQ=255)                                                                cTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCg                                                                                      <  1:272638/150‑1 (MQ=255)                                                                cTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCg                                                                                      >  2:130431/1‑150 (MQ=255)                                                                 tGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCgg                                                                                     >  1:1599064/1‑150 (MQ=255)                                                                  ggCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGa                                                                                    <  1:1192749/150‑1 (MQ=255)                                                                   gcgcTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGCGCAAGCCGGaa                                                                                   <  2:556535/150‑1 (MQ=255)                                                                    cgcTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAg                                                                                  <  2:981993/150‑1 (MQ=255)                                                                     gctgctGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTc                                                                                                                                              <  2:171939/90‑1 (MQ=255)                                                                     gctgctGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTc                                                                                                                                              >  1:171939/1‑90 (MQ=255)                                                                      ctgctgCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGc                                                                                <  1:130431/150‑1 (MQ=255)                                                                         ctgcAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAg                                                                             <  1:830245/150‑1 (MQ=255)                                                                             aGGAAAGCAACCAGCCAAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTggc                                                                         >  1:1005139/1‑150 (MQ=255)                                                                               gAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTg                                                                                                       >  1:1355928/1‑118 (MQ=255)                                                                               gAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTg                                                                                                       <  2:1355928/118‑1 (MQ=255)                                                                               gAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTggcag                                                                       <  1:271810/150‑1 (MQ=255)                                                                                  aGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGAc                                                                                                                                           <  1:1341033/80‑1 (MQ=255)                                                                                  aGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGAc                                                                                                                                           >  2:1341033/1‑80 (MQ=255)                                                                                  aGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCATGGCCTGCTGGGCAAGCCGGAAGGCCAGTggcaggca                                                                    >  1:1454603/1‑150 (MQ=255)                                                                                     aaCCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGc                                                                 >  2:968971/1‑150 (MQ=255)                                                                                        cTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTg                                                              >  1:1639180/1‑150 (MQ=255)                                                                                             cagcagCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAg                                                         >  1:29829/1‑150 (MQ=255)                                                                                             cagcagCCAGGTCGATCAGCACCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAg                                                         >  1:963665/1‑150 (MQ=255)                                                                                               gcagcCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAgcg                                                       >  2:1630289/1‑150 (MQ=255)                                                                                                  gcCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCgg                                                    <  2:246357/150‑1 (MQ=255)                                                                                                   ccAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCggg                                                   <  2:1559567/150‑1 (MQ=255)                                                                                                       gTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATg                                                                                                                             <  1:1911971/73‑1 (MQ=255)                                                                                                       gTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATg                                                                                                                             >  2:1911971/1‑73 (MQ=255)                                                                                                       gTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGaa                                               <  1:1528245/150‑1 (MQ=255)                                                                                                        tCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAg                                              <  2:775570/150‑1 (MQ=255)                                                                                                            tCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGACAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGGAAGCCGGAAGGCCAGTGGCCGGCACGCCTGGACAGCGCGGGTGAAGTgcg                                          >  2:1017833/1‑150 (MQ=255)                                                                                                               gcgcCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGaa                                               >  1:268872/1‑142 (MQ=255)                                                                                                               gcgcCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGaa                                               <  2:268872/142‑1 (MQ=255)                                                                                                               gcgcCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGa                                       >  1:963027/1‑150 (MQ=255)                                                                                                                    tGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCATGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGTACAGCGCGGGTGAAGTGCGCGAGTTGa                                  <  2:1454603/150‑1 (MQ=255)                                                                                                                     gTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGaa                                 <  2:963027/150‑1 (MQ=255)                                                                                                                            ttGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGcc                          >  2:1598426/1‑150 (MQ=255)                                                                                                                                  gATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGcc                    >  1:335036/1‑150 (MQ=255)                                                                                                                                   aTGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCg                   <  2:841931/150‑1 (MQ=255)                                                                                                                                   aTGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCg                   <  1:1444180/150‑1 (MQ=255)                                                                                                                                         tGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTAGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGa             >  1:706858/1‑150 (MQ=255)                                                                                                                                         tGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGGTCCTCGAAGGGGCGT‑TCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGa             >  2:1005930/1‑150 (MQ=255)                                                                                                                                               ccAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCg                   <  2:871086/138‑1 (MQ=255)                                                                                                                                               ccAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCg                   >  1:871086/1‑138 (MQ=255)                                                                                                                                                 aaCCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGAAGCTTGAg     >  2:230934/1‑150 (MQ=255)                                                                                                                                                 aaCCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGAAGCTTGAg     <  1:1386120/150‑1 (MQ=255)                                                                                                                                                    cGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATg              <  2:273657/138‑1 (MQ=255)                                                                                                                                                    cGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATg              >  1:273657/1‑138 (MQ=255)                                                                                                                                                    cGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGAAGCTTGAGGGc  >  1:713722/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                       ATTAGCCGTGGCAAAGGCTGGCTGGCCGATTATGGCAGCGCCCTGCGTGACCAGGCACTGATCCTGGCGCTGCTGCAGGAAAGCAACCTGGCCAGCAGCCAGGTCGATCAGCGCCTGTTCGCCTTGTCGGATGAACTGGCAGCCAACCGTTGGCTGTCGACCCAAGAGCGCAATGCACTGTTCCTCG‑CTGGCCGTGGCCTGCTGGGCAAGCCGGAAGGCCAGTGGCAGGCACGCCTGGACAGCGCGGGTGAAGTGCGCGAGTTGAACAATGCCGAGGCCGGTATGAAGCTTGAGGGC  >  NC_002947/667179‑667475

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