Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,248,398	C→T	6.5%	R558H (CGC→CAC)	PP_3724 ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,248,398	0	C	T	6.5%	183.0 / 2.7	61	R558H (CGC→CAC)	PP_3724	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (2/2);  total (29/32)
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.30e-01

GCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTC  >  NC_002947/4248260‑4248546                                                                                                                                           |                                                                                                                                                       gCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCAcc                                                                                                                                             >  2:69995/1‑150 (MQ=255)  cACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCa                                                                                                                                            <  2:388248/150‑1 (MQ=255)   aCCGATCACCAAGGCCAGGGCAAAGTTCTTCAGAGGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCaa                                                                                                                                           >  1:336752/1‑149 (MQ=255)    ccGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGt                                                                                                                                          <  2:256929/150‑1 (MQ=255)    ccGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGt                                                                                                                                          <  2:221764/150‑1 (MQ=255)       aTCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCa                                                                                                                                       >  1:519635/1‑150 (MQ=255)       aTCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCa                                                                                                                                       >  2:643875/1‑150 (MQ=255)           ccAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCtt                                                                                                                                   >  1:981854/1‑150 (MQ=255)           ccAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCtt                                                                                                                                   >  2:591978/1‑150 (MQ=255)                    gggCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGgcgc                                                                                                                          <  2:519635/150‑1 (MQ=255)                            ttgttTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATcg                                                                                                                  >  2:813796/1‑150 (MQ=255)                                tttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑Ca                                                                                                             <  1:4443/150‑1 (MQ=255)                                 ttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAg                                                                                                            <  2:688272/150‑1 (MQ=255)                                       cTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCaa                                                                                                     <  1:69995/150‑1 (MQ=255)                                        ttGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGt                                                                                                                                          <  1:1123910/114‑1 (MQ=255)                                          gATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCgg                                                                                                                        <  1:410084/130‑1 (MQ=255)                                          gATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCgg                                                                                                                        >  2:410084/1‑130 (MQ=255)                                          gATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGt‑t                                                                                                        >  2:544672/1‑144 (MQ=255)                                          gATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGt‑t                                                                                                        <  1:544672/144‑1 (MQ=255)                                                     tGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGAATCGCG‑CAGGT‑TCAACGCCACAGCACCAg                                                                                       >  2:975446/1‑150 (MQ=255)                                                            ggCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGAGCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTAGGt                                                                                >  1:1150411/1‑150 (MQ=255)                                                                 aGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTtgt                                                                           <  1:813796/150‑1 (MQ=255)                                                                  gggCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTtgtg                                                                          >  1:1194478/1‑150 (MQ=255)                                                                    gcggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCg                                                                        <  1:267728/150‑1 (MQ=255)                                                                      ggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGc                                                                      <  1:591978/150‑1 (MQ=255)                                                                        cTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAg                                                                    <  2:1194478/150‑1 (MQ=255)                                                                               aTGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcgcg                                                             <  1:218293/150‑1 (MQ=255)                                                                                   gCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCgg                                                                                                                        >  2:422841/1‑89 (MQ=255)                                                                                   gCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCgg                                                                                                                        <  1:422841/89‑1 (MQ=255)                                                                                   gCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGa                                                         >  1:1130496/1‑150 (MQ=255)                                                                                   gCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGa                                                         >  2:313949/1‑150 (MQ=255)                                                                                    cTCGAGCAGCGCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGt                                                                                <  2:51693/126‑1 (MQ=255)                                                                                    cTCGAGCAGCGCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGt                                                                                >  1:51693/1‑126 (MQ=255)                                                                                    cTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGAc                                                        <  1:431919/150‑1 (MQ=255)                                                                                          cagcaCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCtt                                                  >  2:1040298/1‑150 (MQ=255)                                                                                                  gATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTtg                                                                            <  1:198475/116‑1 (MQ=255)                                                                                                  gATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCa                                          >  1:362603/1‑150 (MQ=255)                                                                                                  gATCAAACGCAACCTGCGGGGCGCAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTtg                                                                            >  2:198475/1‑116 (MQ=255)                                                                                                       aaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCtt                                                                             >  1:1005694/1‑110 (MQ=255)                                                                                                       aaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCtt                                                                             <  2:1005694/110‑1 (MQ=255)                                                                                                       aaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                      <  2:525162/149‑1 (MQ=255)                                                                                                       aaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                      >  1:525162/1‑149 (MQ=255)                                                                                                        aCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTa                                             <  2:624954/141‑1 (MQ=255)                                                                                                        aCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTa                                             >  1:624954/1‑141 (MQ=255)                                                                                                            aaCCTGCGGGGCGTAGTCAGCGATGAAGCGGCTCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAAcc                                >  1:31710/1‑150 (MQ=255)                                                                                                                 gCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGAcgc                           >  2:305980/1‑150 (MQ=255)                                                                                                                  cGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGAcgcg                          <  2:362603/150‑1 (MQ=255)                                                                                                                     ggCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGa                              >  2:291103/1‑143 (MQ=255)                                                                                                                     ggCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGa                              <  1:291103/143‑1 (MQ=255)                                                                                                                        gTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGggcggc                    <  1:1121119/150‑1 (MQ=255)                                                                                                                           gTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGTCGGCGTc                 <  2:897529/150‑1 (MQ=255)                                                                                                                            tggaagctGCACTGGTGCAGATGCGCCAGTGCATCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt                >  2:721855/8‑150 (MQ=255)                                                                                                                            tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTACAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt                <  2:667881/150‑1 (MQ=255)                                                                                                                              aGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                            >  2:502764/1‑120 (MQ=255)                                                                                                                              aGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                            <  1:502764/120‑1 (MQ=255)                                                                                                                                  aTGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                        >  1:360046/1‑120 (MQ=255)                                                                                                                                  aTGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                        <  2:360046/120‑1 (MQ=255)                                                                                                                                   tGCACTGGTGCAGATGCGCCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGt                                                                                >  2:310432/1‑79 (MQ=37)                                                                                                                                   tGCACTGGTGCAGATGCGCCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGt                                                                                <  1:310432/79‑1 (MQ=37)                                                                                                                                   tGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACgg         <  2:981854/150‑1 (MQ=255)                                                                                                                                   tGCACTGGCGCATCTGCACCAGTGCAGATTCCAGTGCGCGAGC‑‑CGCGACATGTACCTACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACgg         <  2:336752/150‑1 (MQ=255)                                                                                                                                          gcgcATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTc  <  1:1040298/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                       GCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCG‑CAGGT‑TCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTC  >  NC_002947/4248260‑4248546

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