Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,977,578	C→G	6.8%	G1018A (GGC→GCC)	cobN ←	aerobic cobaltochelatase subunit CobN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,977,578	0	C	G	6.8%	190.0 / 3.7	59	G1018A (GGC→GCC)	cobN	aerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base G (2/2);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAG  >  NC_002947/3977437‑3977726                                                                                                                                               |                                                                                                                                                     tcgtggCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTg                                                                                                                                               <  2:133931/149‑1 (MQ=255)   gtcgCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGc                                                                                                                                             <  6:52824/149‑1 (MQ=255)          ccgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACCCCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                      >  8:148391/1‑149 (MQ=255)          ccgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                      <  8:183233/149‑1 (MQ=255)           cgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                      >  1:348334/1‑148 (MQ=255)           cgccgTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGc                                                                                                                                      <  2:348334/148‑1 (MQ=255)                   gCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCgg                                                                                                                              >  4:27543/1‑148 (MQ=255)                       aGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGGACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATc                                                                                                                          >  5:326557/1‑148 (MQ=255)                       aGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCg                                                                                                                                        >  5:230745/1‑134 (MQ=255)                       aGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCg                                                                                                                                        <  6:230745/134‑1 (MQ=255)                          cgccgTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACg                                                                                                                      >  1:290817/1‑149 (MQ=255)                                     aGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGGCTGGCTGGCgg                                                                                                                                           <  8:290738/117‑2 (MQ=255)                                     aGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGGCTGGCTGGCgg                                                                                                                                           >  7:290738/1‑116 (MQ=255)                                       gTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACACGCCAGGCTGGCTGGCgggt                                                                                                                                         >  3:189642/1‑114 (MQ=255)                                       gTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACACGCCAGGCTGGCTGGCgggt                                                                                                                                         <  4:189642/117‑4 (MQ=255)                                         aGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAgc                                                                                                       >  6:21820/1‑149 (MQ=255)                                          gACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAgcg                                                                                                      >  7:121771/1‑149 (MQ=255)                                          gACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAgc                                                                                                       <  3:340909/148‑1 (MQ=255)                                            ccTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGc                                                                                                                                  <  7:32782/119‑1 (MQ=255)                                            ccTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGc                                                                                                                                  >  8:32782/1‑119 (MQ=255)                                             ctctGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAg                                                                                                        <  4:151761/144‑1 (MQ=255)                                             ctctGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAg                                                                                                        >  3:151761/1‑144 (MQ=255)                                                    aGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCcgctc                                                                                            >  2:284419/1‑149 (MQ=255)                                                     ggTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCcgctcg                                                                                           <  3:339906/149‑1 (MQ=255)                                                           cgcgGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTgcgc                                                                                     <  7:148391/149‑1 (MQ=255)                                                            gcgGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGca                                                                                    <  3:27543/149‑1 (MQ=255)                                                                  gtgCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACggg                                                                              >  1:187613/1‑149 (MQ=255)                                                                   tgCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggc                                                                             >  7:103245/1‑149 (MQ=255)                                                                    gCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcg                                                                            >  3:259975/1‑149 (MQ=255)                                                                     ccACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcgg                                                                           >  5:154200/1‑149 (MQ=255)                                                                       acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTACGCGCCCGCTTTGGCACCGAACACCCGCCACCCAGCCTGCCGTGCGGCCTGCTCGGCATCCACGCCTTGACCCTGCACCGCGGCCCGCTCGCTGCCGACACGcgcggc                                                                          >  8:233884/1‑148 (MQ=255)                                                                       acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcggc                                                                          >  6:310916/1‑148 (MQ=255)                                                                       acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGggcggc                                                                          >  2:141069/1‑148 (MQ=255)                                                                       acaAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGcc                                                                         <  2:226324/149‑1 (MQ=255)                                                                         aaaCGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAg                                                                       >  8:234605/1‑149 (MQ=255)                                                                           aCGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAggg                                                                     >  8:70082/1‑149 (MQ=255)                                                                             gCCCGTCGATGGCGTTCTGCACCCCGGCACCAGTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCCTGCGGCCTGCTCGCCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAgggg                                                                    >  1:232786/1‑149 (MQ=255)                                                                               ccGTCGATGGCGTTCTGCACCCCGGCACC‑GTTAGCGCCCGGTTTGGCACCGAACTCCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGtt                                                                 <  3:302458/149‑1 (MQ=255)                                                                                  tCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGATCACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCGACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAg                                                              <  2:232786/149‑1 (MQ=255)                                                                                  tCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAg                                                              <  6:154200/149‑1 (MQ=255)                                                                                  tCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAg                                                              <  7:70082/149‑1 (MQ=255)                                                                                    gATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAgg                                                                      >  1:111326/1‑139 (MQ=255)                                                                                    gATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAgg                                                                      <  2:111326/139‑1 (MQ=255)                                                                                      tGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGTCCCGCTCGCTGCGCACACGggcggc                                                                          <  3:155080/133‑1 (MQ=255)                                                                                      tGGCGTTCTGCACACCGGCACA‑GTAAGCGCCCGGTTTGGCAAAGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCAACCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGcgcggc                                                                          >  4:155080/1‑133 (MQ=255)                                                                                        gCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTc                                                        <  5:21820/149‑1 (MQ=255)                                                                                               gCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCg                                                 >  3:176807/1‑149 (MQ=255)                                                                                                       gCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAgcgc                                         <  5:310916/149‑1 (MQ=255)                                                                                                          cc‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGc                                      >  6:218766/1‑149 (MQ=255)                                                                                                                       ggTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCAc                          >  1:258257/1‑149 (MQ=255)                                                                                                                        gTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCAcg                         >  2:93462/1‑149 (MQ=255)                                                                                                                         tttGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCAcgg                        >  6:161967/1‑149 (MQ=255)                                                                                                                               aCCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGTGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGt                  >  1:229461/1‑149 (MQ=255)                                                                                                                               aCCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGt                  >  7:82201/1‑149 (MQ=255)                                                                                                                                ccGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTc                 <  4:176807/149‑1 (MQ=255)                                                                                                                                    acacCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAAc             <  6:326557/149‑1 (MQ=255)                                                                                                                                      acCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAg           >  7:302322/1‑149 (MQ=255)                                                                                                                                       cccGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAgg          <  1:141069/149‑1 (MQ=255)                                                                                                                                        ccGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGc         >  5:339674/1‑149 (MQ=255)                                                                                                                                          gccagccaGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTgc                                                                                       <  3:116827/69‑1 (MQ=255)                                                                                                                                          gccagccaGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTgc                                                                                       >  4:116827/1‑69 (MQ=255)                                                                                                                                          gcacgccaggcTGGCTGGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAggg                                                                     <  4:315886/83‑1 (MQ=255)                                                                                                                                          gcacgccaggcTGGCTGGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAggg                                                                     >  3:315886/5‑87 (MQ=255)                                                                                                                                               ccagccTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAg  <  2:258257/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     TCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGGTCGTCGCGGCTGTGCCACAAACGCCCGTCGATGGCGTTCTGCACCCCGGCACC‑GTAAGCGCCCGGTTTGGCACCGAACACCCGCCAGCCAGCCTGGCGTGCGGCCTGCTCGGCATCCACGCCTTGAGCCTGCAGCGCGGCCCGCTCGCTGCGCACACGGGCGGCCAGGGGGTTGAGGTCTTCCGGCTCGTCCAGCGCGGCCACGGCCTGCACGGCGGCGTCGAACAGGCGGATCAG  >  NC_002947/3977437‑3977726

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