Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,776,703	C→T	7.2%	G340G (GGC→GGT)	betA‑II →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,776,703	0	C	T	7.2%	205.2 / 3.1	69	G340G (GGC→GGT)	betA‑II	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG  >  NC_002947/5776559‑5776847                                                                                                                                                   |                                                                                                                                                  gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCGCCGATACCGGTGCCa                                                                                                                                               <  1:401681/149‑1 (MQ=255) gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                               <  8:32034/149‑1 (MQ=255) gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                               >  2:197130/1‑149 (MQ=255) gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                               <  2:311941/149‑1 (MQ=255)     gTTGGTGAAAACCTTCAGGACCACCTGGAACTGTCCCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTTGATGGTCAATGGCACCGGTATCGGCGCCAGCaa                                                                                                                                           <  1:259929/149‑1 (MQ=255)           gAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGtt                                                                                                                                     >  6:325313/1‑149 (MQ=255)           gAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGtt                                                                                                                                     >  7:281132/1‑149 (MQ=255)            aaaaCCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTc                                                                                                                                    <  8:231230/149‑1 (MQ=255)                ccTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAgg                                                                                                                                >  5:363960/1‑149 (MQ=255)                     aGGAGCACGTGGAAGTGTACCTTCAGTTCGCCTGCACCCTGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           <  6:270615/149‑1 (MQ=255)                     aGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           <  7:240182/149‑1 (MQ=255)                       gACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGAAGGTTTCGCTATACCCTTCGCTGCTGTGGTGGAACCAGCCGGC‑‑TATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTAGAGGCAggggg                                                                                                                         >  5:385096/1‑149 (MQ=255)                          caccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGtt                                                                                                                      >  4:144/1‑149 (MQ=255)                          caccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGtt                                                                                                                      <  6:363960/149‑1 (MQ=255)                            ccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCa                                                                                                                    >  2:130853/1‑149 (MQ=255)                             cTGGAACTGTACCTGCAGTACGCCTGCACCCAGCTGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCAt                                                                                                                   >  8:90470/1‑149 (MQ=255)                             cTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCAt                                                                                                                   >  1:234150/1‑149 (MQ=255)                               ggAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATcc                                                                                                                 >  4:336373/1‑149 (MQ=255)                                  aCTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGa                                                                                                              <  8:119289/149‑1 (MQ=255)                                   cTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGAc                                                                                                             <  2:234150/149‑1 (MQ=255)                                   cTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGAc                                                                                                             <  8:80144/149‑1 (MQ=255)                                        ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           <  2:90254/130‑1 (MQ=255)                                        ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggg                                                                                                                           >  1:90254/1‑130 (MQ=255)                                         cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCAGGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGt                                                                                                         >  1:120398/1‑149 (MQ=255)                                         cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggggg                                                                                                                         >  8:122886/1‑131 (MQ=255)                                         cTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCggggg                                                                                                                         <  7:122886/131‑1 (MQ=255)                                            cAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCg                                                                                                      <  5:416458/147‑1 (MQ=255)                                                 cGCCTGCACCCAGCTGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCa                                                                                              <  7:90470/149‑1 (MQ=255)                                                       cACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCTGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGc                                                                                        <  5:222700/149‑1 (MQ=255)                                                          ccAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATcc                                                                                                                 <  4:348812/122‑1 (MQ=255)                                                          ccAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATcc                                                                                                                 >  3:348812/1‑122 (MQ=255)                                                          ccAGCCGGCTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCg                                                                                      >  8:145233/1‑148 (MQ=255)                                                               ctgTTTCGCTATACCCGTCGCTGCTGTTGTGGAACAAGCCGGC‑‑GATAGGTGCAGAGTGGATGTTCAATGTCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATc                                                                                <  6:385096/147‑1 (MQ=255)                                                                  cttCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAg                                                                             <  7:145233/148‑1 (MQ=255)                                                                      gCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTAc                                                                          <  1:327068/148‑1 (MQ=255)                                                                      gCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTAc                                                                          >  3:332315/1‑148 (MQ=255)                                                                            cccGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTc                                                                    <  8:319684/148‑1 (MQ=255)                                                                               gTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGt                                                                                           >  8:220497/1‑122 (MQ=255)                                                                               gTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGt                                                                                           <  7:220497/122‑1 (MQ=255)                                                                                   ctgctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGt                                                            <  1:197130/149‑1 (MQ=255)                                                                                    tgctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGTCCGGGGGGTTCATCCGGACTCGTGCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGtt                                                           >  6:278532/1‑149 (MQ=255)                                                                                    tgctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCgggg                                                                                                                          >  2:191095/1‑87 (MQ=255)                                                                                     gctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTg                                                          >  6:290737/1‑149 (MQ=255)                                                                                      ctgTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTg                                                          <  6:76206/148‑1 (MQ=255)                                                                                            tggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTc                                                                    >  8:148308/1‑132 (MQ=255)                                                                                            tggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTc                                                                    <  7:148308/132‑1 (MQ=255)                                                                                            tggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTaa                                                   >  6:108396/1‑149 (MQ=255)                                                                                             ggAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                  >  8:333328/1‑149 (MQ=255)                                                                                                aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      <  6:170893/142‑1 (MQ=255)                                                                                                aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      <  6:164829/142‑1 (MQ=255)                                                                                                aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      >  5:170893/1‑142 (MQ=255)                                                                                                aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                      >  5:164829/1‑142 (MQ=255)                                                                                                aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTa                                                >  6:20747/1‑148 (MQ=255)                                                                                                      cGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGc                                         <  4:186281/149‑1 (MQ=255)                                                                                                       ggC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCt                                        <  3:336373/149‑1 (MQ=255)                                                                                                                    ccGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGt                              <  3:144/148‑1 (MQ=255)                                                                                                                               ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAAt                                   >  3:338008/1‑132 (MQ=255)                                                                                                                               ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAAt                                   <  4:338008/132‑1 (MQ=255)                                                                                                                                tCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                  <  4:107950/116‑1 (MQ=255)                                                                                                                                tCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                                  >  3:107950/1‑116 (MQ=255)                                                                                                                                  aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                       <  1:283009/110‑1 (MQ=255)                                                                                                                                  aaTGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGa                                                       >  2:283009/1‑109 (MQ=255)                                                                                                                                    tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGc                                >  6:240050/1‑130 (MQ=255)                                                                                                                                    tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc              >  6:38654/1‑148 (MQ=255)                                                                                                                                    tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc              >  1:119725/1‑148 (MQ=255)                                                                                                                                       cgccGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAgg          <  5:38654/147‑1 (MQ=255)                                                                                                                                        aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc         >  7:280933/1‑149 (MQ=255)                                                                                                                                          cGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgcg        >  5:214522/1‑148 (MQ=255)                                                                                                                                               tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATg  <  8:281132/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATG  >  NC_002947/5776559‑5776847

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