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%	213.3 / NA	63	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 (31/32);  total (31/32)

GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CAGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352925‑353200                                                                                                                                                |                                                                                                                                        ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑Cggggccc                                                                                                                                   >  2:409401/1‑145 (MQ=255)  gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGc                                                                                                                                  >  8:101885/1‑149 (MQ=255)    gctggccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCCTTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGccc                                                                                                                                >  4:166631/8‑149 (MQ=255)    gctggccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGccc                                                                                                                                >  4:160511/8‑149 (MQ=255)        cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTc                                                                                                                            <  1:75420/149‑1 (MQ=255)        cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTc                                                                                                                            >  6:101104/1‑149 (MQ=255)         ccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTcc                                                                                                                           <  3:166631/147‑1 (MQ=255)         ccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTcc                                                                                                                           <  3:160511/147‑1 (MQ=255)         ccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTcc                                                                                                                           >  5:203591/3‑149 (MQ=255)           gCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTcccc                                                                                                                         >  7:67559/1‑148 (MQ=255)             cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTCGCGGCCCTTCG‑CGGGCACGCCCGCTCCcaca                                                                                                                       <  3:98354/149‑1 (MQ=255)                          gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTgag‑a                                                                                                          <  7:240111/148‑1 (MQ=255)                             gcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑Atctc                                                                                                      <  8:47118/149‑1 (MQ=255)                              ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATctct                                                                                                     >  6:351694/1‑149 (MQ=255)                                        ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGc                                                                                             <  7:101885/147‑1 (MQ=255)                                        ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                           >  6:325287/1‑149 (MQ=255)                                        ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCCCTCCCCCAGGTTTGCGGAG‑ATATCTGTGGGAGCgg                                                                                           >  3:398125/1‑149 (MQ=39)                                              tGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGt                                                                                       <  3:145528/147‑1 (MQ=255)                                                gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGccc                                                                                   <  5:351694/149‑1 (MQ=255)                                                gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGAGAT‑CTCAGCAAACCTGTGGGAGCGGGCGTGccc                                                                                   <  1:409401/149‑1 (MQ=18)                                                 cccgTACACAGTCCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCcg                                                                                  <  4:398125/146‑1 (MQ=255)                                                        acaGTGCCGGAGGCCAAGGCCTACATGGGTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGGTTGCGGAG‑ATCGCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                           >  1:236594/1‑149 (MQ=38)                                                        acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCACCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                           >  6:299521/1‑149 (MQ=255)                                                          aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCCCGCCCCCACGGTTGGTGAG‑ATCTCGGTGGGAGCGGGCGGGCCCGCGCAGGGcc                                                                         >  3:372531/1‑149 (MQ=255)                                                            tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCCCGCCCGCTCCCACACGGTTGCTGAG‑ATGTATGTGGGAGCGGGCGGGGCCGCGAAGGGCcgc                                                                       >  5:25783/1‑149 (MQ=38)                                                               cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                      >  1:296997/1‑147 (MQ=255)                                                                 gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCCCAACCACAGGTTTGGGGAG‑ATGAAGGTGGGAGCGGGCGGGGCCGCGCACGGCCGCGA‑AGc                                                                 >  8:262950/1‑149 (MQ=37)                                                                    gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACTGGTTTGCTGAG‑AGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGc                                                              <  6:203591/149‑1 (MQ=255)                                                                      catggCCTACTTGGGTTGGGGGGGTAGCGTGGTTCGGTGATCCTTTGACGGGGGCTGCTTTGCTGCCCTTCG‑CGGGCACGCCCGCTCCCACTGTTTTGCTTTC‑TTCTCTGTGGGTGGGGGGGTGCCCGCGAAGGGCCGCGA‑AGCGGccc                                                            <  4:372531/146‑1 (MQ=16)                                                                      cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGccc                                                            >  4:146545/1‑149 (MQ=255)                                                                       aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGccc                                                            >  5:322269/1‑148 (MQ=255)                                                                         ggCATACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGAAGCCCTTAG‑CGGGCACGCCCCCTCCCACACGTGTGCTGAG‑ATAGATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGACCCac                                                         >  6:129293/1‑148 (MQ=37)                                                                          gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAAt                                                        >  5:104058/1‑149 (MQ=255)                                                                               cATGGCTTGGGGGGGGATCGGCTTTCGTTTATCAATTGTCCGGGTGTGGTTTGCTGCCGTCGG‑GGGGCCCGCGCGCTCAAACCTGTTTTGTGAG‑ATCTTTGTGGTTGCGGGGGTCCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCtt                                                   <  6:25783/149‑1 (MQ=255)                                                                                aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                               >  3:306438/1‑90 (MQ=255)                                                                                aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                               <  4:306438/90‑1 (MQ=255)                                                                                       cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGt                                            <  5:190778/148‑1 (MQ=255)                                                                                           cgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTc                                       <  1:231196/149‑1 (MQ=255)                                                                                            gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                           <  6:33409/97‑1 (MQ=255)                                                                                            gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                           <  6:27160/97‑1 (MQ=255)                                                                                            gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                           >  5:33409/1‑97 (MQ=255)                                                                                            gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                           >  5:27160/1‑97 (MQ=255)                                                                                             ggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                     <  6:322269/149‑1 (MQ=255)                                                                                               aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCtgt                                                                                                   <  4:301978/86‑1 (MQ=255)                                                                                               aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCtgt                                                                                                   >  3:301978/1‑86 (MQ=255)                                                                                                 cgtcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGcc                                  <  5:325287/148‑1 (MQ=255)                                                                                                        cGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                           <  5:299521/148‑1 (MQ=255)                                                                                                          cTGATCCATTCCCCGGGGCTGCTTTGCAGCCCTGCG‑CGGGCACGCCCGCGCCCACAGGTTTGCTGAC‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGccc                                                            <  7:36947/113‑1 (MQ=21)                                                                                                           tGATCCATTGACCGGGGCTGCTTTGCATCCCTTCGTCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                            <  1:43393/67‑1 (MQ=18)                                                                                                           tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                            >  2:43393/1‑66 (MQ=255)                                                                                                           tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCCCTCCCACAGGTTTGGTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGt                       >  2:350899/1‑149 (MQ=255)                                                                                                             aTCCATTGACCGTGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACTGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                      <  5:129293/148‑1 (MQ=255)                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCCCAGGGTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGC‑AGCGGGCCCAATTGGTTCAACCGTGGTTCTTGCCGCGCAGCGCGTggga                   >  1:293914/1‑145 (MQ=37)                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGGCCCAATTGCTTAAACCGTGGTTGTTGCCGAGCAGCGCGTGGTa                   >  2:338670/1‑147 (MQ=255)                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACCCCCCCAACCACAGGTTTTCTGAG‑ATAAAAGTTGGAGCCGGCCTTCCCGCGAAGGGCCCCGA‑AGCGGCCCCAAATTATTAAAACGTGGTTAATTACGAGCAGCCCCTGGTaa                  >  2:147644/1‑148 (MQ=255)                                                                                                                      ccGGGGCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg            <  6:104058/149‑1 (MQ=255)                                                                                                                         gggCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                   <  7:235937/139‑1 (MQ=255)                                                                                                                         gggCTGCTTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                   >  8:235937/1‑139 (MQ=255)                                                                                                                               ctTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGACAGCGGCCCCAATTGCTTCAACCGTGGTTc                                       <  1:246910/114‑1 (MQ=255)                                                                                                                               ctTTGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGGTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGGGGTTc                                       >  2:246910/1‑113 (MQ=255)                                                                                                                                  tGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACCGGTTTCCTCTC‑ATCGCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  2:293914/147‑1 (MQ=38)                                                                                                                                  tGCAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTAGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  1:360912/147‑1 (MQ=255)                                                                                                                                    cAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                   <  3:301205/112‑1 (MQ=255)                                                                                                                                    cAGCCCTTCG‑CGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                   >  4:301205/1‑112 (MQ=255)                                                                                                                                                |                                                                                                                                        GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCG‑CAGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGA‑AGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352925‑353200

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