Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,508,409	G→T	9.1%	I124I (ATC→ATA) ‡	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,409	0	G	T	9.1%	196.7 / 9.3	66	I124I (ATC→ATA) ‡	cysG	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base T (3/3);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

GTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAC  >  NC_002947/4508262‑4508555                                                                                                                                                    |                                                                                                                                                    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                                                                                                              <  7:280146/64‑1 (MQ=255)                                                                                                                           gggCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGAt                                                                                                              >  8:280146/1‑64 (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                    >  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                                                              >  1:336746/1‑88 (MQ=255)                                                                                                                                                 aGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTagg                                                              <  2:336746/90‑3 (MQ=255)                                                                                                                                                  gCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAc  <  7:123831/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GTTGACGTCCGGGTACAGGGATTTGACCTTGTGCCGGAAGCGTGCAGCCAGCCCGGCCAGCTCACCATAGGCCGAAGGGATCCACGCCTCGAGCTTGGCTCGAATCAGGCGTGCCAGTACCGGGGCGTCACCGCCGCTGGAAACAGCGATTACCA‑GCGGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAC  >  NC_002947/4508262‑4508555

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