Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,508,402	A→T	9.5%	S127T (TCC→ACC)	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,402	0	A	T	9.5%	170.9 / 8.2	63	S127T (TCC→ACC)	cysG	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base T (3/3);  total (28/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 = 8.49e-01

CGACGCTGGTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGC  >  NC_002947/4508254‑4508549                                                                                                                                                     |                                                                                                                                                     cgACGCTGGTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGa                                                                                                                                                      <  8:98211/149‑1 (MQ=255)   aCGCTGGTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGaaa                                                                                                                                                    >  3:211435/1‑149 (MQ=255)         gTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGc                                                                                                                                                <  5:21028/147‑1 (MQ=255)         gTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGa                                                                                                                                              <  1:185464/149‑1 (MQ=255)             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                                                                                                    >  6:54408/1‑149 (MQ=255)                                                  gTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCAcc                                                                                                    >  7:17703/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                                                           >  6:266710/1‑148 (MQ=255)                                                                                            aCGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCt                                                           >  5:144215/1‑148 (MQ=255)                                                                                               ccTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAgg                                                       <  5:266710/149‑1 (MQ=255)                                                                                                     gCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTgc                                                 <  8:17703/149‑1 (MQ=255)                                                                                                     gCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTgc                                                 <  2:233857/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                                                                                     <  4:351575/73‑1 (MQ=255)                                                                                                                                             gccgcTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGgc                                                                                     >  3:351575/1‑73 (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)                                                                                                                                                     |                                                                                                                                                     CGACGCTGGTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGC  >  NC_002947/4508254‑4508549

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