Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,468,261	G→C	10.2%	intergenic (‑25/+98)	PP_2161 ← / ← PP_2162	hypothetical protein/phosphonoacetaldehyde phosphonohydrolase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,468,261	0	G	C	10.2%	166.1 / 9.7	59	intergenic (‑25/+98)	PP_2161/PP_2162	hypothetical protein/phosphonoacetaldehyde phosphonohydrolase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (30/23);  new base C (3/3);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACCAGGCAG  >  NC_002947/2468117‑2468407                                                                                                                                                 |                                                                                                                                                    cgcCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCaa                                                                                                                                                  <  7:912/148‑1 (MQ=255)   ccTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGc                                                                                                                                                >  5:217644/1‑148 (MQ=255)     tCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAt                                                                                                                                             >  8:276716/1‑149 (MQ=255)       gATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCATGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAt                                                                                                                                             <  5:270832/147‑1 (MQ=255)             gcgcgcCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCTGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCATTGatgg                                                                                                                                       <  7:276716/147‑4 (MQ=255)             gcgcgcCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACctct                                                                                                                                     >  3:139389/1‑149 (MQ=255)                      ttCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGt                                                                                                                              <  5:316402/147‑1 (MQ=255)                            cAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTTTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAgg                                                                                                                      <  7:24230/149‑1 (MQ=255)                              gggCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAttg                                                                                                                                           >  7:149593/1‑124 (MQ=255)                              gggCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGCCAAGCCAttg                                                                                                                                           <  8:149593/126‑3 (MQ=255)                                gcggcCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCt                                                                                                                  >  4:272392/1‑149 (MQ=255)                                 cggcCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTg                                                                                                                 <  3:305581/149‑1 (MQ=255)                                  ggcCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTgg                                                                                                                <  3:272392/149‑1 (MQ=255)                                    ccAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGcc                                                                                                              >  2:18402/1‑149 (MQ=255)                                    ccAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGcc                                                                                                              >  4:275534/1‑149 (MQ=255)                                      aGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTg                                                                                                            >  1:335466/1‑149 (MQ=255)                                        ctctTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTc                                                                                                          <  1:18402/149‑1 (MQ=255)                                        ctctTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTc                                                                                                          <  5:175089/149‑1 (MQ=255)                                         tcttctTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGt                                                                                                           >  6:308757/1‑147 (MQ=255)                                         tcttctTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGt                                                                                                           <  5:308757/147‑1 (MQ=255)                                          cttcttCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTATCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTcat                                                                                                        >  6:347557/1‑147 (MQ=255)                                             cttcCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGtcatca                                                                                                      <  5:347557/148‑1 (MQ=255)                                             cttcCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTc                                                                                                       <  6:35473/147‑1 (MQ=255)                                             cttcCGAAAGCGGGACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCa                                                                                                      <  2:122094/148‑1 (MQ=255)                                                 cGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAgg                                                                                                    <  8:44619/146‑1 (MQ=255)                                                 cGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAgg                                                                                                    >  7:44619/1‑146 (MQ=255)                                                       cGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTAc                                                                                           <  6:217644/149‑1 (MQ=255)                                                         gTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACct                                                                                         >  1:339078/1‑149 (MQ=255)                                                                  ggTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGAt                                                                               <  2:219159/149‑1 (MQ=255)                                                                            cTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATc                                                                     >  7:317527/1‑149 (MQ=255)                                                                             tCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATcc                                                                    >  1:155845/1‑149 (MQ=255)                                                                               gcgcAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCcgc                                                                  >  7:6264/1‑149 (MQ=255)                                                                                  cAAGCTTTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGc                                                                           <  3:206418/137‑1 (MQ=255)                                                                                  cAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGc                                                                           >  4:206418/1‑137 (MQ=255)                                                                                    aGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTc                                                               >  2:131020/1‑147 (MQ=255)                                                                                        gTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACg                                                         >  4:151770/1‑149 (MQ=255)                                                                                         ttCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACgg                                                        <  3:275534/149‑1 (MQ=255)                                                                                                 ctctTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGt                                                                >  5:200954/1‑133 (MQ=255)                                                                                                 ctctTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGt                                                                <  6:200954/133‑1 (MQ=255)                                                                                                      gCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAggatc                                           >  7:208412/1‑146 (MQ=255)                                                                                                         aTTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTGAGCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGtttt                                         <  1:65545/149‑1 (MQ=255)                                                                                                                     gcaTTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATCTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGcc                            <  5:128219/149‑1 (MQ=255)                                                                                                                       aTTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGCCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGATGTATCAGGGTTTTTCACCCTTGCTCCgg                          >  5:32350/1‑149 (MQ=255)                                                                                                                           gCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCt                      >  8:309214/1‑149 (MQ=255)                                                                                                                                 ctgcATCTATGGCTTGGCAAGCGATGGAACTGTGCGAGGTGGCATAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTTGGCg                 <  6:32350/145‑1 (MQ=255)                                                                                                                                    cATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGt             >  7:32161/1‑149 (MQ=255)                                                                                                                                    cATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGt             >  4:28877/1‑149 (MQ=255)                                                                                                                                       ccaTGGCTTGCCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTcca          >  3:239146/3‑149 (MQ=255)                                                                                                                                       cAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTcca          >  8:175731/1‑149 (MQ=255)                                                                                                                                        aaTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTccac         >  3:288354/1‑149 (MQ=255)                                                                                                                                        aaTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTccac         >  8:260645/1‑149 (MQ=255)                                                                                                                                        aaTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGAAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTccac         <  2:339078/149‑1 (MQ=255)                                                                                                                                         aTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTgc                               <  6:294570/126‑1 (MQ=255)                                                                                                                                         aTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTgc                               >  5:294570/1‑126 (MQ=255)                                                                                                                                         aTGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTccac         >  7:13500/1‑148 (MQ=255)                                                                                                                                            gCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcagg     >  6:150189/1‑149 (MQ=255)                                                                                                                                            gCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcagg     >  4:129133/1‑149 (MQ=255)                                                                                                                                             cTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcaggc    <  4:139389/149‑1 (MQ=255)                                                                                                                                             cTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcagg     >  3:149530/1‑148 (MQ=255)                                                                                                                                                ggCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACcaggcag  <  2:155845/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CGCCTCGATCTGGCGCGCCTCTTCGTGCAGGGCGGCCAGCTCTTCTTCCGAAAGCGGTACGTTGCGGTCCAGTTGCTCGCGCAAGCTGTTCAGGCGCTCTTGCAATTCGCGGGCAGGCATTGGCTATTCTCCATCAATGGCTTGGCAAGCCATGGACCTCTGCGACGTGGCAAAGGTTCTGGCCTGTCTTCAGGGTAATCTACCTGCTG‑GCGATTTGCATGATCCGCGTCAAGACGGCTGTATCAGGGTTTTTCACCCTTGCGCCGGCGCTGGGCGATGTCCACCAGGCAG  >  NC_002947/2468117‑2468407

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