Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	680,295	C→G	5.7%	N394K (AAC→AAG) ‡	PP_0582 →	thiolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	680,295	0	C	G	5.7%	227.1 / 4.1	70	N394K (AAC→AAG) ‡	PP_0582	thiolase family protein
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base G (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

AGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAG  >  NC_002947/680147‑680442                                                                                                                                                       |                                                                                                                                                     aGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAAc                                                                                                                                                      <  3:390151/149‑1 (MQ=255)     cTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTg‑g                                                                                                                                                 <  6:219375/149‑1 (MQ=255)       gggAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑Gc                                                                                                                                                >  5:344240/1‑148 (MQ=255)          aaGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAg                                                                                                                                            <  2:3655/149‑1 (MQ=255)            gATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCt                                                                                                                                          >  2:151796/1‑149 (MQ=255)              tGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGt                                                                                                                                        >  8:38624/1‑149 (MQ=255)               gCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCaa                                                                                                                                             >  1:412737/1‑143 (MQ=255)               gCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCaa                                                                                                                                             <  2:412737/143‑1 (MQ=255)                  gATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGc                                                                                                                                    <  3:240855/149‑1 (MQ=255)                        tGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAAcggc                                                                                                                              >  4:236709/1‑149 (MQ=255)                        tGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAAcggc                                                                                                                              >  4:293102/1‑149 (MQ=255)                        tGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAAcggc                                                                                                                              >  1:307026/1‑149 (MQ=255)                          ccGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAAcggcgg                                                                                                                            <  8:75045/149‑1 (MQ=255)                          ccGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAAcggcgg                                                                                                                            <  1:125349/149‑1 (MQ=255)                          ccGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAAcggcgg                                                                                                                            >  2:317839/1‑149 (MQ=255)                             tGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGt                                                                                                                          >  6:6006/1‑148 (MQ=255)                                aCGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAg                                                                                                                       <  2:307026/148‑1 (MQ=255)                                       ggCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCgg                                                                                                               >  7:358437/1‑149 (MQ=255)                                         ccTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAGGTTG‑GCCaatatc                                                                                                                                         >  8:265766/1‑117 (MQ=255)                                         ccTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAGGTTG‑GCCaatatc                                                                                                                                         <  7:265766/121‑5 (MQ=255)                                          cTGGTGGCGCCGGTTGGCGGGATCGCCCGCAGCAAACTCAATGTCAGGGGCAGTTCGCGTGCGGCCGGGCACCCGTTTTCGGC‑AAACGGTCAGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGc                                                                                                             <  1:317839/149‑1 (MQ=255)                                           tGGAGGCGCCGCTTGGCGCGATCGACCGCAGGAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTg                                                                                                                                      >  2:30050/1‑122 (MQ=255)                                           tGGAGGCGCCGCTTGGCGCGATCGACCGCAGGAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTg                                                                                                                                      <  1:30050/122‑1 (MQ=255)                                             gAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGAt                                                                                                         >  2:89622/1‑149 (MQ=255)                                                gcgcCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGAttt                                                                                                       <  8:358437/148‑1 (MQ=255)                                                  gccgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCCGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGa                                                                                                    >  4:246737/1‑149 (MQ=255)                                                    cgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTGTGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGAt                                                                                                   <  1:89622/149‑1 (MQ=255)                                                            gcgATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCgg                                                                                           >  4:246389/1‑148 (MQ=255)                                                              gATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGTGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAg                                                                                        <  3:246737/149‑1 (MQ=255)                                                               aTCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGTGTAAGGGGCGCGGGCTGATTTCGAGCTGTGCGGCAgg                                                                                       >  7:137278/1‑149 (MQ=255)                                                                  gACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGaa                                                                                                   <  7:172120/134‑2 (MQ=255)                                                                  gACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGaa                                                                                                   >  8:172120/1‑133 (MQ=255)                                                                   aCCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGt                                                                                                                                        <  5:303656/96‑1 (MQ=255)                                                                   aCCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGt                                                                                                                                        >  6:303656/1‑96 (MQ=255)                                                                       cagcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCaa                                                                                                                                  <  8:390097/98‑1 (MQ=255)                                                                       cagcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCaa                                                                                                                                  >  7:390097/1‑98 (MQ=255)                                                                         gcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                             <  6:344240/149‑1 (MQ=255)                                                                         gcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                             >  5:378921/1‑149 (MQ=255)                                                                         gcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                             <  7:38624/149‑1 (MQ=255)                                                                           aaaCTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGt                                                                           >  5:293877/1‑149 (MQ=255)                                                                             aCTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGa                                                                         <  3:293102/149‑1 (MQ=255)                                                                               tCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                       >  6:404377/1‑149 (MQ=255)                                                                                     tCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAtc                                                                 <  3:246389/149‑1 (MQ=255)                                                                                        ggggCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAtcgtc                                                              <  3:236709/149‑1 (MQ=255)                                                                                          ggCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCTGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGa                                                            <  5:404377/149‑1 (MQ=255)                                                                                           gCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                       >  1:229465/1‑137 (MQ=255)                                                                                           gCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                       <  2:229465/137‑1 (MQ=255)                                                                                                   cTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAg                                                   <  6:293877/149‑1 (MQ=255)                                                                                                        ggcCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAg                                                                                >  2:294244/1‑115 (MQ=255)                                                                                                        ggcCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAg                                                                                <  1:294244/115‑1 (MQ=255)                                                                                                         gcCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTc                                             <  5:6006/149‑1 (MQ=255)                                                                                                              cccaccGGTTGCCGCAAACGGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGa                                                                         >  5:347206/5‑117 (MQ=255)                                                                                                              cccaccGGTTGCCGCAAACGGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGa                                                                         <  6:347206/113‑1 (MQ=255)                                                                                                               ccaccGGTTGCCGCAAACGGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAggg                                        >  3:95142/4‑149 (MQ=255)                                                                                                                  ccGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAgg                                         <  4:76666/144‑1 (MQ=255)                                                                                                                  ccGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAgg                                         >  3:76666/1‑144 (MQ=255)                                                                                                                  ccGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAAc                                    >  4:36589/1‑149 (MQ=255)                                                                                                                  ccGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAAc                                    >  4:208353/1‑149 (MQ=255)                                                                                                                        gCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCa                                >  4:377803/1‑147 (MQ=255)                                                                                                                                 cGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGAGCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGa                      <  1:151796/149‑1 (MQ=255)                                                                                                                                   gt‑gtGCGGATATTGGCCAACCTG‑GCAAAGCTGTTGGCAACGGCGGGTAAGTGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGt                    <  8:137278/149‑1 (MQ=255)                                                                                                                                        gCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCt                <  3:377803/149‑1 (MQ=255)                                                                                                                                          ggATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTc               >  4:4146/1‑148 (MQ=255)                                                                                                                                          ggATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTc               >  1:256232/1‑148 (MQ=255)                                                                                                                                             gcttGGCCAGGTTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTAc                                                <  6:341172/110‑1 (MQ=255)                                                                                                                                             gcttGGCCAGGTTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTAc                                                >  5:341172/3‑112 (MQ=255)                                                                                                                                                tGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                       >  6:405847/1‑86 (MQ=255)                                                                                                                                                tGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                       <  5:405847/86‑1 (MQ=255)                                                                                                                                                   ccAACCTG‑GCCAAGCTCTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATATCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTAATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATcatcata     >  8:106987/1‑148 (MQ=255)                                                                                                                                                     aaCCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATc     >  6:201864/1‑147 (MQ=255)                                                                                                                                                      aCCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAg  <  3:208353/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     AGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGC‑AACCGGT‑GGGCGGATATTGGCCAACCTG‑GCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAG  >  NC_002947/680147‑680442

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