Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,919,040	Δ1 bp	6.3%	coding (1439/1557 nt)	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,919,040	0	T	.	6.3%	190.8 / 7.3	63	coding (1439/1557 nt)	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base . (2/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACCACGCCAAGCCATAG  >  NC_002947/4918892‑4919185                                                                                                                                                      |                                                                                                                                                    gCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCt                                                                                                                                                     >  2:280879/1‑149 (MQ=255) gCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCt                                                                                                                                                     >  7:171898/1‑149 (MQ=255)  cccTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTg                                                                                                                                                    >  3:220211/1‑149 (MQ=255)    cTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTgcg                                                                                                                                                  >  1:143729/1‑149 (MQ=255)         gcgcATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTgcgc‑gcgc                                                                                                                                            <  2:384748/149‑1 (MQ=255)                  gCTGGCCTTCGACCTGTCCCCGCTCGGCCACCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGC‑AGGCGCGC‑GCGCAGCGCGccaagctgg                                                                                                                                  <  6:334036/149‑8 (MQ=255)                   cTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGc                                                                                                                                  >  2:268353/1‑149 (MQ=255)                         ttCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCgggg                                                                                                                            >  6:136680/1‑149 (MQ=255)                         ttCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCgggg                                                                                                                            >  5:233640/1‑149 (MQ=255)                         ttCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGC‑AGGCGCGC‑GCGCAGCGCGccaagctgg                                                                                                                                  <  1:358935/142‑8 (MQ=255)                         ttCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGC‑AGGCGCGC‑GCGCAGCGCGccaagctgg                                                                                                                                  >  2:358935/1‑135 (MQ=255)                          tCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGc                                                                                                                           >  7:39274/1‑149 (MQ=255)                           cGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCggg                                                                                                                             <  1:268353/147‑1 (MQ=255)                             aCCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTg                                                                                                                        <  5:106594/149‑1 (MQ=255)                              ccTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTgg                                                                                                                       <  1:124103/149‑1 (MQ=255)                                tGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGCCCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTgc                                                                                                                       >  2:259224/1‑146 (MQ=255)                                  tCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGaa                                                                                                                     <  3:917/147‑1 (MQ=255)                                  tCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGaa                                                                                                                     <  5:20261/147‑1 (MQ=255)                                  tCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGGTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGt                                                                                                                   >  4:305359/1‑149 (MQ=255)                                    cccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTc                                                                                                                  <  3:65522/148‑1 (MQ=255)                                    cccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCg                                                                                                                 >  7:212411/1‑149 (MQ=255)                                    cccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCg                                                                                                                 >  6:194991/1‑149 (MQ=255)                                     ccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCgg                                                                                                                <  8:171898/149‑1 (MQ=255)                                      agCTCGGCCCCCTTCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCAAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGTCCCGGGGCCTGGAAGTCGGc                                                                                                               <  6:248045/148‑1 (MQ=255)                                        cTCGGCCCCCTGCAAGTGCAGGCTCAACTGAGCGAGGGCCGCCTGAGCGGGCAGGTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCAGCTGGCGCGGGGCCTGGATGTCGGCGa                                                                                                             >  4:185489/1‑149 (MQ=255)                                         tCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGAc                                                                                                            >  2:197857/1‑149 (MQ=255)                                          cGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGAcc                                                                                                           >  8:40104/1‑149 (MQ=255)                                            gCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTg                                                                                                         <  2:292869/149‑1 (MQ=255)                                            gCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTg                                                                                                         <  1:197857/149‑1 (MQ=255)                                              ccccTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTTGGCGACCTGGa                                                                                                       <  3:65370/149‑1 (MQ=255)                                              ccccTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGa                                                                                                       <  2:199954/149‑1 (MQ=255)                                                 cTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGTGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CGCAGCTTGGCGGGCTGCGC‑GCGCGCCTGCTGGCCGGGGGCCTGGAAGTCGGCGACCTGGAGTg                                                                                                    <  3:185489/149‑1 (MQ=255)                                                             gCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGt                                                                                                     <  1:380429/136‑1 (MQ=255)                                                                   cTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑Cccgcc                                                                                 <  8:39274/149‑1 (MQ=255)                                                                    tGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑Cccgccg                                                                                >  4:311045/1‑149 (MQ=255)                                                                     gAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑Cccgccg                                                                                >  8:222508/1‑148 (MQ=255)                                                                          aGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAgggg                                                                          <  4:220211/149‑1 (MQ=255)                                                                                   tGAGCGGGCAGCTATGGGCGGAGCGTGAGCAGACCGCGCAGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAgggg                                                                          >  4:246394/1‑140 (MQ=255)                                                                                   tGAGCGGGCAGCTATGGGCGGAGCGTGAGCAGACCGCGCAGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAgggg                                                                          <  3:246394/140‑1 (MQ=255)                                                                                     aGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACg                                                               <  3:220275/149‑1 (MQ=255)                                                                                       cGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGa                                                             <  8:150092/149‑1 (MQ=255)                                                                                        gggCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGAt                                                            <  1:280879/149‑1 (MQ=255)                                                                                        gggCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGAt                                                            <  6:333872/149‑1 (MQ=255)                                                                                         ggCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCCGCCGCAGGGGCCACGCACACGGAt                                                            <  1:204899/149‑1 (MQ=255)                                                                                          gCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGTTGGCCCGGGGTCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATgg                                                          <  1:118651/149‑1 (MQ=255)                                                                                           cAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGa                                                         <  5:136680/149‑1 (MQ=255)                                                                                                    gCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTg                                                >  1:352858/1‑149 (MQ=255)                                                                                                    gCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTg                                                >  6:318333/1‑149 (MQ=255)                                                                                                          cgcgAGCAGACGGAGCGGTTGATCGACA‑CCGAGCTTGGCGCGCTGCGC‑GCGCGGCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCGACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTgg                                           <  2:143729/148‑1 (MQ=255)                                                                                                           gcgAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGAc                                         >  1:178214/1‑149 (MQ=255)                                                                                                               gCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGAc                                         <  6:286639/145‑1 (MQ=255)                                                                                                               gCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGAc                                         >  5:286639/1‑145 (MQ=255)                                                                                                                    ccGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCa                                <  7:263819/149‑1 (MQ=255)                                                                                                                     cgcgcgGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCAt                               <  6:233640/149‑1 (MQ=255)                                                                                                                      tcgcgGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTCCTGGCCCGGGGCCAGGAACTCGGCGACCTGGAGTGCCACCCGGGCAG‑CCCGCCGCAGGGGCCACGCACACGGAGCGAACAACGCTGCGTGGACGAAACCGCAt                               >  6:377435/2‑148 (MQ=255)                                                                                                                              gATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATa                       <  8:212411/148‑1 (MQ=255)                                                                                                                               aTCGACA‑CCCAGCTTGGCGCGCTGAGC‑GCGCGCCTGCTGGCCCGGGGCCAGGAACTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATAc                      >  5:311981/1‑148 (MQ=255)                                                                                                                                     a‑cCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGAcacca               >  2:205767/1‑149 (MQ=255)                                                                                                                                     a‑cCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAAACGCATGAACGATACGAcacca               >  5:117325/1‑149 (MQ=255)                                                                                                                                        ccAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGCCCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGGGGAGGAAACCGCATGAACGATACGCCACCACg             >  6:263091/1‑149 (MQ=255)                                                                                                                                        ccAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACCACg             >  7:145151/1‑149 (MQ=255)                                                                                                                                          aGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGACACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGCTACGACACCACGc            >  5:271146/1‑148 (MQ=255)                                                                                                                                              aggCGCGC‑GCGCAGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCCGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACCACGCCAAGc       >  8:171999/2‑149 (MQ=255)                                                                                                                                               ggCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCGACCCGGGCAT‑CCCGCCGCAGGGGCGACGCACACTGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACCACGCCAAGc       <  6:311981/148‑1 (MQ=255)                                                                                                                                                   cgcTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGTTTGACGAAACCGCATGAACGATACGACACCACGCCAAGCCATAg  <  2:352858/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACA‑CCCAGCTTGGCGCGCTGCGC‑GCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAT‑CCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACCACGCCAAGCCATAG  >  NC_002947/4918892‑4919185

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