Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,776,696	G→A	7.2%	G338D (GGT→GAT)	betA‑II →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,776,696	0	G	A	7.2%	203.2 / 4.4	69	G338D (GGT→GAT)	betA‑II	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base A (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 = 7.97e-01

CACGACCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCG  >  NC_002947/5776548‑5776841                                                                                                                                                       |                                                                                                                                                   cACGACCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCgg                                                                                                                                                    >  7:319684/1‑149 (MQ=255)  aCGACCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCgg                                                                                                                                                    >  6:416458/1‑148 (MQ=255)     aCCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATc                                                                                                                                                <  4:87318/149‑1 (MQ=255)            gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCGCCGATACCGGTGCCa                                                                                                                                         <  1:401681/149‑1 (MQ=255)            gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                         <  2:311941/149‑1 (MQ=255)            gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                         <  8:32034/149‑1 (MQ=255)            gggCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCa                                                                                                                                         >  2:197130/1‑149 (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                                                                                                                <  6:363960/149‑1 (MQ=255)                                     caccTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGtt                                                                                                                >  4:144/1‑149 (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                                                                    >  3:332315/1‑148 (MQ=255)                                                                                 gCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTAc                                                                    <  1:327068/148‑1 (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                                                >  5:164829/1‑142 (MQ=255)                                                                                                           aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                >  5:170893/1‑142 (MQ=255)                                                                                                           aCCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGAt                                                <  6:164829/142‑1 (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                             <  4:338008/132‑1 (MQ=255)                                                                                                                                          ttCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAAt                             >  3:338008/1‑132 (MQ=255)                                                                                                                                           tCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                            >  3:107950/1‑116 (MQ=255)                                                                                                                                           tCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAAc                                            <  4:107950/116‑1 (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)                                                                                                                                                       |                                                                                                                                                   CACGACCTTCCGGGCGTTGGTGAAAACCTTCAGGACCACCTGGAACTGTACCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGC‑‑GATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAG‑TTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCG  >  NC_002947/5776548‑5776841

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