Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,402,947	G→C	6.2%	R69R (CGC→CGG)	PP_3881 ←	terminase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,402,947	0	G	C	6.2%	196.4 / 2.3	63	R69R (CGC→CGG)	PP_3881	terminase large subunit
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base C (2/2);  total (30/33)
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

ACCGCTCTCGACGTAGACGACCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCACCGGAGTGC  >  NC_002947/4402803‑4403082                                                                                                                                                  |                                                                                                                                        aCCGCTCTCGACGTAGACGACCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGa                                                                                                                                      <  8:210866/148‑1 (MQ=255)     ctctCGACGTCGACGACCCGGAAGCGGCCATAACCGTAACTGCACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGaa                                                                                                                                 <  2:145758/149‑1 (MQ=255)     ctctCGACGTAGACGACCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACCATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGaa                                                                                                                                 <  2:323900/149‑1 (MQ=255)     ctctCCACGTAGACGACCCGGAAGCGGCTATAACCGTCACTGGCCTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGaa                                                                                                                                 <  3:272151/149‑1 (MQ=255)      tctcGACGTAGACGACCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTCCGATAGAAg                                                                                                                                >  4:13706/1‑149 (MQ=255)      tctcGACGTAGACGACCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTCCGATAGAAg                                                                                                                                >  4:8773/1‑149 (MQ=255)                 acgacCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCttt                                                                                                                     <  8:47871/149‑1 (MQ=255)                   gacCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTttgt                                                                                                                   >  5:317161/1‑149 (MQ=255)                   gacCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTttgt                                                                                                                   >  2:11143/1‑149 (MQ=255)                     cccGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAg                                                                                                                         >  7:4522/1‑141 (MQ=255)                     cccGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAg                                                                                                                         <  8:4522/141‑1 (MQ=255)                      ccGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCtt                                                                                                                      <  2:21128/143‑1 (MQ=255)                      ccGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCtt                                                                                                                      >  1:21128/1‑143 (MQ=255)                        ggAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGcc                                                                                                              >  8:287225/1‑149 (MQ=255)                                cgATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCaa                                                                                                      <  7:833/149‑1 (MQ=255)                                   tAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCaa                                                                                                      >  4:317918/1‑146 (MQ=255)                                   tAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCaa                                                                                                      <  3:317918/146‑1 (MQ=255)                                      ccGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGc                                                                                                                               >  4:238630/1‑118 (MQ=255)                                      ccGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGc                                                                                                                               <  3:238630/118‑1 (MQ=255)                                       cGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGTCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAg                                                                                               <  2:67283/149‑1 (MQ=255)                                       cGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAg                                                                                               >  5:124097/1‑149 (MQ=255)                                       cGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAg                                                                                               >  1:165218/1‑149 (MQ=255)                                       cGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCAAGCGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGATTAGCTTTGATGGCCGCATACAAGTCCAAg                                                                                               >  3:44447/1‑149 (MQ=255)                                         tCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGt                                                                                             <  6:160388/149‑1 (MQ=255)                                                aCTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTttgtt                                                                                                                  >  3:64900/1‑121 (MQ=255)                                                aCTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTttgtt                                                                                                                  <  4:64900/121‑1 (MQ=255)                                                  ttCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCcgc                                                                                     >  7:16515/1‑148 (MQ=255)                                                     cACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCtt                                                                                 >  3:34467/1‑149 (MQ=255)                                                      aCCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTc                                                                                >  5:84643/1‑149 (MQ=255)                                                        ccGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCgg                                                                              <  3:110454/149‑1 (MQ=255)                                                        ccGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCgg                                                                              <  6:317161/149‑1 (MQ=255)                                                          gAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTG‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGcc                                                                            >  7:57735/1‑149 (MQ=255)                                                             gATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCt                                                                          >  3:38968/1‑148 (MQ=255)                                                                   gcCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTCCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCa                                                                   <  3:8773/149‑1 (MQ=255)                                                                   gcCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTCCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCa                                                                   <  3:13706/149‑1 (MQ=255)                                                                    ccGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGACGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCATtactc                                                                     >  4:111487/1‑146 (MQ=255)                                                                           tAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCg                                                            >  7:36903/1‑148 (MQ=255)                                                                           tAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGc                                                           <  4:38968/149‑1 (MQ=255)                                                                                    gCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGgcgc                                                  >  8:182740/1‑149 (MQ=255)                                                                                     ccAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGgc                                                    <  6:124097/146‑1 (MQ=255)                                                                                           cagcagCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGcgcacg                                           <  2:165218/149‑1 (MQ=255)                                                                                                gcTCGAACGGTAGCCCTTCAAACTCTCCACCGTTCAAACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCAtt                                       >  1:58923/1‑149 (MQ=255)                                                                                                gcTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCAttt                                      <  8:36903/149‑1 (MQ=255)                                                                                                 cTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTc                                     <  1:11143/149‑1 (MQ=255)                                                                                                    gAACGGTAGCCCTTCAAACTCTCGACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTtcctcc                                                                    <  2:188757/115‑1 (MQ=255)                                                                                                    gAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCATGTCAGCCCGCGCTTCGGCCCTTCCTCCa                                                                   <  1:194814/116‑1 (MQ=255)                                                                                                    gAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTtcctcc                                                                    >  1:188757/1‑115 (MQ=255)                                                                                                    gAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGaa                                  <  2:58923/149‑1 (MQ=255)                                                                                                      aCGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAat                                 <  8:57735/148‑1 (MQ=255)                                                                                                        ggTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATc                              >  8:307522/1‑149 (MQ=255)                                                                                                        ggTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATc                              <  4:34467/149‑1 (MQ=255)                                                                                                            gCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTgg                           <  7:182740/148‑1 (MQ=255)                                                                                                               cTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCg                       <  8:16515/149‑1 (MQ=255)                                                                                                               cTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCg                       >  7:176238/1‑149 (MQ=255)                                                                                                               cTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCg                       >  8:251547/1‑149 (MQ=255)                                                                                                                        ctcCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCAAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCACGCAGATGGCGCGCGCACGCATTTCGAATATCTTGCCCGTCAACACg               >  6:120076/1‑148 (MQ=255)                                                                                                                         tcCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCt             >  2:301317/1‑149 (MQ=255)                                                                                                                           caacGTTG‑AACTTGATGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCTCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGAAACACGCt             <  5:120076/144‑1 (MQ=255)                                                                                                                           cACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTTGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCa           >  3:237562/1‑149 (MQ=255)                                                                                                                              cGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCACCg        <  7:30607/149‑1 (MQ=255)                                                                                                                                  c‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCACCGGAg     <  7:287225/148‑1 (MQ=255)                                                                                                                                     aCTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCACCGGAGTGc  <  8:176238/149‑1 (MQ=255)                                                                                                                                     aCTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTAGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCACCGGAGTg   >  1:57844/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                        ACCGCTCTCGACGTAGACGACCCGGAAGCGGCGATAACCGTCACTGGACTTCCACCCGAAGATGCTGCCGACGATAAATTTCTGCCAGGGCAGCAGCTCGAACGGTAGCCCTTCAAACTCTCCACCGTTC‑AACTTGAGGACGGTGCGATAGAAGCGAATAGCTTTGTTGGCCGCATCCAAGTCCCAGGTCAGCCCGCGCTTCGGCCCTTCCTCCAGATCCCGCAGATGGCGCGCGCACGCATTTCGAATATCTGGCCCGGCAACACGCTCACCGGAGTGC  >  NC_002947/4402803‑4403082

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