Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,174,709	A→C	8.7%	intergenic (‑71/+48)	glmU ← / ← atpC	bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,174,709	0	A	C	8.7%	191.9 / 7.3	67	intergenic (‑71/+48)	glmU/atpC	bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base C (3/3);  total (34/33)
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

GGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGA  >  NC_002947/6174562‑6174856                                                                                                                                                    |                                                                                                                                                     gggTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAg                                                                                                                                                     <  6:210625/149‑1 (MQ=255)       aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAc                                                                                                                                                >  7:154044/1‑148 (MQ=255)        acacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTAccc                                                                                                                                              >  1:299725/1‑149 (MQ=255)         cacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCt                                                                                                                                             <  2:299725/149‑1 (MQ=255)         cacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                      <  6:28143/140‑1 (MQ=255)         cacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                      >  5:28143/1‑140 (MQ=255)          acTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCtt                                                                                                                                            >  7:134045/1‑149 (MQ=255)              tGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                         <  6:720/148‑1 (MQ=255)                ggCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACt                                                                                                                                       >  2:280214/1‑148 (MQ=255)                 gCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCTTTTTCttt                                                                                                                                     <  2:310113/149‑3 (MQ=255)                      gcCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                                  <  3:139777/147‑1 (MQ=255)                      gcCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                                  >  4:139777/1‑147 (MQ=255)                                    ggTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCa                                                                                                                  <  8:134045/149‑1 (MQ=255)                                        ccTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCg                                                                                                              >  7:183214/1‑149 (MQ=255)                                         cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                             >  3:18295/1‑149 (MQ=255)                                         cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                             >  8:172207/1‑149 (MQ=255)                                                 gcggcgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAAAAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTAATCAATCGCAACAGAAGTTATAACCCGGGGCCGAt                                                                                                      >  3:256158/1‑148 (MQ=255)                                                   ggcgAGAATAACTATATCTAGTGACATTACTGCCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACATAAGTTATCAGCCggg                                                                                                            <  4:295704/140‑1 (MQ=255)                                                   ggcgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCggg                                                                                                            >  3:295704/1‑140 (MQ=255)                                                    gcgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACt                                                                                                  >  7:86730/1‑149 (MQ=255)                                                     cgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGGCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACtt                                                                                                 <  2:86621/149‑1 (MQ=255)                                                             aCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                         >  5:285252/1‑149 (MQ=255)                                                             aCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                         <  4:18295/149‑1 (MQ=255)                                                                 tatCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                                  <  5:140970/104‑1 (MQ=255)                                                                 tatCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                                  >  6:140970/1‑104 (MQ=255)                                                                 tatCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAt                                                                                     <  7:172207/149‑1 (MQ=255)                                                                   tCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAttt                                                                                   <  8:117161/149‑1 (MQ=255)                                                                          aCATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAGGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                            >  6:122992/1‑149 (MQ=255)                                                                          aCATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGa                                                                             <  8:154044/148‑1 (MQ=255)                                                                           cATTACTGGCTACCATTCCTGGGCGGGCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                            <  4:218384/148‑1 (MQ=255)                                                                           cATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACg                                                                           <  8:86730/149‑1 (MQ=255)                                                                           cATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTAAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACg                                                                           >  7:213575/1‑149 (MQ=255)                                                                                 tGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATc                                                                                                                   <  7:35023/103‑1 (MQ=255)                                                                                 tGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATc                                                                                                                   >  8:35023/1‑103 (MQ=255)                                                                                     tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCt                                                                 >  2:220889/1‑149 (MQ=255)                                                                                       ccAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTgg                                                               <  1:280214/149‑1 (MQ=255)                                                                                        cAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTggg                                                              >  7:162984/1‑149 (MQ=255)                                                                                         aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATc                                                                                                                   >  4:223942/1‑95 (MQ=255)                                                                                         aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGc                                                             >  5:285497/1‑149 (MQ=255)                                                                                          aTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGc                                                                                 >  4:204808/1‑128 (MQ=255)                                                                                          aTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGc                                                                                 <  3:204808/128‑1 (MQ=255)                                                                                             cTGGGCGGTCAGTGACAGACCGGGGGAGAAAACAGAAAAAGAACAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                         <  8:178143/149‑1 (MQ=255)                                                                                               gggCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑c                                                      <  7:200778/149‑1 (MQ=255)                                                                                                  cGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTc                                                   <  6:285497/149‑1 (MQ=255)                                                                                                   ggTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGc                                                          <  1:254896/142‑1 (MQ=255)                                                                                                   ggTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGc                                                          >  2:254896/1‑142 (MQ=255)                                                                                                   ggTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCg                                                  >  5:91698/1‑149 (MQ=255)                                                                                                     tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                         >  4:43438/1‑141 (MQ=255)                                                                                                     tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                         <  3:43438/141‑1 (MQ=255)                                                                                                     tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGc                                                <  8:22067/149‑1 (MQ=255)                                                                                                         tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGa                                            >  2:163533/1‑149 (MQ=255)                                                                                                           aaaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                         <  7:107794/55‑1 (MQ=255)                                                                                                           aaaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                         >  8:107794/1‑55 (MQ=255)                                                                                                                           aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCg                          >  1:261063/1‑149 (MQ=255)                                                                                                                              agaaaaagaaaaAGGGTAGCCGAGGCTCCCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGTGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAg                       >  7:281011/1‑149 (MQ=255)                                                                                                                              agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAg                       >  4:186035/1‑149 (MQ=255)                                                                                                                                      taaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGt                      <  1:99726/141‑1 (MQ=255)                                                                                                                                      taaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGt                      >  2:99726/2‑142 (MQ=255)                                                                                                                                       aaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTAATTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTCCGG‑CCTAGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTcc              >  3:164744/1‑149 (MQ=255)                                                                                                                                       aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAGGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTcc              <  5:122992/149‑1 (MQ=255)                                                                                                                                          agGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc            <  1:220889/149‑1 (MQ=255)                                                                                                                                           gggTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCgcgc          >  3:299943/1‑149 (MQ=255)                                                                                                                                            ggTAGCCGAGGCTACCCTTTTACTAAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCACCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGcc         >  3:263097/1‑149 (MQ=255)                                                                                                                                                gCCGTGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGTGGGCGATTACTTGCCTTTGCGCATTTGCTGGACGGTAGGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGATGTCCGCGCCTTTc     <  8:183214/149‑1 (MQ=255)                                                                                                                                                gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTc     <  4:256158/149‑1 (MQ=255)                                                                                                                                                gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTc     <  4:299943/149‑1 (MQ=255)                                                                                                                                                gCCGAGGCTACCCTATTACTCAATCGCAACAGAAGTTATCAGCAGGGGCCGAATACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTCCGG‑CCTCGGCCAGACGTGCGGAGGCGCCGCCGTAGTCGAAGTCCGCGCCTTTc     >  5:254003/1‑149 (MQ=255)                                                                                                                                                   gAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGa  <  2:245056/149‑1 (MQ=255)                                                                                                                                                   gAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGa  <  8:162984/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGG‑CCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGA  >  NC_002947/6174562‑6174856

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