Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,253,433	G→C	6.4%	intergenic (‑79/+28)	PP_1096 ← / ← metG	electron transport complex protein/methionine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,253,433	0	G	C	6.4%	201.6 / 4.0	63	intergenic (‑79/+28)	PP_1096/metG	electron transport complex protein/methionine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base C (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

GCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCGGG  >  NC_002947/1253288‑1253579                                                                                                                                                  |                                                                                                                                                    gCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCg                                                                                                                                                  >  3:172488/1‑149 (MQ=255)    gTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGAGTAATTGTGGGAGCGGGCAGGCCGcga                                                                                                                                               >  8:284230/1‑147 (MQ=255)     tCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTcc                                                                                                                                              >  1:74050/1‑149 (MQ=255)      cGGAGCAGGCACAGGTGATTGATCAGTGCTGAGCTGATCAGGACCAGGATGTAGTCATTCATGAGAGGACGATAGCCGGCAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGCCTGCCCGCTcc                                                                                                                                              <  5:187964/148‑1 (MQ=255)              gCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGc                                                                                                                                     <  2:46386/149‑1 (MQ=255)                acaGGTGATTGATCAGTGCGGCGCTGATCATGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGccc                                                                                                                                   >  6:5740/1‑149 (MQ=255)                acaGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGccc                                                                                                                                   >  4:253288/1‑149 (MQ=255)                acaGGTGATTGATCAGTGCGGCCCTGATCAGGACCAGGATGTAGTCATTCATGATATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGccc                                                                                                                                   <  3:14685/149‑1 (MQ=255)                        ttgatCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCg                                                                                                                           <  2:173454/149‑1 (MQ=255)                        ttgatCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCg                                                                                                                           <  7:303186/149‑1 (MQ=255)                         tgatCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCgg                                                                                                                          >  6:332645/1‑149 (MQ=255)                           atCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGt                                                                                                                         >  1:218790/1‑148 (MQ=255)                               gTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACtt                                                                                                                    <  2:218790/149‑1 (MQ=255)                                  cggcgCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCa                                                                                                                                >  4:128831/1‑134 (MQ=255)                                  cggcgCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCa                                                                                                                                <  3:128831/134‑1 (MQ=255)                                  cggcgCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGa                                                                                                                  <  3:115746/148‑1 (MQ=255)                                    gcgcTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATgc                                                                                                               <  1:194193/149‑1 (MQ=255)                                         gATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctgg                                                                                                          >  5:83923/1‑149 (MQ=255)                                          aTCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctggc                                                                                                         >  2:135200/1‑149 (MQ=255)                                          aTCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctg                                                                                                           <  4:357086/147‑1 (MQ=255)                                             aGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGctggcctg                                                                                                      >  3:35879/1‑149 (MQ=255)                                               gACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCg                                                                                                                           >  6:122419/1‑126 (MQ=255)                                               gACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGGTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCg                                                                                                                           <  5:122419/126‑1 (MQ=255)                                                  cAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTg                                                                                                 >  3:129498/1‑149 (MQ=255)                                                      aTGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgcgc                                                                                             <  7:251315/149‑1 (MQ=255)                                                      aTGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgcgc                                                                                             <  4:129498/149‑1 (MQ=255)                                                      aTGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgcgc                                                                                             <  5:238369/149‑1 (MQ=255)                                                           gTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTg                                                                                        >  5:260592/1‑149 (MQ=255)                                                              attcatGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑cc                                                                                    <  6:260592/149‑1 (MQ=255)                                                               ttcatGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCg                                                                                   >  8:274243/1‑149 (MQ=255)                                                                tcatGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCgg                                                                                  >  3:189306/1‑149 (MQ=255)                                                                  atGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGTAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGc                                                                                >  7:126847/1‑149 (MQ=255)                                                                  atGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCggg                                                                                 <  5:5740/148‑1 (MQ=255)                                                                  atGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGc                                                                                >  7:81616/1‑149 (MQ=255)                                                                       aTGACGATAGCCGGGAAGCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGGTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTcagc                                                                           <  6:346731/149‑1 (MQ=255)                                                                            gATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTa                                                                      <  4:172488/149‑1 (MQ=255)                                                                              tAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGa                                                                    <  3:253288/149‑1 (MQ=255)                                                                                     gAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATttcttc                                                             <  8:243233/149‑1 (MQ=255)                                                                                          ggggCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGc                                                        <  6:83923/149‑1 (MQ=255)                                                                                          ggggCGATTGTCCGGGTTTGACATGAAGGGGACTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGATACTTGAAGCGCTGGCCTGGCTTGGCCACGCTGT‑CCGGGCACAGCAGGTAGATTTCTTCACCGc                                                        >  4:355766/1‑149 (MQ=255)                                                                                               gATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAggg                                                   >  8:332355/1‑149 (MQ=255)                                                                                               gATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTACCACAATGGCCCGCAGGGCGGTTACTTGATTCGCTGCCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAggg                                                   >  5:109820/1‑149 (MQ=255)                                                                                                    tCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCgg                                               >  2:297533/1‑148 (MQ=255)                                                                                                    tCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCgg                                               >  3:200897/1‑148 (MQ=255)                                                                                                          tttGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCAt                                        <  2:74050/149‑1 (MQ=255)                                                                                                             gACATGAAGGGGCCTGCGTAATTGTGGGAGCGTGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGa                                      <  1:292589/148‑1 (MQ=255)                                                                                                             gACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGa                                      <  7:274243/148‑1 (MQ=255)                                                                                                                 tGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGTCCCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCAt                                  <  1:294305/149‑1 (MQ=255)                                                                                                                   aaGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGc                                <  6:143972/148‑1 (MQ=255)                                                                                                                   aaGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGc                                <  6:137890/148‑1 (MQ=255)                                                                                                                      gggCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTc                            >  6:221291/1‑149 (MQ=255)                                                                                                                      gggCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAACGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTc                            <  5:332645/149‑1 (MQ=255)                                                                                                                        gCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCg                           >  1:191869/1‑148 (MQ=255)                                                                                                                          cTGCGTAATTGTGGGAGCTGGCTGGAGGCTCCCACAATGGCCCGGAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCTGac                        <  8:126847/149‑1 (MQ=255)                                                                                                                             cGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACg                     >  7:63221/1‑149 (MQ=255)                                                                                                                                     tgGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAgcg             >  5:304854/1‑149 (MQ=255)                                                                                                                                      gggAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGCCCATGCCTTCGGACACGCCAAAgcg             >  6:369471/1‑148 (MQ=255)                                                                                                                                      gggAGCGGCCTGCCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAgcg             >  8:199247/1‑148 (MQ=255)                                                                                                                                       ggAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGCCCCTGCCTTCGGCCACGCCAAAGCGCa           >  2:124043/1‑149 (MQ=255)                                                                                                                                        gAGCGGCCTGCCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgc                                                                                               >  7:330376/1‑65 (MQ=37)                                                                                                                                        gAGCGGCCTGCCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGgc                                                                                               <  8:330376/65‑1 (MQ=37)                                                                                                                                                aGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCCGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCggg  >  6:258545/1‑149 (MQ=255)                                                                                                                                                aGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCggg  <  5:221291/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCTGTCGGAGCAGGCACAGGTGATTGATCAGTGCGGCGCTGATCAGGACCAGGATGTAGTCATTCATGAGATGACGATAGCCGGGAATCGGGGCGATTGTCCGGGTTTGACATGAAGGGGCCTGCGTAATTGTGGGAGCGGGCAGGCCGCTCCCACAATGGCCCGCAGGGCGGTTACTTGATGCGCTGGCCTGGCTTGGCGCCGCTGT‑CCGGGCTCAGCAGGTAGATTTCTTCACCGCCAGGGCCGGCAGCCATGACCATGCCTTCGGACACGCCAAAGCGCATCTTCCGGG  >  NC_002947/1253288‑1253579

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