Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,496,997:1	+C	11.3%	coding (70/1212 nt)	cobL →	precorrin‑6Y C(5,15)‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,496,997	1	.	C	11.3%	212.1 / 25.3	62	coding (70/1212 nt)	cobL	precorrin‑6Y C(5,15)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (28/27);  new base C (3/4);  total (31/31)
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.01e-01

TGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGTGTG  >  NC_002947/5496851‑5497139                                                                                                                                                     |                                                                                                                                                tgctgcACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGgc‑‑g                                                                                                                                               >  1:58163/1‑149 (MQ=255)           tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCctgct                                                                                                                                      <  2:58163/147‑1 (MQ=255)           tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTgg                                                                                                                                    <  5:248216/149‑1 (MQ=255)            aCCTGTAGGAGCGGGCTTGCCCGCGAAGGGGC‑GCTAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTgg                                                                                                                                    <  4:25134/148‑1 (MQ=255)              cTGTAGGAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGgcg                                                                                                                                 >  5:165891/1‑149 (MQ=255)                   ggAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgc                                                                                                                              <  8:8433/149‑3 (MQ=255)                    gAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGcc                                                                                                                           >  3:37009/1‑149 (MQ=255)                      gCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCgg                                                                                                                         <  6:306218/149‑1 (MQ=255)                      gCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGACCCGGCCAGCGCCCAGCAGGgc                                                                                                                                  >  5:64303/1‑142 (MQ=255)                      gCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGACCCGGCCAGCGCCCAGCAGGgc                                                                                                                                  <  6:64303/142‑1 (MQ=255)                           cgtGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAACATGGCCGACTGAAATACCTGACACCCTCGCTGACACTAGTAGGCATCGCCGAAGCCGGCGTCCGCGGCCTGCGCAACCAGCCCCGGC‑‑GCCCCCTGCCGCGCGCTGGCCGGGTCtt                                                                                                                    <  2:127705/147‑1 (MQ=255)                              gCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                 >  8:110854/1‑149 (MQ=255)                                ccGTGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGCCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCa                                                                                                               <  1:98029/149‑1 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgccg                                                                                                                            >  3:126760/1‑134 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgccg                                                                                                                            <  4:126760/138‑5 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgcc                                                                                                                             <  5:145646/137‑4 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCCAGCGCCCAGCAGGgcgcgcc                                                                                                                             >  6:145646/1‑134 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                 >  1:35041/1‑147 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                 >  4:270916/1‑147 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                 >  8:269290/1‑147 (MQ=255)                                ccGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCa                                                                                                               >  1:224163/1‑149 (MQ=255)                                          gc‑gcGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCgcgc                                                                                                       >  1:229636/1‑147 (MQ=255)                                             gcgAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGc                                                                                                   >  3:169558/1‑149 (MQ=255)                                                      ccccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgc                                                                                          <  6:165891/149‑1 (MQ=255)                                                       cccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCAGCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgct                                                                                         <  7:269290/149‑1 (MQ=255)                                                        ccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGc                                                                                                                >  5:251284/1‑125 (MQ=255)                                                        ccAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGc                                                                                                                <  6:251284/125‑1 (MQ=255)                                                         cAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgctgc                                                                                       <  7:110854/149‑1 (MQ=255)                                                          aaaaaTGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgctgc                                                                                       >  6:129615/1‑148 (MQ=255)                                                          aaaaaTGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGACGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTtgctgc                                                                                       >  5:118765/1‑148 (MQ=255)                                                            aaaTGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCAc                                                                                    <  2:229636/149‑1 (MQ=255)                                                                 gCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCTGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACTTTGc                                                                               <  3:61333/149‑1 (MQ=255)                                                                 gCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGc                                                                               <  5:129615/149‑1 (MQ=255)                                                                    gATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTg                                                                            <  2:35041/149‑1 (MQ=255)                                                                         aaaTACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGc                                                                       <  2:41253/149‑1 (MQ=255)                                                                          aaTACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGA‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCAAGCTGGCCTTGCTGCCACGTTGCGTGACCGCt                                                                      >  8:19915/1‑149 (MQ=255)                                                                           aTACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCt                                                                      >  6:48087/1‑148 (MQ=255)                                                                                   cacCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtg                                                             <  3:89463/149‑1 (MQ=255)                                                                                     cccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggc                                                           <  4:37009/149‑1 (MQ=255)                                                                                     cccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggc                                                           <  3:267607/149‑1 (MQ=255)                                                                                      ccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggcc                                                          >  6:185639/1‑149 (MQ=255)                                                                                      ccTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTtggcc                                                          >  8:256135/1‑149 (MQ=255)                                                                                         ggCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                 >  4:266363/1‑91 (MQ=255)                                                                                         ggCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCgg                                                                                                                 <  3:266363/91‑1 (MQ=255)                                                                                                 gtagtaGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGccc                                               <  2:227072/149‑1 (MQ=255)                                                                                                     tagGCATCGACGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTtc                                           <  6:118765/149‑1 (MQ=255)                                                                                                         cATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGc                                       <  7:256135/149‑1 (MQ=255)                                                                                                          aTCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGccc                                                    >  1:90133/1‑135 (MQ=255)                                                                                                          aTCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGccc                                                    <  2:90133/135‑1 (MQ=255)                                                                                                             ggCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTgg                                    <  3:270916/148‑1 (MQ=255)                                                                                                             ggCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGc                                   >  1:40602/1‑149 (MQ=255)                                                                                                                      ggCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCTCACTGCTGGGCGCTGGACGGGTCTTCGGCAGCGCTCGCCAGCTTGCCTTGCTGCCACGTTGCGTGAGCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCt                          <  7:19915/149‑1 (MQ=255)                                                                                                                        cTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTgg                        >  8:98627/1‑149 (MQ=255)                                                                                                                            aGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCt                    >  1:170384/1‑149 (MQ=255)                                                                                                                              cGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCCGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTc                  >  8:77705/1‑149 (MQ=255)                                                                                                                                    gggCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCgcggc             >  2:94980/1‑148 (MQ=255)                                                                                                                                        aaaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGcc        >  2:1791/1‑149 (MQ=255)                                                                                                                                        aaaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGcc        <  6:29182/149‑1 (MQ=255)                                                                                                                                        aaaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGcc        >  5:275668/1‑149 (MQ=255)                                                                                                                                         aaCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCAAGCCg       <  2:40602/149‑1 (MQ=255)                                                                                                                                          aCAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGCACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCgg      >  7:239114/1‑149 (MQ=255)                                                                                                                                             ggCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGgtgt   >  6:276067/1‑149 (MQ=255)                                                                                                                                              gCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGtgtg  >  2:156646/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                TGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGC‑GCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGC‑‑GCGCCCTGCTGGGCGCTGGCCGGGTCTTCGGCAGCCCGCGCCAGCTGGCCTTGCTGCCACGTTGCGTGACCGCTGAACAGTTGGCCTGGCCCAGCCCGTTCTCGCTGGCCCCAGTGCTGGCCCTTCGCGGCGAGCCGGTGTG  >  NC_002947/5496851‑5497139

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