Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,774,751	A→T	5.1%	L254Q (CTG→CAG)	PP_5521 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,774,751	0	A	T	5.1%	241.5 / 2.7	78	L254Q (CTG→CAG)	PP_5521	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base T (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GCAAAGTGTGGGTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGG  >  NC_002947/2774607‑2774896                                                                                                                                                 |                                                                                                                                                  gCAAAGTGTGGGTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCt                                                                                                                                               <  6:165701/149‑1 (MQ=255)    aaGTGTGGGTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGAt                                                                                                                                            <  8:302428/149‑1 (MQ=255)         tgGGTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATgtcgt                                                                                                                                       >  5:366620/1‑149 (MQ=255)           ggTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTcc                                                                                                                                     >  5:62878/1‑149 (MQ=255)           ggTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTcc                                                                                                                                     >  1:197119/1‑149 (MQ=255)            gTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCg                                                                                                                                    >  1:363517/1‑149 (MQ=255)             tATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGGCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCgg                                                                                                                                   <  7:347695/149‑1 (MQ=255)              aTTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGTTCCATGCTCAGCCCCAGGCTGATGTCGTCCggg                                                                                                                                  >  1:170504/1‑149 (MQ=255)                  aCTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCTGGCTGATGTCGTCCGGGTCgg                                                                                                                              >  4:10309/1‑149 (MQ=255)                  aCTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCTGGCTGATGTCGTCCGGGTCgg                                                                                                                              >  6:397238/1‑149 (MQ=255)                  aCTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCgg                                                                                                                              >  3:153207/1‑149 (MQ=255)                  aCTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCgg                                                                                                                              <  1:405098/149‑1 (MQ=255)                   cTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGt                                                                                                                             >  1:221569/1‑149 (MQ=255)                   cTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCACGCTGATGTCGTCCGGGTCGGt                                                                                                                             <  4:113036/149‑1 (MQ=255)                      cccGGGCCGATGTGCGCCAGCGGCGTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGTTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTg                                                                                                                            <  2:170504/147‑1 (MQ=255)                       ccGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAACCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGccc                                                                                                                         >  5:266828/1‑149 (MQ=255)                         gggCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTc                                                                                                                       <  8:41776/149‑1 (MQ=255)                         gggCCGATGTGCGCCAGCGGCTT‑‑GTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCgg                                                                                                                     <  6:376475/149‑1 (MQ=255)                          ggCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCg                                                                                                                      >  7:82189/1‑149 (MQ=255)                           gCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCgg                                                                                                                     <  2:363517/149‑1 (MQ=255)                             cGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGcc                                                                                                                   >  5:110773/1‑149 (MQ=255)                               aTGTGCGCCAGCGGCTTGAGTCGTTCCTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAg                                                                                                                 <  6:169799/149‑1 (MQ=255)                                       cAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTc                                                                                                           >  2:147543/1‑147 (MQ=255)                                       cAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTc                                                                                                           >  2:141734/1‑147 (MQ=255)                                             cTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCTGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCAcgc                                                                                                   <  5:397238/149‑1 (MQ=255)                                             cTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCTGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCAcgc                                                                                                   <  3:10309/149‑1 (MQ=255)                                             cTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCCGCTCATCGAGCAcgc                                                                                                   >  7:255888/1‑149 (MQ=255)                                                gAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGt                                                                                                 <  6:366620/148‑1 (MQ=255)                                                 aGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGc                                                                                               >  5:268354/1‑149 (MQ=255)                                                 aGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGc                                                                                               <  4:153207/149‑1 (MQ=255)                                                        cgtCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCctgctg                                                                                         <  6:105089/148‑1 (MQ=255)                                                          tCCATAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGGCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTgg                                                                                      <  3:871/149‑1 (MQ=255)                                                          tCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTgg                                                                                      <  5:7656/149‑1 (MQ=255)                                                          tCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTgg                                                                                      >  1:326140/1‑149 (MQ=255)                                                           ccAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCtt                                                                                                            <  6:243495/126‑1 (MQ=255)                                                           ccAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCtt                                                                                                            >  5:243495/1‑126 (MQ=255)                                                            cAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTgg                                                                                      >  5:302070/1‑147 (MQ=255)                                                            cAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTgg                                                                                      >  2:262403/1‑147 (MQ=255)                                                            cAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGgca                                                                                    <  3:360725/149‑1 (MQ=255)                                                            cAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTCCCTGCTGTgg                                                                                      >  4:360725/1‑147 (MQ=255)                                                              gAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGgcag                                                                                   <  5:152005/148‑1 (MQ=255)                                                              gAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGgcag                                                                                   >  6:152005/1‑148 (MQ=255)                                                               aaCAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCctgctgccacag                                                                                   >  8:245767/1‑141 (MQ=255)                                                                 cAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAgg                                                                               <  2:197119/149‑1 (MQ=255)                                                                    gATGTCGGCTGCCTTGCTCACCCACGACGGTCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGc                                                                              >  7:344569/1‑147 (MQ=255)                                                                    gATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGCCACAGCAGGcac                                                                            <  7:245767/149‑3 (MQ=255)                                                                    gATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCCCAGGTTGCGCCCTGGTGAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGa                                                                            <  1:338104/149‑1 (MQ=255)                                                                      tGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGt                                                                          <  1:141734/149‑1 (MQ=255)                                                                        tCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTc                                                                         <  7:17496/148‑1 (MQ=255)                                                                        tCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCg                                                                        >  4:200650/1‑149 (MQ=255)                                                                          ggCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGt                                                                          <  1:147543/145‑1 (MQ=255)                                                                            cTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCAcc                                                                    <  6:110773/149‑1 (MQ=255)                                                                                cTTTCTCACCCACGACGATCAGCCTGCCATCCGCCGACACCGCACCATTCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGa                                                                <  8:255888/149‑1 (MQ=255)                                                                                   gCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCAcgc                                                                                                   <  3:163455/111‑1 (MQ=255)                                                                                   gCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCAcgc                                                                                                   >  4:163455/1‑111 (MQ=255)                                                                                    cTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTg                                                             <  8:82189/148‑1 (MQ=255)                                                                                    cTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTg                                                             <  1:262403/148‑1 (MQ=255)                                                                                       aCCCACGACGGTCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGc                                                           <  8:344569/147‑1 (MQ=255)                                                                                           acgacgATCAGCCGGTCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAc                                                     >  2:64707/1‑149 (MQ=255)                                                                                           acgacgATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAc                                                     >  4:181427/1‑149 (MQ=255)                                                                                               cgaTCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAg                                                 <  6:268354/149‑1 (MQ=255)                                                                                                 aTCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGa                                                >  3:202777/1‑148 (MQ=255)                                                                                                     gCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTa                                                      >  7:372476/1‑138 (MQ=255)                                                                                                     gCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTa                                                      <  8:372476/138‑1 (MQ=255)                                                                                                                gccgACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGc                                                           >  2:350194/1‑122 (MQ=255)                                                                                                                  cgACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGTGTGCct                                                                                             <  3:324478/86‑1 (MQ=255)                                                                                                                  cgACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGTGTGCct                                                                                             >  4:324478/1‑86 (MQ=255)                                                                                                                  cgACACCGCACCATGCGCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACca                              <  6:62878/149‑1 (MQ=255)                                                                                                                   gACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACcag                             >  5:242156/1‑149 (MQ=255)                                                                                                                   gACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACcag                             <  3:181427/149‑1 (MQ=255)                                                                                                                    acacCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACcag                             >  4:221623/1‑148 (MQ=255)                                                                                                                      acCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACcagcag                          >  6:103475/1‑149 (MQ=255)                                                                                                                             catgctccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGa                                                >  8:208066/1‑120 (MQ=255)                                                                                                                             catgctccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGa                                                <  7:208066/120‑1 (MQ=255)                                                                                                                                  tccatgctcAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGcc              >  4:53166/1‑149 (MQ=255)                                                                                                                                             ccccAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAgg   <  2:326140/149‑1 (MQ=255)                                                                                                                                             ccccAGGCTGATGTCGTCCGGGGCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAgg   <  4:218087/149‑1 (MQ=255)                                                                                                                                              cccATGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAggg  <  2:221569/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GCAAAGTGTGGGTATTCACTGCCCGGGCCGATGTGCGCCAGCGGCTTGAGTCGTTCGTCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCGGTGCCCTCGGCCAGCTCTTCGAGCACGCGTGCCTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGG  >  NC_002947/2774607‑2774896

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