Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,950,984	C→A	10.7%	R658L (CGT→CTT)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,950,984	0	C	A	10.7%	181.4 / 9.8	63	R658L (CGT→CTT)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base A (3/4);  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 = 9.96e-01

CTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGA  >  NC_002947/5950839‑5951127                                                                                                                                                   |                                                                                                                                                cTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGcg                                                                                                                                              >  6:274185/1‑148 (MQ=255) cTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGc                                                                                                                                               <  2:81621/148‑1 (MQ=255)  tGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGcgt                                                                                                                                             >  6:12997/1‑147 (MQ=255)    gTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAATGTCCATGCGTGGGTATTGTCGGTCATGTAGCCTTCCAGCTCCTTCTTCACCCGCATCAACAGTTTCTTTCCCTTCCACCGGGAAGCAGGTCTCGACACGCtt                                                                                                                                             <  1:33718/148‑1 (MQ=255)    gTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGt                                                                                                                                           >  2:77751/1‑149 (MQ=255)    gTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCgtgt                                                                                                                                           >  1:196869/1‑149 (MQ=255)     ttGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTc                                                                                                                                          <  1:119794/149‑1 (MQ=255)     ttGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTc                                                                                                                                          >  3:258261/1‑149 (MQ=255)     ttGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTc                                                                                                                                          >  6:239838/1‑149 (MQ=255)        ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGa                                                                                                                                        <  5:274185/148‑1 (MQ=255)         cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAga                                                                                                                                      <  2:196869/149‑2 (MQ=255)           gTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgacc                                                                                                                                    <  5:12997/149‑4 (MQ=255)                  tgctgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                                             <  4:258261/149‑1 (MQ=255)                   gctgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcca                                                                                                                            <  4:263702/149‑1 (MQ=255)                       cgcACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGCATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTTCGCtccatc                                                                                                                          <  8:132171/147‑1 (MQ=255)                          aCGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                      >  5:61514/1‑132 (MQ=255)                          aCGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                      <  6:61514/132‑1 (MQ=255)                             tAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgc                                                                                                                  <  1:89369/149‑1 (MQ=255)                                            tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGa                                                                                                            <  4:104548/140‑1 (MQ=255)                                            tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGa                                                                                                            >  3:104548/1‑140 (MQ=255)                                            tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAgg                                                                                                        >  8:184044/1‑144 (MQ=255)                                               aaGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTg                                                                                                >  3:20452/1‑149 (MQ=255)                                                 ggTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACACGCAGCAACAGTTTC‑TTGCCTTCCAACGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGCTCTGct                                                                                              >  8:140220/1‑149 (MQ=255)                                                  gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGc                                                                                               <  5:95017/147‑1 (MQ=255)                                                  gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGc                                                                                               >  6:95017/1‑147 (MQ=255)                                                    ccAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGAGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGctct                                                                                            >  5:25731/1‑148 (MQ=255)                                                              gTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCAtt                                                                                 <  6:25731/149‑1 (MQ=255)                                                                 ttGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                                               <  5:239838/148‑1 (MQ=255)                                                                 ttGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCACCGCCATTGAg                                                                              >  2:205245/1‑149 (MQ=255)                                                                      ggTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAg                                                                          <  4:20452/148‑1 (MQ=255)                                                                      ggTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAg                                                                          >  1:239578/1‑148 (MQ=255)                                                                       gTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTa                                                                        <  1:86140/149‑1 (MQ=255)                                                                            gTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAaca                                                                   >  1:311759/1‑149 (MQ=255)                                                                             tAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTg                                                                                >  7:152057/1‑135 (MQ=255)                                                                             tAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTg                                                                                <  8:152057/135‑1 (MQ=255)                                                                              aGCCTTCCAGCCCCTTCTTCACCCGCAGCAACACTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGATTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGACATTGAGGAAGTAGAACAcc                                                                 >  8:8451/1‑149 (MQ=255)                                                                                  ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                       <  5:182122/75‑1 (MQ=255)                                                                                  ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                       >  6:182122/1‑75 (MQ=255)                                                                                         tCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTc                                                      >  2:189400/1‑149 (MQ=255)                                                                                           cttcttCACCCGCAGCAACCGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                                                    <  7:173580/149‑1 (MQ=255)                                                                                           cttcttCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                                                    <  7:184044/149‑1 (MQ=255)                                                                                             tcttcACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAg                                                   >  2:71097/1‑148 (MQ=255)                                                                                              cttcACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTg                                                                                <  6:220469/118‑1 (MQ=255)                                                                                              cttcACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAgg                                                  >  3:231188/1‑148 (MQ=255)                                                                                              cttcACCCGCAGCAACAGTTTC‑TTGCCTTCCACCCGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTg                                                                                >  5:220469/1‑118 (MQ=255)                                                                                                    ccGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAgg                                                                             <  8:223586/115‑1 (MQ=255)                                                                                                    ccGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAgg                                                                             >  7:223586/1‑115 (MQ=255)                                                                                                    ccGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGc                                           <  2:307897/149‑1 (MQ=255)                                                                                                           aaCAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTACTCGAGGAAGCGGCCAATGAt                                    >  6:117473/1‑149 (MQ=255)                                                                                                              aGTTTC‑TTGCCTTCCACCGGGAAGTAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGa                                 <  1:204944/149‑1 (MQ=255)                                                                                                              aGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGa                                 <  4:231188/149‑1 (MQ=255)                                                                                                                     ttGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCAc                           >  5:120113/1‑149 (MQ=255)                                                                                                                       gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACa                                                                                                          >  6:140213/1‑68 (MQ=255)                                                                                                                       gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACa                                                                                                          <  5:140213/68‑1 (MQ=255)                                                                                                                        ccTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCt                         <  2:239578/148‑1 (MQ=255)                                                                                                                          ttCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGa                      >  4:227763/1‑149 (MQ=255)                                                                                                                              aCCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCGCCTGGAtgtt                  >  1:148072/1‑149 (MQ=255)                                                                                                                                cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttgt                <  7:140220/149‑1 (MQ=255)                                                                                                                                   gAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGc              <  1:276157/148‑1 (MQ=255)                                                                                                                                      gCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGaaa          <  2:311759/149‑1 (MQ=255)                                                                                                                                           tctcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGACGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAccccc     <  1:71097/149‑1 (MQ=255)                                                                                                                                             tcGACACGGTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACACGCGTGTGCTCGAGGAAGCAGCCAATGAACGAGCGCACCTGGATGTTGTGCGAAACCCCCg    >  4:195314/1‑148 (MQ=255)                                                                                                                                               gACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCACCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGa  <  1:205245/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                CTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTC‑TTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGA  >  NC_002947/5950839‑5951127

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