Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	680,307	G→A	6.1%	L398L (CTG→CTA)	PP_0582 →	thiolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	680,307	0	G	A	6.1%	196.1 / 2.4	65	L398L (CTG→CTA)	PP_0582	thiolase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base A (2/2);  total (26/39)
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.68e-01

TGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAGCTGT  >  NC_002947/680160‑680446                                                                                                                                                      |                                                                                                                                            tGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGt                                                                                                                                            >  7:194786/1‑149 (MQ=255)   ccGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGTGGATATTGGCCAACCTGGCCAAGCTGtt                                                                                                                                           <  1:213044/148‑1 (MQ=255)   ccGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGtt                                                                                                                                           <  2:53092/148‑1 (MQ=255)   ccGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTg                                                                                                                                          <  6:288473/149‑1 (MQ=255)    cGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGTGGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTg                                                                                                                                          <  1:288749/149‑1 (MQ=255)    cGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTgg                                                                                                                                         <  6:3214/149‑1 (MQ=255)              cGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCggg                                                                                                                               <  8:307878/149‑1 (MQ=255)                  aaCGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAg                                                                                                                           <  2:239134/149‑1 (MQ=255)                    ctgTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCTGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAAGGGCGGGTAAggg                                                                                                                         <  6:298435/147‑1 (MQ=255)                      gTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGgc                                                                                                                       >  3:13675/1‑149 (MQ=255)                          ggCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCgg                                                                                                                   >  3:7433/1‑149 (MQ=255)                            ccTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGc                                                                                                                 <  8:194786/149‑1 (MQ=255)                               ggAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAggg                                                                                                                         <  6:374015/138‑1 (MQ=255)                               ggAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAggg                                                                                                                         >  5:374015/1‑138 (MQ=255)                                  ggCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGAttt                                                                                                           <  8:247159/149‑1 (MQ=255)                                      ccgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGAt                                                                                                       <  5:220844/149‑1 (MQ=255)                                       cgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATc                                                                                                      <  7:222251/149‑1 (MQ=255)                                         cTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCtg                                                                                                    >  2:365104/1‑149 (MQ=255)                                          ttGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCTGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCtgt                                                                                                   <  8:139700/149‑1 (MQ=255)                                            ggCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGTCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCtgt                                                                                                   <  8:187290/147‑1 (MQ=255)                                              cgcgATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTATTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCgg                                                                                               >  3:215784/1‑149 (MQ=255)                                              cgcgATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTATTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCgg                                                                                               >  3:220441/1‑149 (MQ=255)                                                       ccGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGc                                                                                              <  1:160712/141‑1 (MQ=255)                                                       ccGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGc                                                                                              >  2:160712/1‑141 (MQ=255)                                                        cGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAAGCCGGTGGGCGGAAATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGc                                                                                      >  1:294510/1‑149 (MQ=255)                                                            gcaAACTCAATGTCAGGGGTAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                                 >  3:68489/1‑149 (MQ=255)                                                            gcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTATTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                                 <  4:220441/149‑1 (MQ=255)                                                            gcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTATTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                                 <  4:215784/149‑1 (MQ=255)                                                             caAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGAt                                                                                                       <  8:230985/126‑1 (MQ=255)                                                             caAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGAt                                                                                                       >  7:230985/1‑126 (MQ=255)                                                             caAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCtgt                                                                                                   >  2:98858/1‑130 (MQ=255)                                                             caAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCtgt                                                                                                   <  1:98858/130‑1 (MQ=255)                                                              ataCTCACTGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGCTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGt                                                                               <  8:314154/147‑1 (MQ=255)                                                                aCTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGa                                                                             >  7:135026/1‑149 (MQ=255)                                                                aCTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGa                                                                             >  5:37207/1‑149 (MQ=255)                                                                  tCAATGTCAGGGGCAGTTCGCTGGCGGCAG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTggct                                                                           >  6:340881/1‑146 (MQ=255)                                                                    aaTGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                           <  4:13675/147‑1 (MQ=255)                                                                    aaTGTCAGGGGCAGTTCGCCGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGAcc                                                                           <  8:57467/147‑1 (MQ=255)                                                                     aTGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGtt                                                                        <  8:135026/149‑1 (MQ=255)                                                                        tCAGGGGTAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAt                                                                      <  4:68489/148‑1 (MQ=255)                                                                           ggggCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAtcgtc                                                                  <  2:294510/149‑1 (MQ=255)                                                                             ggCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGa                                                                <  2:193519/149‑1 (MQ=255)                                                                                  ttCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCt                                                           <  4:7433/149‑1 (MQ=255)                                                                                     gCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTg                                                          >  6:40079/1‑147 (MQ=255)                                                                                      cTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATGCTCGAGCGCTGAAg                                                       >  8:116604/1‑149 (MQ=255)                                                                                         gcggcCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTAc                                                    >  6:257172/1‑149 (MQ=255)                                                                                                    aCCCGTTTGCGGCAA‑CCGGTGGGCGGATTTTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGGGCGGGCTGATTTCGATCTGTGCGGCTGGCGGGGAGGGGGTGACCGTTTTGGTCGAGGGCTGAAGTACCTCCAGGGt                                           >  6:375472/1‑148 (MQ=255)                                                                                                     cccGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTaa                                         <  7:304052/149‑1 (MQ=255)                                                                                                      ccGTTTGCTGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGAGTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAAc                                        >  2:5218/1‑149 (MQ=255)                                                                                                             cGGCAA‑CCGGTTGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCaga                                  <  2:303820/148‑1 (MQ=255)                                                                                                              ggCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAgg                                >  7:316237/1‑149 (MQ=255)                                                                                                                           gCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAg                                                               >  2:341682/1‑106 (MQ=255)                                                                                                                           gCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAg                                                               <  1:341682/106‑1 (MQ=255)                                                                                                                           gCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCt                    <  6:37207/149‑1 (MQ=255)                                                                                                                             ggATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAg                                              >  2:178664/1‑121 (MQ=255)                                                                                                                             ggATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAg                                              <  1:178664/121‑1 (MQ=255)                                                                                                                               atatTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTa                                                     >  5:215314/1‑112 (MQ=255)                                                                                                                               atatTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTa                                                     <  6:215314/112‑1 (MQ=255)                                                                                                                               agcttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGc                      >  5:193035/4‑143 (MQ=255)                                                                                                                               agcttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGc                      <  6:193035/140‑1 (MQ=255)                                                                                                                                gagttgCCCACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCGGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATc               <  8:251209/144‑1 (MQ=255)                                                                                                                                  ttGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCa              <  7:296314/148‑1 (MQ=255)                                                                                                                                    ggCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATg            <  8:316237/148‑1 (MQ=255)                                                                                                                                          ccTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAGc     <  5:40079/149‑1 (MQ=255)                                                                                                                                               ccAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAGCTGt  >  3:255829/1‑147 (MQ=255)                                                                                                                                                      |                                                                                                                                            TGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCG‑GGCACCCGTTTGCGGCAA‑CCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGATCCAGCTGT  >  NC_002947/680160‑680446

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