Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	353,067	A→G	intergenic (+34/+84)	hisF → / ← cbcV	imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	353,067	0	A	G	100.0%	195.8 / NA	59	intergenic (+34/+84)	hisF/cbcV	imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/24);  total (35/24)

TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCAC  >  NC_002947/352921‑353210                                                                                                                                                    |                                                                                                                                                 ttCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                                <  3:146028/149‑1 (MQ=255)      gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                              <  1:153522/146‑1 (MQ=255)      gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                             <  4:150437/147‑1 (MQ=255)      gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                             <  4:369953/147‑1 (MQ=255)       aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGcc                                                                                                                                          >  7:242338/1‑149 (MQ=255)         gggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCg                                                                                                                                        >  4:258495/1‑149 (MQ=255)             aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                    >  5:305883/1‑149 (MQ=255)             aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGATGATTTGCAGCCCTTAGCGGGCACGccccaacc                                                                                                                                    >  8:63989/1‑144 (MQ=255)                 cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCc‑a                                                                                                                                 >  4:172132/1‑147 (MQ=255)                    gaagatGCTGGCCGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCATTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCC‑ACAAGt                                                                                                                            >  1:418748/6‑149 (MQ=38)                         tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCt                                                                                                                       <  6:305883/149‑1 (MQ=255)                              gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAg                                                                                                                              <  1:317142/137‑1 (MQ=255)                              gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAg                                                                                                                              >  2:317142/1‑137 (MQ=255)                                ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGCAGGCCAAGGCCTACATCGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGCGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCC‑ACAGGTTTGCTGAGAtct                                                                                                                >  7:101152/1‑149 (MQ=39)                                     gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCtgtg                                                                                                           <  3:172132/149‑1 (MQ=255)                                      cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTgg                                                                                                          <  8:160079/149‑1 (MQ=255)                                          ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAg                                                                                                       >  1:111742/1‑149 (MQ=255)                                                tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                >  6:442142/1‑149 (MQ=255)                                                      gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                          >  5:327081/1‑149 (MQ=255)                                                           cacaGTGCCGGAGGCCAAGGCCTACATGGATTA‑GCGACGGATCGACGTTCGCTGATCCATTGACCGGGGCAGCTTTGCAACCCTTCGCGGGACCGCCCGCTACC‑CCAAGTATGCTGAGATCTATGTGGGAGCGGGCGTGCCCGCGaa                                                                                       >  6:121848/1‑147 (MQ=255)                                                            acaGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTAGCTGATCAATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCA‑ACAGGTTTTCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGACggg                                                                                    >  8:342603/1‑149 (MQ=38)                                                              aGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGGTTGCTGAGCTGTCTGTGGGAGCGGGCGTGCCCGCGCAGGGcc                                                                                  >  1:176318/1‑149 (MQ=38)                                                                  ccGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCACCCAACAGG‑TTGGAGAGGGATAGGTGGGAGCCGGCGTGCCCGGGCACGGCcgcg                                                                               >  7:160384/1‑148 (MQ=255)                                                                     gAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCC‑ACACGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCCCGAAGc                                                                           >  3:334255/1‑149 (MQ=39)                                                                       ggCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                         <  8:242338/149‑1 (MQ=255)                                                                         ccAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                                >  2:50964/1‑77 (MQ=255)                                                                         ccAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                                <  1:50964/77‑1 (MQ=255)                                                                         ccAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCT‑CC‑ACAGGTTTGCTGAGATCTCAGTGGGAGCGGGCGTGCCCCCGAAGGGCCGCGAACCGGccc                                                                      >  2:18435/1‑149 (MQ=39)                                                                          cAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                      >  8:328567/1‑149 (MQ=255)                                                                          cAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCCCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATATCTGTGGGACCGGGCGTGCCCCCGAAGGGCCCCGAAGCGGccc                                                                      >  6:322785/1‑149 (MQ=38)                                                                            aGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGGTTGCTGAGAGGGCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGgccc                                                                     >  7:342125/1‑148 (MQ=38)                                                                                cTACTTGGCTGC‑GCGGGGGCTCGTGGTTCGCTGATCCATTGACCGGGGGTCCTTTGCAGCCCTTGGCGGGCACGCCGGCTCCC‑ACAGGTTTGCTGATATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTg                                                                <  8:86352/149‑1 (MQ=37)                                                                                             gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                    >  3:67911/1‑149 (MQ=255)                                                                                             gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTACTTCAACCGTgg                                                    <  7:138127/149‑1 (MQ=255)                                                                                             gcgcgGGAGCGTCGTGCGCTCATCCATTGACCGGGGCTGCTTTCCAGCCCTGCGCGGCCACGCCCGCGCCC‑ACCGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                     <  3:258495/148‑1 (MQ=37)                                                                                                    aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                          >  2:51322/1‑104 (MQ=255)                                                                                                    aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                          <  1:51322/104‑1 (MQ=255)                                                                                                    aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                             >  7:57765/1‑149 (MQ=255)                                                                                                              gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcgc                                   >  1:12139/1‑149 (MQ=255)                                                                                                                  aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                                <  2:111742/148‑1 (MQ=255)                                                                                                                      aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                     >  1:374262/1‑123 (MQ=255)                                                                                                                      aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                     <  2:374262/123‑1 (MQ=255)                                                                                                                          aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTACACCGTGGTTCTTTCCGAGCAGCGCGTGGTAAAGCtg                       >  8:302809/1‑148 (MQ=255)                                                                                                                           gaagggCTGCAAAGCAGCCCCTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGt                                 <  3:311912/135‑1 (MQ=38)                                                                                                                           ccGTGGCTGCTTTGCAGCCGTTAGCGTGCACGCCCGCTCCC‑ACAGGTTTGCTGTGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                      <  4:334255/149‑1 (MQ=38)                                                                                                                           ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcg                                                                                           <  3:327136/80‑1 (MQ=255)                                                                                                                           ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcg                                                                                           >  4:327136/1‑80 (MQ=255)                                                                                                                           ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                      <  7:234279/149‑1 (MQ=255)                                                                                                                                gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATc                 <  7:342603/149‑1 (MQ=255)                                                                                                                                gctgctTTGCAGCCCTGCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATc                 <  7:373456/149‑1 (MQ=255)                                                                                                                                  tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ATAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGGTCTTGCCGCGCAGCGCGTGGTTAAGCTCGGTATccc               >  7:336167/1‑149 (MQ=38)                                                                                                                                  tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc               >  1:63418/1‑149 (MQ=255)                                                                                                                                  tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc               >  8:250860/1‑149 (MQ=255)                                                                                                                                  ggttttGCAGCCCTTCGCGGCCACGCCCTCTCCC‑ACTGTTTTGCTGAGTTCTCTGTGGGCGCGGGGGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc               <  8:160384/146‑1 (MQ=255)                                                                                                                                   gctTTGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                      <  3:280500/141‑1 (MQ=255)                                                                                                                                       tGCAGCCCTTCGCGGGCACGCCCGCTCCC‑ACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt                    <  7:392626/139‑1 (MQ=255)                                                                                                                                       tGCAGCCCTTCGCGGGCACGCCCCCGCCC‑ACAGGTTTGGGGAGATGTCTGTGGGGGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCCAATGGTTCAACCGGGGGTCTTGCCGCGCCGCGCGGGGGTAAGCGCggg                    >  8:392626/1‑138 (MQ=19)                                                                                                                                         cAGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCGGAGATGTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGGAAAGCTCGGTATCCCCGAGTAt        >  5:167933/1‑149 (MQ=255)                                                                                                                                          aGCCCTTCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGCGCAGCGCGTGGTAAAGCTCGGTATGCCCGAGTATc       >  5:123604/1‑149 (MQ=255)                                                                                                                                           gCCCTTTGCGGTCACGCCCGCTCCG‑ACTGGTTTGCTGAGATCTGTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGTAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcc      <  8:101152/149‑1 (MQ=38)                                                                                                                                               ttCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccac  >  2:82453/1‑149 (MQ=255)                                                                                                                                               ttCGCGGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccac  >  4:51090/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTC‑GCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCC‑ACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCAC  >  NC_002947/352921‑353210

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