Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	66,208	Δ1 bp	6.7%	coding (1266/1653 nt)	betA‑I →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	66,208	0	G	.	6.7%	186.2 / 6.7	60	coding (1266/1653 nt)	betA‑I	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base . (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACC  >  NC_002947/66062‑66355                                                                                                                                                   |                                                                                                                                                     ggTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGAATTGCg                                                                                                                                                    <  1:112964/149‑1 (MQ=255)   tGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCggg                                                                                                                                                  <  7:164521/149‑1 (MQ=255)               tCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGTCCAGACGCGGGAGCATCGACATCTGCAGCACCGACATGAACAGAACGCCGCGGATCGACCCCAATTACATCAGTGCTCAGCAGGACTTGCGGGTTGATGCCGATg                                                                                                                                      >  7:223232/1‑149 (MQ=255)                  gTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTg                                                                                                                                                      <  5:168724/130‑1 (MQ=255)                  gTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTg                                                                                                                                                      >  6:168724/1‑130 (MQ=255)                      ccGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATg                                                                                                                                      <  6:36056/142‑1 (MQ=255)                      ccGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATg                                                                                                                                      >  5:36056/1‑142 (MQ=255)                      ccGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCg                                                                                                                               >  1:197930/1‑149 (MQ=255)                        ggCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGcc                                                                                                                             <  2:187601/149‑1 (MQ=255)                                cTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCgg                                                                                                                     <  1:135725/149‑1 (MQ=255)                                   ccTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGTATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGgc                                                                                                                    <  2:197930/147‑1 (MQ=255)                                    cTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCa                                                                                                                 <  2:140265/149‑1 (MQ=255)                                      cGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATc                                                                                                               >  1:235557/1‑149 (MQ=255)                                                 cTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGcccc                                                                                                    >  5:84559/1‑149 (MQ=255)                                                       ccGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCt                                                                                              >  8:100942/1‑149 (MQ=255)                                                       ccGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCt                                                                                              >  3:63893/1‑149 (MQ=255)                                                          gCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTcgc                                                                                           <  6:214628/149‑1 (MQ=255)                                                               ccGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCAGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCtt                                                                                      >  1:198531/1‑149 (MQ=255)                                                                gccgGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCtt                                                                                      <  4:42903/146‑1 (MQ=255)                                                                 gcgGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTc                                                                                    <  2:235557/149‑1 (MQ=255)                                                                  cgGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTc                                                                                    >  2:155788/1‑148 (MQ=255)                                                                         atagACATATTCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTTCTCCGCAGGAATTGCGGTTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCaa                                                                            <  8:223232/146‑1 (MQ=255)                                                                         aTCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGCCCCCAATTACCTCAGTCCTCCGCAGGACTTGCCGCTTGCTGCCGATGCCATTCCCCTTACCAGGCGCATCGTCCAGCCCCCTGCCCACGCCGCTTTCGAGccc                                                                              >  1:124615/1‑147 (MQ=255)                                                                         aTCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCaa                                                                            >  2:130380/1‑149 (MQ=255)                                                                             aCATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTc                                                                                             >  7:159134/1‑128 (MQ=255)                                                                             aCATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTc                                                                                             <  8:159134/128‑1 (MQ=255)                                                                                    cagcaCCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGc                                                                 <  6:84559/149‑1 (MQ=255)                                                                                     agcaCCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCtt                                                                                      <  6:4447/127‑1 (MQ=255)                                                                                     agcaCCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCtt                                                                                      >  5:4447/1‑127 (MQ=255)                                                                                        aCCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCa                                                                                                          >  4:239498/1‑104 (MQ=255)                                                                                        aCCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCa                                                                                                          <  3:239498/104‑1 (MQ=255)                                                                                          cGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGcc                                                          <  4:63893/149‑1 (MQ=255)                                                                                              aTGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGgc                                                      >  7:100166/1‑149 (MQ=255)                                                                                                     gCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAatat                                                                     <  4:123612/128‑1 (MQ=255)                                                                                                     gCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAatat                                                                     >  3:123612/1‑128 (MQ=255)                                                                                                       aCGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAg                                             >  1:134831/1‑149 (MQ=255)                                                                                                         gccgcGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCaaa                                              >  2:173243/1‑146 (MQ=255)                                                                                                          ccgcGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCATGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGa                                          <  1:130380/149‑1 (MQ=255)                                                                                                            gcgGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTT‑‑GGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAgg                                        >  4:77117/1‑147 (MQ=255)                                                                                                            gcgGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTT‑‑GGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAgg                                        >  4:80905/1‑147 (MQ=255)                                                                                                              ggATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑gg                                                             <  5:174592/126‑1 (MQ=255)                                                                                                              ggATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑gg                                                             >  6:174592/1‑126 (MQ=255)                                                                                                              ggATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGaa                                      <  2:134831/149‑1 (MQ=255)                                                                                                                    aCCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGa                                          >  8:46798/1‑139 (MQ=255)                                                                                                                    aCCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGa                                          <  7:46798/139‑1 (MQ=255)                                                                                                                    aCCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTCCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTg                                >  6:163996/1‑149 (MQ=255)                                                                                                                     ccccAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGt                               >  5:224997/1‑149 (MQ=255)                                                                                                                       ccAATTACCTCAGTGCTCCGCAGGACTT‑‑GGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGa                           <  3:80905/149‑1 (MQ=255)                                                                                                                       ccAATTACCTCAGTGCTCCGCAGGACTT‑‑GGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGa                           <  3:77117/149‑1 (MQ=255)                                                                                                                        cAATTACCTCAGTGCTCCGCATGACTTTCGGGTTGCTGCCGATTCCATTCGCCTTACCCGGCGGATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGCGCCCAAGGAATATTTTCCTGGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTc                             <  1:173243/149‑1 (MQ=255)                                                                                                                        cAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGAcc                                  >  5:89494/1‑143 (MQ=255)                                                                                                                        cAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGAcc                                  <  6:89494/143‑1 (MQ=255)                                                                                                                         aaTTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTc                             >  2:102557/1‑147 (MQ=255)                                                                                                                           ttACCTCAGTGCTCCTCAGGACTTGGGGGTTGCTGCCGTTGCCATTCGCCTTACCCGGCGCATCGTCCAGGGCGCTGCCCTCGCCGCTTTCGAGCCCACGGTATTTTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGATGACCTGTTCGAgg                         <  2:124615/149‑1 (MQ=255)                                                                                                                               cTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGc                     >  7:178750/1‑149 (MQ=255)                                                                                                                               cTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGc                     >  5:13832/1‑149 (MQ=255)                                                                                                                                  aGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCt                    <  8:178750/147‑1 (MQ=255)                                                                                                                                  aGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCt                    <  6:13832/147‑1 (MQ=255)                                                                                                                                         cGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGc                 >  3:33799/1‑143 (MQ=255)                                                                                                                                           cAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCgg         >  2:12420/1‑149 (MQ=255)                                                                                                                                            aGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGc        >  7:113227/1‑149 (MQ=255)                                                                                                                                                cTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGcacca    <  2:198531/149‑1 (MQ=255)                                                                                                                                                   gCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCaccacc  <  7:100942/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGGCCAGCCGCGGGCGCATCGACATCTGCAGCACCGACATGAACAGCACGCCGCGGATCGACCCCAATTACCTCAGTGCTCCGCAGGACTTGCGGGTTGCTGCCGATGCCATTCGCCTTACCCGGCGCATCGTCCAGGCCCCTGCCCTCGCCGCTTTCGAGCCCAAGGAATATTTGCC‑GGGCCCGGCGCTGCAAAGCGAGGAAGACCTGTTCGAGGCCGCTGGCAAGATCGGCACCACC  >  NC_002947/66062‑66355

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