Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,166,934	Δ1 bp	6.9%	coding (586/1212 nt)	creA ←	creatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,166,934	0	G	.	6.9%	177.6 / 6.8	58	coding (586/1212 nt)	creA	creatinase
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base . (2/2);  total (31/27)
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.60e-01

GATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGAT  >  NC_002947/4166785‑4167078                                                                                                                                                       |                                                                                                                                                  gATCTCGCCCTTGCCCACACGGCGGCTGGTCACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTCGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAg                                                                                                                                                   <  1:804347/150‑1 (MQ=255)  aTCTCGCACTTGCCCACCCGGCGGCTGGTGACTGCGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                                  <  2:11221/150‑1 (MQ=255)   tctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAg                                                                                                                                                   >  1:346307/1‑148 (MQ=255)   tctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAg                                                                                                                                                   <  2:346307/148‑1 (MQ=255)   tctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTg                                                                                                                                                 >  2:124825/1‑150 (MQ=255)    ctcGCCCGTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGGCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCAGCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGc                                                                                                                                                <  2:135134/150‑1 (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                                                                                                                                                  >  2:326137/1‑121 (MQ=255)                               gACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                                  <  1:326137/121‑1 (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                                                                                       >  2:1144092/1‑150 (MQ=255)                                                            gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                       >  1:146451/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                                                                             >  2:1139811/1‑150 (MQ=255)                                                                      gAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCAc                                                                             >  1:35183/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)                                                                                                                                                       |                                                                                                                                                  GATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTC‑GTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGAT  >  NC_002947/4166785‑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: ‑