Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,508,413	A→T	8.9%	V123E (GTA→GAA)	cysG ←	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,508,413	0	A	T	8.9%	188.6 / 5.8	65	V123E (GTA→GAA)	cysG	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base T (3/3);  total (30/35)
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

ACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAC  >  NC_002947/4508266‑4508555                                                                                                                                                    |                                                                                                                                                aCGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTAc                                                                                                                                                <  2:194110/149‑1 (MQ=255)   gTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTAc                                                                                                                                                >  4:87946/1‑147 (MQ=255)    tCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑g                                                                                                                                            <  7:47242/149‑1 (MQ=255)        ggTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGTAATCGCTGTTTCCA‑gcggc                                                                                                                                        <  3:124653/149‑1 (MQ=255)        ggTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑gcggc                                                                                                                                        <  4:211435/149‑1 (MQ=255)             aGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACg                                                                                                                                   <  8:95439/149‑1 (MQ=255)                      aCCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGAt                                                                                                                          >  7:267325/1‑149 (MQ=255)                        cTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCg                                                                                                                        >  1:124953/1‑149 (MQ=255)                            tgCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCgg                                                                                                                    <  8:273968/149‑1 (MQ=255)                             gCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCATGCGGCGAACGGTCCACGATCGCCg                                                                                                                     <  7:228141/148‑1 (MQ=255)                              ccGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGa                                                                                                                  <  2:124953/149‑1 (MQ=255)                                 gAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCggg                                                                                                                   >  2:109749/1‑145 (MQ=255)                                      gTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCAcc                                                                                                          >  7:17703/1‑149 (MQ=255)                                      gTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCAcc                                                                                                          >  6:54408/1‑149 (MQ=255)                                      gTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCAcc                                                                                                          >  6:108521/1‑149 (MQ=255)                                       tGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCg                                                                                                         >  2:369513/1‑149 (MQ=255)                                          agccagccCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑gc                                                                                                                                           <  7:145121/112‑1 (MQ=255)                                          agccagccCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑gc                                                                                                                                           >  8:145121/1‑112 (MQ=255)                                                 ccGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGc                                                                                               >  1:233857/1‑149 (MQ=255)                                                     ccAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGgc                                                                                           >  3:7471/1‑149 (MQ=255)                                                        gCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGa                                                                                                               <  2:345949/126‑1 (MQ=255)                                                        gCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCAGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGa                                                                                                               >  1:345949/1‑126 (MQ=255)                                                          tCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGAt                                                                                                                          <  5:56873/113‑1 (MQ=38)                                                          tCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGAt                                                                                                                          >  6:56873/1‑113 (MQ=38)                                                                     cGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGa                                                                           >  2:268258/1‑149 (MQ=255)                                                                       aaGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGAcc                                                                         >  4:360284/1‑149 (MQ=255)                                                                          ggATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGc                                                                      <  3:26333/149‑1 (MQ=255)                                                                                aCGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCt                                                                 >  5:144215/1‑148 (MQ=255)                                                                                aCGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCt                                                                 >  6:266710/1‑148 (MQ=255)                                                                                   ccTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAgg                                                             <  5:266710/149‑1 (MQ=255)                                                                                         gCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTgc                                                       <  2:233857/149‑1 (MQ=255)                                                                                         gCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTgc                                                       <  8:17703/149‑1 (MQ=255)                                                                                               cTCGAATCAGGCGTGCCAGGACCGGGTCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCg                                                 <  7:176596/149‑1 (MQ=255)                                                                                                  gAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGcc                                              <  6:144215/149‑1 (MQ=255)                                                                                                  gAATCAGGCGTGCCAGTACCGGGGCGGCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCgg                                                <  8:267325/147‑1 (MQ=255)                                                                                                         gCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCt                                                         <  4:30464/131‑1 (MQ=255)                                                                                                         gCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCt                                                         >  3:30464/1‑131 (MQ=255)                                                                                                          cGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTg                                      <  1:268258/149‑1 (MQ=255)                                                                                                           gTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCATGCTAAGGGCCTGCGCATCGGCCGACACCTgg                                     <  5:108521/149‑1 (MQ=255)                                                                                                           gTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTgg                                     <  1:369513/149‑1 (MQ=255)                                                                                                              ccAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCa                                  <  3:87946/149‑1 (MQ=255)                                                                                                                  tACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCa                              <  5:199888/149‑1 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGAt                                                                                                              >  8:280146/1‑64 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGAt                                                                                                              <  7:280146/64‑1 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGc                                                                      <  6:103371/104‑1 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGc                                                                      >  5:103371/1‑104 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAgg                                                             >  3:203888/1‑113 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAgg                                                             <  4:203888/113‑1 (MQ=255)                                                                                                                       gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCa                         <  4:7471/149‑1 (MQ=255)                                                                                                                         gCGTCACCGCCGCTAGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGTCCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAgg                       <  3:61900/149‑1 (MQ=255)                                                                                                                         gCGTCACCGCCACTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAgg                       >  8:123831/1‑149 (MQ=255)                                                                                                                            tCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAgg                       <  3:162201/146‑1 (MQ=255)                                                                                                                            tCACAGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGCAACATCACCGTACACAGGCCCGGCGCATCCACCACAATCAACGCCAGCCTAAGGGCCTGCGCATAGGCCGACAACTGTGCATACAGCCCAgg                       >  4:162201/1‑146 (MQ=255)                                                                                                                               ccgccgCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtcgt                  >  8:125753/1‑148 (MQ=255)                                                                                                                               ccgccgCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtcgt                  <  7:125753/148‑1 (MQ=255)                                                                                                                                 gccgcTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtc                    >  6:191504/1‑144 (MQ=255)                                                                                                                                 gccgcTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtc                    <  5:191504/144‑1 (MQ=255)                                                                                                                                 gccgcTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGgc                                                                                           >  3:351575/1‑73 (MQ=255)                                                                                                                                 gccgcTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGgc                                                                                           <  4:351575/73‑1 (MQ=255)                                                                                                                                  ccgcTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGt              >  2:262196/1‑149 (MQ=255)                                                                                                                                        gAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTcgccgc        >  5:43582/1‑149 (MQ=255)                                                                                                                                         aaaCAGCGATTACCA‑GCGGCGAACGGTCAACGATCTCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCTGCAGGCTAAGGGCCTGCTCCACGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGAcgccgc        >  8:125380/1‑148 (MQ=255)                                                                                                                                             aGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTagg                                                              <  2:336746/90‑3 (MQ=255)                                                                                                                                             aGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTagg                                                              >  1:336746/1‑88 (MQ=255)                                                                                                                                              gCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAc  <  7:123831/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                ACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAC  >  NC_002947/4508266‑4508555

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