Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,103,203	T→C	5.7%	T15A (ACG→GCG)	PP_3611 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,103,203	0	T	C	5.7%	211.3 / 2.7	70	T15A (ACG→GCG)	PP_3611	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (2/2);  total (36/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 = 1.89e-01

CCTTCTTCTCGTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGA  >  NC_002947/4103056‑4103348                                                                                                                                                      |                                                                                                                                                  ccTTCTTCTCGTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTg                                                                                                                                                  >  6:199680/1‑149 (MQ=255)  cttcttCTCGTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTg                                                                                                                                                  <  5:425810/148‑1 (MQ=255)   ttcttcTCGTCAGGGCCTTTCACCGTGACTTGTCGTGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTACCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTgg                                                                                                                                                 >  1:268810/1‑148 (MQ=255)   ttcttcTCGTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTgg                                                                                                                                                 >  3:224798/1‑148 (MQ=255)           gTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAggg                                                                                                                                           <  4:332904/146‑1 (MQ=255)           gTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTa                                                                                                                                        >  4:100163/1‑149 (MQ=255)             cAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACa                                                                                                                                      <  4:224798/149‑1 (MQ=255)              aGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAg                                                                                                                                     <  5:199680/149‑1 (MQ=255)              aGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCA‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTa                                                                                                                                        >  5:390122/1‑146 (MQ=255)                    tttCACCGTGACTTGTCGGGTGGCCAGGTCGATAGGCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAc                                                                                                                                  >  1:61106/1‑147 (MQ=255)                      tCACCGTGACTTGTCGTGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagc                                                                                                                              <  2:268810/148‑1 (MQ=255)                      tCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagc                                                                                                                              <  3:100163/148‑1 (MQ=255)                      tCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagc                                                                                                                              <  5:368674/148‑1 (MQ=255)                        agtgTGACTTGTCGGGTGTCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAg                                                                                                                                     <  8:170730/136‑1 (MQ=255)                        aCCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagcga                                                                                                                            <  4:261111/148‑1 (MQ=255)                        aCCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGACTGTACCGTGAAGCTGGCAGCAGAGGCAAGCGTGGCAAGGGTTACAGCagatcggg                                                                                                                            >  7:170730/1‑141 (MQ=255)                         ccGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagcga                                                                                                                            >  5:372775/1‑147 (MQ=255)                          cGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagcga                                                                                                                            <  7:77329/146‑1 (MQ=255)                          cGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCAcgagcga                                                                                                                            >  8:77329/1‑146 (MQ=255)                           gTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGagag                                                                                                                         >  5:254143/1‑148 (MQ=255)                            tGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTa                                                                                                                       <  6:390122/149‑1 (MQ=255)                            tGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTa                                                                                                                       >  7:182855/1‑149 (MQ=255)                              aCTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTAc                                                                                                                      >  8:81532/1‑148 (MQ=255)                               cTTGTCGGGTGGCCAGGTCGATA‑GCCAGTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTAc                                                                                                                      >  1:4275/1‑148 (MQ=255)                                   tCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATg                                                                                                                 <  7:81532/148‑1 (MQ=255)                                      ggTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGt                                                                                                             >  1:220037/1‑149 (MQ=255)                                        tGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGACCGAGAGTACGAATGGGGtt                                                                                                            >  8:184390/1‑148 (MQ=255)                                                  gATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGc                                                                                                 >  3:216682/1‑149 (MQ=255)                                                        cc‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTAc                                                                                                                      >  4:284035/1‑123 (MQ=255)                                                        cc‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTAc                                                                                                                      <  3:284035/123‑1 (MQ=255)                                                         c‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCt                                                                                           <  2:220037/149‑1 (MQ=255)                                                         c‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCt                                                                                           <  4:370050/149‑1 (MQ=255)                                                         c‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGCCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCt                                                                                           >  1:54408/1‑149 (MQ=255)                                                           gTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGCGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCtg                                                                                          >  3:206649/1‑149 (MQ=255)                                                               ccTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGAt                                                                                      >  6:164135/1‑149 (MQ=255)                                                                cTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGAt                                                                                      >  5:15302/1‑148 (MQ=255)                                                                   tgtgGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGCGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTcc                                                                                             >  5:417524/1‑138 (MQ=255)                                                                   tgtgGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGCGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTcc                                                                                             >  5:411819/1‑138 (MQ=255)                                                                   tgtgGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGCGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTcc                                                                                             <  6:411819/138‑1 (MQ=255)                                                                   tgtgGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGCGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTcc                                                                                             <  6:417524/138‑1 (MQ=255)                                                                     tggtgACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGACAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTc                                                                                <  8:51843/149‑1 (MQ=255)                                                                       gtgACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGc                                                                                                     >  7:22306/1‑126 (MQ=255)                                                                       gtgACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGc                                                                                                     <  8:22306/126‑1 (MQ=255)                                                                          aCATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGaa                                                                            <  4:216682/148‑1 (MQ=255)                                                                            aTGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGc                                                                                                                                    >  5:347010/1‑90 (MQ=255)                                                                            aTGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGc                                                                                                                                    <  6:347010/90‑1 (MQ=255)                                                                            aTGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAg                                                                         >  6:82149/1‑149 (MQ=255)                                                                              ggCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTAACGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGCTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGgc                                                                       >  1:246841/1‑149 (MQ=255)                                                                                 tttCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGcgcg                                                                    >  7:231/1‑149 (MQ=255)                                                                                     cAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGag                                                                >  3:65339/1‑148 (MQ=255)                                                                                           cTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGagcctttc                                                          <  4:65339/149‑8 (MQ=255)                                                                                             cGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCAc                                                        >  7:124076/1‑149 (MQ=255)                                                                                                  ccGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCt                                                   <  6:372775/149‑1 (MQ=255)                                                                                                     aTGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGAc                                                <  2:4275/149‑1 (MQ=255)                                                                                                        gggATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACaa                                              <  6:254143/148‑1 (MQ=255)                                                                                                             aCCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAg                                        >  4:152965/1‑149 (MQ=255)                                                                                                                   aGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTa                                  <  5:164135/149‑1 (MQ=255)                                                                                                                        gTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCTTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGgc                             >  5:239599/1‑149 (MQ=255)                                                                                                                             gTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCt                        >  7:295853/1‑149 (MQ=255)                                                                                                                                 aGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGc                     <  8:231/148‑1 (MQ=255)                                                                                                                                   cTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCgg                  <  3:152965/149‑1 (MQ=255)                                                                                                                                    tGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGc                 >  1:217196/1‑149 (MQ=255)                                                                                                                                      gcagcaGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGcc                <  8:182855/148‑1 (MQ=255)                                                                                                                                          cagAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCtctc           >  3:407514/1‑149 (MQ=255)                                                                                                                                           agagGCCAGCGTGGCAAGCGTTACAGCACGAGCGAGAGTACGAATGGTGTTCATAGGCGTGCTTCCCTGTGATCAAGTCAGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCg          <  4:206649/149‑1 (MQ=255)                                                                                                                                             ctgCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCgg         <  5:7188/146‑1 (MQ=255)                                                                                                                                             ctgCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCgg         <  3:66637/146‑1 (MQ=255)                                                                                                                                               gCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTa                                  >  3:162449/1‑121 (MQ=255)                                                                                                                                               gCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTa                                  <  4:162449/121‑1 (MQ=255)                                                                                                                                                    agTGGCAAGGGGTACAGCACGAGCGAGAGTACGATTGGCGTTCATAGGCGTGCTTCCCTCTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGCCGa  <  2:246841/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCTTCTTCTCGTCAGGGCCTTTCACCGTGACTTGTCGGGTGGCCAGGTCGATA‑GCC‑GTTACCTGTGTGGTGACATGGCTTTCCAGCGCCTCGCTCCCGATGGGGATACCGGCAGCCTGTACCGTGAAGCTGGCAGCAGAGGCCAGCGTGGCAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGA  >  NC_002947/4103056‑4103348

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