Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,508,406	A→C	9.4%	A125A (GCT→GCG)	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,406	0	A	C	9.4%	185.4 / 10.1	64	A125A (GCT→GCG)	cysG	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base C (3/3);  total (29/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.13e-01

GTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGC  >  NC_002947/4508262‑4508549                                                                                                                                                 |                                                                                                                                                 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:108521/1‑149 (MQ=255)                                          gTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCAcc                                                                                                    >  6:54408/1‑149 (MQ=255)                                          gTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCAcc                                                                                                    >  7:17703/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                                                                                                                    >  6:56873/1‑113 (MQ=38)                                                              tCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGAt                                                                                                                    <  5:56873/113‑1 (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                                                 <  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                                                   >  3:30464/1‑131 (MQ=255)                                                                                                             gCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCt                                                   <  4:30464/131‑1 (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            <  7:125753/148‑1 (MQ=255)                                                                                                                                   ccgccgCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtcgt            >  8:125753/1‑148 (MQ=255)                                                                                                                                     gccgcTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtc              <  5:191504/144‑1 (MQ=255)                                                                                                                                     gccgcTGGTAATCGCTGTTTCCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGgtc              >  6:191504/1‑144 (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)                                                                                                                                                 |                                                                                                                                                 GTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGC  >  NC_002947/4508262‑4508549

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