Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,166,929:1	+C	6.6%	coding (591/1212 nt)	creA ←	creatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,166,929	1	.	C	6.6%	243.8 / 10.7	61	coding (591/1212 nt)	creA	creatinase
Reads supporting (aligned to +/- strand):  ref base . (30/27);  new base C (2/2);  total (32/29)
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.21e-01

GGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGAT  >  NC_002947/4166784‑4167078                                                                                                                                                   |                                                                                                                                                       ggATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCa                                                                                                                                                    <  2:746503/150‑1 (MQ=255)  gATCTCGCCCTTGCCCACACGGCGGCTGGTCACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTCGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAg                                                                                                                                                   <  1:804347/150‑1 (MQ=255)   aTCTCGCACTTGCCCACCCGGCGGCTGGTGACTGCGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                                  <  2:11221/150‑1 (MQ=255)    tctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAg                                                                                                                                                   <  2:346307/148‑1 (MQ=255)    tctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAg                                                                                                                                                   >  1:346307/1‑148 (MQ=255)    tctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTg                                                                                                                                                 >  2:124825/1‑150 (MQ=255)     ctcGCCCGTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGGCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCAGCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGc                                                                                                                                                <  2:135134/150‑1 (MQ=255)      tcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTATTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCa                                                                                                                                                         <  2:375189/141‑1 (MQ=255)      tcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTATTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCa                                                                                                                                                         >  1:375189/1‑141 (MQ=255)          ccTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTAGGCGCCACGTGTGTCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCAccc                                                                                                                                           >  2:33927/1‑149 (MQ=255)              gCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTa                                                                                                                                      >  2:1223944/1‑150 (MQ=255)                ccACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACt                                                                                                                                    >  1:936945/1‑150 (MQ=255)                ccACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTAGGTGGCA‑TGCAGTGCCACCTC‑GTACt                                                                                                                                    >  2:590383/1‑150 (MQ=255)                  acacGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCg                                                                                                                                  >  2:1100884/1‑150 (MQ=255)                    acGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCggg                                                                                                                                <  1:1223944/150‑1 (MQ=255)                       gcggcTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCTATCAGTTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCGTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCAc                                                                                                                             <  1:943530/150‑1 (MQ=255)                             ggTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTc                                                                                                                       <  1:551048/150‑1 (MQ=255)                                gACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                                  <  1:326137/121‑1 (MQ=255)                                gACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                                  >  2:326137/1‑121 (MQ=255)                                gACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCACGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACg                                                                                                                    >  1:608319/1‑150 (MQ=255)                                     ggTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCACTGCA‑TGCCAcccag                                                                                                                                         <  2:853764/125‑4 (MQ=255)                                     ggTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCACTGCA‑TGCCAcccag                                                                                                                                         >  1:853764/1‑122 (MQ=255)                                          tgtgCGCCATCTGTGTTGAGCCCCGACTGGAACCAGGTCCAGGTGTCCATCAGCGCGGCGTCGGGGAAGCGTTTGGCAAGTTCGCTCACCATCCCCTGGGTGGCA‑TGCAGTGCCACCTT‑GTACCGGTGCACCTGCTCACGCAGGGCTTCg                                                                                                          <  1:33927/150‑1 (MQ=255)                                          tgtgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCCTCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTGCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCg                                                                                                          <  2:786171/150‑1 (MQ=255)                                                ccATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACt                                                                                                    >  1:835316/1‑150 (MQ=255)                                                      gtgtTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgccgcc                                                                                              >  2:705281/1‑150 (MQ=255)                                                       tgtTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGAc                                                                                                        <  2:403153/139‑1 (MQ=255)                                                       tgtTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCGC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGAc                                                                                                        >  1:403153/1‑139 (MQ=255)                                                         ttGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACCGGGTCCCCTGCTCACGCAGGGCTTCGACCACTGccgccccg                                                                                           >  1:977906/1‑150 (MQ=255)                                                             gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                       >  1:146451/1‑150 (MQ=255)                                                             gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                       >  2:1144092/1‑150 (MQ=255)                                                              cccGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGAt                                                                                      >  2:1178374/1‑150 (MQ=255)                                                               ccGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATg                                                                                     <  1:972253/150‑1 (MQ=255)                                                                cGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCAACTGATCACGCAGGGATTCCGCCATTGCCTCCACGCCCATgg                                                                                    >  2:915121/1‑149 (MQ=255)                                                                 ggACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTc                                                                                   <  2:936945/150‑1 (MQ=255)                                                                  gACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCg                                                                                  <  1:1144092/150‑1 (MQ=255)                                                                      ggAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTAGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCa                                                                              <  1:590383/150‑1 (MQ=255)                                                                       gAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCAc                                                                             >  1:35183/1‑150 (MQ=255)                                                                       gAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCAc                                                                             >  2:1139811/1‑150 (MQ=255)                                                                            aggtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCtgtg                                                                        >  2:673224/1‑150 (MQ=255)                                                                                  aggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGt                                                                  <  2:835316/150‑1 (MQ=255)                                                                                  aggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGt                                                                  <  1:1178374/150‑1 (MQ=255)                                                                                         cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                           >  1:508854/1‑150 (MQ=255)                                                                                           tCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATaa                                                         >  2:345892/1‑150 (MQ=255)                                                                                                tCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGcc                                                    >  2:678419/1‑150 (MQ=255)                                                                                                 cGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCt                                                   <  1:705281/150‑1 (MQ=255)                                                                                                         gggAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCg                                           <  1:1111054/150‑1 (MQ=255)                                                                                                               cGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCTGCCACCTGTGCACCGTTACGGATAAGTGCCTGCTCCTCGGCCGAc                                     <  2:146451/150‑1 (MQ=255)                                                                                                                gTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACt                                    >  1:138460/1‑150 (MQ=255)                                                                                                                ctttGGCAGTATCGCCCCCAATCGCCTCGGTCCTC‑AGCAGTGCCCCCAG‑TTACTGGCGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACt                                    <  2:977906/149‑1 (MQ=255)                                                                                                                 tttGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGc                                         >  1:188781/1‑144 (MQ=255)                                                                                                                 tttGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGc                                         <  2:188781/144‑1 (MQ=255)                                                                                                                         tttCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgc                           <  2:508854/150‑1 (MQ=255)                                                                                                                         tttCGCGCACCATCGCCTGGGTGGCA‑TGCAGTCCAACGTCACAACCAAGGCTTC‑GCCCACGCAGGGCTTTGAGCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgc                           >  1:542808/1‑150 (MQ=255)                                                                                                                                 acaaTCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCGTGCTCCTTGGCCCACTTGATCATGCGCATACGc                   >  1:598755/4‑150 (MQ=255)                                                                                                                                       gCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCAc             >  1:7014/1‑150 (MQ=255)                                                                                                                                       gCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGCGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCAc             >  2:1133064/1‑150 (MQ=255)                                                                                                                                        ccTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                     >  1:970009/1‑141 (MQ=255)                                                                                                                                        ccTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                     <  2:970009/141‑1 (MQ=255)                                                                                                                                         cTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACgg           <  2:608319/150‑1 (MQ=255)                                                                                                                                            ggTGGCACTGCA‑TGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgcg                          <  2:1217255/132‑1 (MQ=255)                                                                                                                                            ggTGGCACTGCA‑TGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgcg                          >  1:1217255/1‑132 (MQ=255)                                                                                                                                                gCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCg    <  1:1100884/150‑1 (MQ=255)                                                                                                                                                  a‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGAt  <  1:1195182/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGAT  >  NC_002947/4166784‑4167078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑