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%	217.4 / NA	64	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 (30/34);  total (30/34)

TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAT  >  NC_002947/352921‑353204                                                                                                                                                    |                                                                                                                                          ttCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                         >  2:49248/1‑149 (MQ=255)    tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGCTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGccc                                                                                                                                      >  5:1686/1‑147 (MQ=255)      gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                      <  2:26526/147‑1 (MQ=255)            cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCtg                                                                                                                              >  5:8451/1‑148 (MQ=255)             aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                              >  1:69581/1‑148 (MQ=255)             aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                              <  2:69581/148‑1 (MQ=255)                         tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                 <  8:233226/149‑1 (MQ=255)                           cTGGCGGCCAGCATCTTCCACTTAGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                               <  5:233505/149‑1 (MQ=255)                              gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                              <  3:236121/147‑1 (MQ=255)                                  ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATctct                                                                                                        <  6:77671/149‑1 (MQ=255)                                          ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                      >  4:116224/1‑143 (MQ=255)                                          ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                      <  3:116224/143‑1 (MQ=255)                                           tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGACCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                        >  1:274527/1‑125 (MQ=255)                                           tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                        <  2:274527/124‑1 (MQ=255)                                            caaCTTTGGCGATTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGc                                                                                                <  4:66070/145‑1 (MQ=255)                                              aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTggg                                                                                                   >  1:152219/1‑142 (MQ=255)                                              aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTggg                                                                                                   <  2:152219/142‑1 (MQ=255)                                               cTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCg                                                                                           >  4:59994/1‑149 (MQ=255)                                                tttGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACTTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGCGACCCCCGCTACCACAGGTGTGCTGACATCACGGTGGGAGCGGGcct                                                                                          >  5:232165/1‑147 (MQ=255)                                                  tGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                          <  5:304033/147‑1 (MQ=255)                                                      gAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGCTCTCTGTGGGAGCGGGCGTGCCcgc                                                                                    >  7:309467/1‑149 (MQ=255)                                                      gAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                    >  7:214666/1‑149 (MQ=255)                                                         tACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCCGAGATCTCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                 <  1:49248/149‑1 (MQ=255)                                                            acaGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                    <  8:75068/127‑1 (MQ=255)                                                            acaGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                    >  7:75068/1‑127 (MQ=255)                                                            acaGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                              >  5:189699/1‑149 (MQ=255)                                                            acaGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                              >  1:80472/1‑149 (MQ=255)                                                                tGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                            <  2:57773/131‑1 (MQ=255)                                                                tGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                            >  1:57773/1‑131 (MQ=255)                                                                tGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgc                                                                          >  8:301427/1‑149 (MQ=255)                                                                ggCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATGGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgc                                                                          <  3:34992/148‑1 (MQ=255)                                                                  ccGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                   <  7:89953/138‑1 (MQ=255)                                                                  ccGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                   >  8:89953/1‑138 (MQ=255)                                                                   cGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCGGAGATGTCTGTGGGAGCGGGCGGGCCCGCGAAGGGCcgcg                                                                         >  4:50454/1‑147 (MQ=255)                                                                     gAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                     >  7:34445/1‑149 (MQ=255)                                                                         ccAAGG‑CCTACATGGCTTGGCGCGGGTTCGTCGTTCGCTTATCCATTGACGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACTGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                           <  4:296165/139‑1 (MQ=38)                                                                         ccAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                           >  3:296165/1‑139 (MQ=255)                                                                          cAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                >  6:399869/1‑149 (MQ=255)                                                                                ccTACCTGGCTTCGCGCGGGATCGTCGTGCGCTGATCCATTCCCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                           <  6:1686/149‑1 (MQ=38)                                                                                 cTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACACCCCCAACCAAACAATTTCGGGGATGACAGTGGGAACGGGCGGGCCCCCGAAAGGGCGCGCACCGGGCCCAAATg                                                          >  4:224108/1‑149 (MQ=255)                                                                                    cATGGCTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGc                                                         <  3:50454/147‑1 (MQ=255)                                                                                        gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTTAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaa                                                    <  8:214666/148‑1 (MQ=255)                                                                                               gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCg                                                                    <  4:184452/125‑1 (MQ=255)                                                                                               gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCg                                                                    >  3:184452/1‑125 (MQ=255)                                                                                                    aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                       <  7:301427/149‑1 (MQ=255)                                                                                                    aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                       >  5:403432/1‑149 (MQ=255)                                                                                                       gtcttTCGGTGTTCGATTGACGGGGGGTGCTTTGCAGCCTTTCGCGGGCAGGCCCGCTCCCACTGGTTTGCTTCCATCTCTGGGGGAGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                    <  6:232165/149‑1 (MQ=255)                                                                                                             cGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcg                              >  2:303347/1‑149 (MQ=255)                                                                                                              gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcgc                             >  2:319969/1‑149 (MQ=255)                                                                                                               cTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                            <  1:319969/149‑1 (MQ=255)                                                                                                                     cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                       >  3:330815/1‑148 (MQ=255)                                                                                                                      aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                      >  3:187498/1‑68 (MQ=255)                                                                                                                      aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                      <  4:187498/68‑1 (MQ=255)                                                                                                                      aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                     >  4:135728/1‑149 (MQ=255)                                                                                                                           ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTTAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                <  8:255967/149‑1 (MQ=255)                                                                                                                           ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                <  3:59994/149‑1 (MQ=255)                                                                                                                           ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                <  2:80472/149‑1 (MQ=255)                                                                                                                             ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt              <  8:34445/149‑1 (MQ=255)                                                                                                                             ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt              <  6:93679/149‑1 (MQ=255)                                                                                                                              gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                <  3:17757/146‑1 (MQ=255)                                                                                                                              gggCGGCTTGCCTGCCCTGCGCGGCCACGCCCGCTCCCACTGTTTGCCTCCTATCTCTGTGGGCGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTa             <  7:346146/149‑1 (MQ=25)                                                                                                                                        gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTCATATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt    <  3:224108/148‑1 (MQ=255)                                                                                                                                         cAGCCCTTCGCGGGCACGCCCGCTACCACAGGTTTGCGGAGATCTATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGGTTCAACCGTGGGTCTTGCCGAGCAGCGCGTGGTTAAGCACGGTATCCCCGAGGAt  >  8:182809/1‑149 (MQ=37)                                                                                                                                               ttCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                    >  2:141004/1‑45 (MQ=255)                                                                                                                                               ttCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                    <  1:141004/45‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                          TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGG‑CCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAT  >  NC_002947/352921‑353204

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