Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	665,907	A→G	5.8%	S925G (AGC→GGC)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	665,907	0	A	G	5.8%	256.9 / 4.5	86	S925G (AGC→GGC)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (48/33);  new base G (3/2);  total (51/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TTGTCGGCGGTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAAC  >  NC_002947/665760‑666053                                                                                                                                                     |                                                                                                                                                   ttgtCGGCGGTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCa                                                                                                                                                    <  4:227167/148‑1 (MQ=255) ttgtCGGCGGTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAg                                                                                                                                                   >  6:249220/1‑149 (MQ=255)  tgtCGGCGGTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGc                                                                                                                                                  >  6:133394/1‑149 (MQ=255)  atctggcggTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGCTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGc                                                                                                                                                  <  7:357654/145‑1 (MQ=255)      ggcggTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTgg                                                                                                                                              >  6:363004/1‑149 (MQ=255)       gcggTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGATGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGc                                                                                                                                                  >  8:357654/1‑144 (MQ=255)         ggTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGt                                                                                                                                           >  7:231899/1‑149 (MQ=255)         ggTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGt                                                                                                                                           >  8:294113/1‑149 (MQ=255)             gACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCgg                                                                                                                                       <  6:233454/149‑1 (MQ=255)             gACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCgg                                                                                                                                       <  6:228083/149‑1 (MQ=255)             gACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCgg                                                                                                                                       <  5:133394/149‑1 (MQ=255)             gACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCgg                                                                                                                                       >  4:85840/1‑149 (MQ=255)                          tCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATgcggc                                                                                                                           >  7:310487/1‑148 (MQ=255)                          tCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATgcggc                                                                                                                           >  3:296317/1‑148 (MQ=255)                              cATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATgg                                                                                                                      >  1:227707/1‑149 (MQ=255)                              cATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATgg                                                                                                                      <  8:337872/149‑1 (MQ=255)                               aTCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGc                                                                                                                     >  5:19556/1‑149 (MQ=255)                                tCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGc                                                                                                                     >  5:200638/1‑148 (MQ=255)                                  aCCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGcc                                                                                                                    >  8:42057/1‑147 (MQ=255)                                  aCCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCaa                                                                                                                  <  5:363004/149‑1 (MQ=255)                                    cGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAgg                                                                                                                >  6:201341/1‑149 (MQ=255)                                      aTTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAgggg                                                                                                              <  6:262043/149‑1 (MQ=255)                                       tttCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCACCTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGt                                                                                                                                      <  3:66600/124‑1 (MQ=255)                                           aaaaCCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGc                                                                                                           >  6:38566/1‑147 (MQ=255)                                                    gACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCaa                                                                                                >  1:38701/1‑149 (MQ=255)                                                      cgatttCGACAGCCTGTTCTTGTTCAAGGCCTATGGTGACGATCAATTGGTCATCTACTGCGAGTTGATGTAGTACGGCCAAGGCCGCC‑TGGCCAGCGTGGCGTTCTGTGTCGATGAGGGGTTGGCCAAGTGGGGCAAGCGAGCCAact                                                                                              <  8:310487/145‑2 (MQ=255)                                                        cGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAacac                                                                                             >  4:377755/1‑148 (MQ=255)                                                         gTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTa                                                                                           >  6:66595/1‑149 (MQ=255)                                                         gTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTa                                                                                           >  7:173934/1‑149 (MQ=255)                                                          ttCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGAt                                                                                                                        >  1:359554/1‑119 (MQ=255)                                                          ttCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGAt                                                                                                                        <  2:359554/119‑1 (MQ=255)                                                          ttCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTAc                                                                                          >  4:102243/1‑149 (MQ=255)                                                             gCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACgg                                                                                        <  5:249220/148‑1 (MQ=255)                                                             gCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGt                                                                                       <  8:192186/149‑1 (MQ=255)                                                              ccAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGt                                                                                       >  5:39642/1‑148 (MQ=255)                                                               cAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACg                                                                                         >  7:347152/1‑145 (MQ=255)                                                               cAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACg                                                                                         <  8:347152/145‑1 (MQ=255)                                                                   cTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCa                                                                                  >  7:32336/1‑148 (MQ=255)                                                                    tGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCa                                                                                  >  2:5584/1‑147 (MQ=255)                                                                      ttCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATa                                                                                <  7:296520/147‑1 (MQ=255)                                                                      ttCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGt                                                                              >  3:191781/1‑149 (MQ=255)                                                                          ggCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGcc                                                                          >  4:89752/1‑149 (MQ=255)                                                                           gcgcAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGccc                                                                         <  7:42057/149‑1 (MQ=255)                                                                            cgcAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCTACCCAACACTACGGTAACCATAGTCGCCca                                                                        <  6:39642/149‑1 (MQ=255)                                                                               aaGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGTCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCaa                                                                                                >  7:267837/1‑122 (MQ=255)                                                                               aaGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGTCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCaa                                                                                                <  8:267837/122‑1 (MQ=255)                                                                                  gCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCagag                                                                  >  8:48534/1‑149 (MQ=255)                                                                                    cTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGcc                                                                >  4:221174/1‑149 (MQ=255)                                                                                       tGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTg                                                              <  3:89752/148‑1 (MQ=255)                                                                                         gTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCg                                                           <  8:173934/149‑1 (MQ=255)                                                                                         gTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAACGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCg                                                           >  3:374913/1‑149 (MQ=255)                                                                                             cGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGa                                                       <  6:196143/149‑1 (MQ=255)                                                                                              gATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGAc                                                      >  6:206713/1‑149 (MQ=255)                                                                                              gATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGAc                                                      <  7:294113/149‑1 (MQ=255)                                                                                                       gACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGgacga                                             >  1:110688/1‑149 (MQ=255)                                                                                                          aTCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAg                                                                               <  5:320392/112‑1 (MQ=255)                                                                                                          aTCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAg                                                                               >  6:320392/1‑112 (MQ=255)                                                                                                            cTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAgg                                        >  5:245142/1‑149 (MQ=255)                                                                                                               cGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAgg                                        <  5:66595/146‑1 (MQ=255)                                                                                                               cGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGt                                      <  7:48534/148‑1 (MQ=255)                                                                                                                 gCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGCCGCTGGACGACAAGGGTGAg                                   >  8:303199/1‑149 (MQ=255)                                                                                                                 gCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAg                                   <  5:201341/149‑1 (MQ=255)                                                                                                                  ccAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAgg                                  >  8:281762/1‑149 (MQ=255)                                                                                                                  acttttGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAgg                                  <  8:279052/145‑1 (MQ=255)                                                                                                                   cAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAggg                                 >  4:108110/1‑149 (MQ=255)                                                                                                                         aTCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAgg                            <  7:281762/148‑1 (MQ=255)                                                                                                                         aTCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGGGACGCTGGACGACAAGGGTGAGGGCCAGGc                           >  8:296683/1‑149 (MQ=255)                                                                                                                           ccaGGCTGGCC‑AGGCGGCCTTGGCCGGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAgg                                  <  7:344076/138‑1 (MQ=255)                                                                                                                           ccaGGCTGGCC‑AGGCGGCCTTGGCCGGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAgg                                  >  8:344076/3‑140 (MQ=255)                                                                                                                           cGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTa                         >  4:150832/1‑149 (MQ=255)                                                                                                                             aGGCTGGCC‑AGGCGGCCTTGGCCGGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACg                       >  2:42604/1‑149 (MQ=255)                                                                                                                             aGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACg                       >  1:156715/1‑149 (MQ=255)                                                                                                                                  ggccAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGa                  <  4:296317/149‑1 (MQ=255)                                                                                                                                  ggccAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGACTGCCGGTGACGCTGGACGCCAAGGGTGAGGGCCAGGCTACGGGCGa                  >  8:357429/1‑149 (MQ=255)                                                                                                                                   gccAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGa                  >  1:125484/1‑148 (MQ=255)                                                                                                                                    caccgaacgcc‑AGGCTGGCCAGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACAACAAGGGTGAGGGCCAGGCTACGGTc                    <  1:23752/138‑1 (MQ=255)                                                                                                                                    caccgaacgcc‑AGGCTGGCCAGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACAACAAGGGTGAGGGCCAGGCTACGGTc                    >  2:23752/8‑145 (MQ=255)                                                                                                                                     cAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACAt               >  6:335204/1‑149 (MQ=255)                                                                                                                                     cAAGGCCGCC‑TGGCCAGCCTGGCGTTCGCTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACAt               >  1:219880/1‑149 (MQ=255)                                                                                                                                        ggCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCt                                                   >  6:83654/1‑110 (MQ=255)                                                                                                                                        ggCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCt                                                   <  5:83654/110‑1 (MQ=255)                                                                                                                                        ggCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATccc            >  5:267942/1‑149 (MQ=255)                                                                                                                                         gccgcc‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAATGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATcccc           <  4:374913/149‑1 (MQ=255)                                                                                                                                           cgcc‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGa         <  5:198617/149‑1 (MQ=255)                                                                                                                                             cc‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACt       <  6:245142/149‑1 (MQ=255)                                                                                                                                                  gCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCaa   <  5:206713/149‑1 (MQ=255)                                                                                                                                                   ccAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAAc  >  6:327681/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TTGTCGGCGGTGGACGTGGGCATCCTCAACATCACCGATTTCAAAACCCCCGACCCGTTCGCCAGCCTGTTCGGGCGCAAGGCCTATGGTGCCGATCAACTGGACATCTACGGCCAGTTGATCGAGGCCGGCCAAGGCCGCC‑TGGCCAGCCTGGCGTTCGGTGGCGATGCGGCGATGGCCAAGGGGGGCAAGCGACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAAC  >  NC_002947/665760‑666053

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