Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	33,975	C→T	6.2%	R138C (CGT→TGT)	czrSA →	metal‑sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	33,975	0	C	T	6.2%	193.7 / 2.0	64	R138C (CGT→TGT)	czrSA	metal‑sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (2/2);  total (33/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 = 4.71e-01

CTCACAGAATCAGACTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTGGCCATTGGCTGGTGGACT  >  NC_002947/33830‑34118                                                                                                                                                  |                                                                                                                                                 cTCACAGAATCAGACTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGtt                                                                                                                                               >  5:273613/1‑149 (MQ=255)    acaGAATCAGACTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTgg                                                                                                                                            <  6:273613/149‑1 (MQ=255)             gACTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGaa                                                                                                                                   >  4:159138/1‑149 (MQ=255)              aCTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACTAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAAc                                                                                                                                  <  2:61781/149‑1 (MQ=255)              aCTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGaa                                                                                                                                   >  5:247073/1‑148 (MQ=255)                  cccTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACagag                                                                                                                              <  1:185967/149‑1 (MQ=255)                  cccTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCACAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACagag                                                                                                                              <  4:273762/149‑1 (MQ=255)                     tGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGt                                                                                                                            >  3:6802/1‑148 (MQ=255)                       ctctCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTc                                                                                                                           <  6:247073/147‑1 (MQ=255)                       ctctCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCa                                                                                                                          <  4:6802/148‑1 (MQ=255)                       cgctcGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTc                                                                                                                           <  4:77318/145‑1 (MQ=255)                            gATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGc                                                                                                                        <  1:294733/145‑1 (MQ=255)                            gATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGCACCGTCCTCAGCTGAGTGATTCTCGAATAATGCGTTTAGTCGATGGAACAGACGTCAGc                                                                                                                        >  2:294733/1‑145 (MQ=255)                               tGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCtat                                                                                                                  <  3:15717/148‑1 (MQ=255)                                   tCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGAc                                                                                                             >  4:283056/1‑149 (MQ=255)                                    cGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGAc                                                                                                             >  8:253125/1‑148 (MQ=255)                                      aaaCAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAt                                                                                                          <  3:159138/149‑1 (MQ=255)                                       aaCAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑g                                                                                                        >  2:246322/1‑149 (MQ=255)                                              cAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGt                                                                                                 >  2:116022/1‑149 (MQ=255)                                               aaTCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGTACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTgg                                                                                                                                            >  1:227842/1‑106 (MQ=255)                                               aaTCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGTACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTgg                                                                                                                                            <  2:227842/106‑1 (MQ=255)                                                           cacaAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAgg                                                                                     <  1:116022/148‑1 (MQ=255)                                                             caAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATc                                                                                             <  3:41947/138‑1 (MQ=255)                                                             caAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATc                                                                                             >  4:41947/1‑138 (MQ=255)                                                              aaGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCgct                                                                                 >  8:303348/1‑149 (MQ=255)                                                                 tttCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCgctgct                                                                              >  6:76638/1‑149 (MQ=255)                                                                   tCAGCCGGCGGCCCAGCCCAAAGAGCACGAAAGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GCACCAGTCATCCAGTGAGGCGCTGCTc                                                                             >  1:197073/1‑148 (MQ=255)                                                                     aGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGGGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGAt                                                                          >  5:331798/1‑149 (MQ=255)                                                                      gCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATg                                                                         >  8:105159/1‑149 (MQ=255)                                                                      gCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATg                                                                         <  7:303348/149‑1 (MQ=255)                                                                       ccGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATc                                                                                             >  3:148447/1‑128 (MQ=255)                                                                       ccGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATc                                                                                             <  4:148447/128‑1 (MQ=255)                                                                         ggcggcCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGt                                                                                         <  2:97033/130‑1 (MQ=255)                                                                         ggcggcCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGt                                                                                         >  1:97033/1‑130 (MQ=255)                                                                         ggcggcCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGt                                                                      <  3:283056/149‑1 (MQ=255)                                                                           cggcCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTg                                                                    <  5:76638/149‑1 (MQ=255)                                                                                  gcccaAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGc                                                                                                                                          >  5:319848/1‑73 (MQ=255)                                                                                  gcccaAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGc                                                                                                                                          <  6:319848/73‑1 (MQ=255)                                                                                     caAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGTGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGAc                                                                                                     >  4:119795/1‑106 (MQ=255)                                                                                     caAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGTGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGAc                                                                                                     <  3:119795/106‑1 (MQ=255)                                                                                     caAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGTGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGAc                                                                                                     >  6:212256/1‑106 (MQ=255)                                                                                     caAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGTGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGAc                                                                                                     <  5:212256/106‑1 (MQ=255)                                                                                      aaaGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGc                                                         >  6:319984/1‑149 (MQ=255)                                                                                            cACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTtgat                                                   >  6:88357/1‑149 (MQ=255)                                                                                             aCGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCATCGGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTtgat                                                   <  1:246322/149‑1 (MQ=255)                                                                                               gAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTgatgat                                                >  5:187733/1‑149 (MQ=255)                                                                                                    gCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTccc                                           <  6:187733/149‑1 (MQ=255)                                                                                                     cGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCa                                                                                                           <  8:271605/85‑1 (MQ=255)                                                                                                     cGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCa                                                                                                           >  7:271605/1‑85 (MQ=255)                                                                                                     cGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTcccc                                          >  5:86209/1‑149 (MQ=255)                                                                                                       tGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCg                                        >  2:27154/1‑149 (MQ=255)                                                                                                               aGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTc                                <  7:105159/149‑1 (MQ=255)                                                                                                                   ccGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCAc                            >  2:310195/1‑149 (MQ=255)                                                                                                                         tctGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCt                                                                                                                    >  1:581/1‑56 (MQ=255)                                                                                                                         tctGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCt                                                                                                                    <  2:581/56‑1 (MQ=255)                                                                                                                         tctGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCt                      >  8:78872/1‑149 (MQ=255)                                                                                                                         tctGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCt                      >  1:172487/1‑149 (MQ=255)                                                                                                                             ctgAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTgg                   <  1:310195/148‑1 (MQ=255)                                                                                                                                       ctcGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTGGCCATTGGCtgg        >  5:36036/1‑149 (MQ=255)                                                                                                                                       ctcGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTGGCCATTGGCt          <  2:172487/147‑1 (MQ=255)                                                                                                                                          gAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTGGCCATTGGCtggtg      <  5:319984/148‑1 (MQ=255)                                                                                                                                             tCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTgg                                                           >  4:111261/1‑92 (MQ=255)                                                                                                                                             tCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTgg                                                           <  3:111261/92‑1 (MQ=255)                                                                                                                                             tCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTGGCCATTGGCTGGTGGACt  <  7:253125/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CTCACAGAATCAGACTCCCCTGCTCTCGATTGGCTCGAAACAGGTCAATCTGGAGGTGCACAAGTTTCAGCCGGCGGCCCAGCCCAAAGAGCACGAATGGCGTGATACGCAGGACCGTCCTCTGCTGAGTGCTTCTCGAATCATGCGTTTGGGCGATGGAACAGAGGTCAGCGTCTATATGACCAT‑GGACCAGTCATCCAGTGAGGCGCTGCTCGATGCGTTGCTGACCTGGGCCTTGATGATGTCCCCCGTCATCCTCGCACTGATCTTGGCCATTGGCTGGTGGACT  >  NC_002947/33830‑34118

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