Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,507	A→T	6.9%	E308V (GAG→GTG) ‡	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,507	0	A	T	6.9%	198.5 / 3.6	69	E308V (GAG→GTG) ‡	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTC  >  NC_002947/4605359‑4605652                                                                                                                                                     |                                                                                                                                                  caAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGa                                                                                                                                                   <  7:227514/149‑1 (MQ=255)          ccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                            >  7:73969/1‑147 (MQ=255)              aGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTc                                                                                                                                      <  5:121169/149‑1 (MQ=255)                gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGa                                                                                                                                    <  8:357462/149‑1 (MQ=255)                gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGa                                                                                                                                    <  8:351833/149‑1 (MQ=255)                 tGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGat                                                                                                                                   <  3:392944/149‑2 (MQ=37)                   ggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCt                                                                                                                                 <  7:162571/149‑1 (MQ=255)                      aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGaa                                                                                                                              <  2:197238/149‑1 (MQ=255)                      aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTACACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCCGCCCAGGTCGAGCTg                                                                                                                                >  2:200062/1‑147 (MQ=255)                              ccATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAg                                                                                                                        >  1:329017/1‑147 (MQ=255)                                aTTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGc                                                                                                                    >  2:212263/1‑149 (MQ=255)                                 tttCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGc                                                                                                                    >  5:411064/1‑148 (MQ=255)                                 tttCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGc                                                                                                                    >  5:45885/1‑148 (MQ=255)                                     aCGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGTCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAgg                                                                                                                <  8:24931/148‑1 (MQ=255)                                     aCGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgc                                                                                                               <  6:45885/149‑1 (MQ=255)                                        cgTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCa                                                                                                            <  1:325764/149‑1 (MQ=255)                                          tCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGt                                                                                                          >  8:227316/1‑149 (MQ=255)                                             aCTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCg                                                                                                       <  7:227316/149‑1 (MQ=255)                                                 ccAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCtt                                                                                                    <  3:181357/148‑1 (MQ=255)                                                  cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                                  <  1:359833/149‑1 (MQ=255)                                                       acCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcg                                                                                                              <  3:200308/132‑1 (MQ=255)                                                       acCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcg                                                                                                              >  4:200308/1‑132 (MQ=255)                                                         ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCa                                                                                                                                          <  5:306486/102‑1 (MQ=255)                                                         ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCa                                                                                                                                          >  6:306486/1‑102 (MQ=255)                                                          cTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCTAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                          >  3:33015/1‑149 (MQ=255)                                                           ttCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                          >  8:7935/1‑148 (MQ=255)                                                             cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                        <  2:329017/148‑1 (MQ=255)                                                                gTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCCAGCTGAAGCCCAGAAGCGAGACGCAGTCCGCTTACGAAATCTTCAgtg                                                                                     >  4:207729/1‑148 (MQ=255)                                                                          aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGt                                                                           >  6:279765/1‑148 (MQ=255)                                                                              tttCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                       <  5:279765/148‑1 (MQ=255)                                                                                       gTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGa                                                              >  2:298859/1‑148 (MQ=255)                                                                                           tGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcgc                                                                                                             <  5:356373/97‑1 (MQ=255)                                                                                           tGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcgc                                                                                                             >  6:356373/1‑97 (MQ=255)                                                                                            gCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                         >  5:189002/1‑148 (MQ=255)                                                                                             cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                                <  1:369876/124‑1 (MQ=255)                                                                                             cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                                >  2:369876/1‑124 (MQ=255)                                                                                             cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                         >  3:353360/1‑147 (MQ=255)                                                                                             cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                         <  4:353360/147‑1 (MQ=255)                                                                                             cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaa                                                          <  6:125302/146‑1 (MQ=255)                                                                                             cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCa                                                       <  6:411064/149‑1 (MQ=255)                                                                                              aCCCCGTCTGTCATCTTGTTGAATCACAGTGCCGATCCCATCGACCTGGCCATCGTTGCTGACCATGTCGAGCTGAAGCCCAGAAGCGATGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAAccc                                                       <  3:207729/148‑2 (MQ=255)                                                                                                cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGaa                                                                             >  6:300861/1‑124 (MQ=255)                                                                                                cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                                     >  3:403644/1‑148 (MQ=255)                                                                                                cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                                     >  6:185505/1‑148 (MQ=255)                                                                                                        tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAg                                             <  4:403644/148‑1 (MQ=255)                                                                                                        tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCGGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGt                                            <  1:200062/149‑1 (MQ=255)                                                                                                          aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                          <  1:197129/101‑1 (MQ=255)                                                                                                          aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                          >  2:197129/1‑101 (MQ=255)                                                                                                          aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAc                                          >  8:265601/1‑149 (MQ=255)                                                                                                                 tCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCTAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGa                                    <  4:33015/148‑1 (MQ=255)                                                                                                                    aGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTATCGGTGAGCAc                                >  2:186455/1‑149 (MQ=255)                                                                                                                        caGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGc                            <  7:7935/149‑1 (MQ=255)                                                                                                                             ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGa                         <  2:249525/147‑1 (MQ=255)                                                                                                                             ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGa                         >  1:249525/1‑147 (MQ=255)                                                                                                                             ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGAc                        <  1:212263/148‑1 (MQ=255)                                                                                                                             ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACt                       >  6:331592/1‑149 (MQ=255)                                                                                                                                    cATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGcc                <  1:298859/149‑1 (MQ=255)                                                                                                                                     aTCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGccc               >  4:35358/1‑149 (MQ=255)                                                                                                                                      tCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACATGCGGACTTTTCGCCCt              >  5:344991/1‑149 (MQ=255)                                                                                                                                        gACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt            <  8:73969/149‑1 (MQ=255)                                                                                                                                         aCCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                       >  5:24922/1‑89 (MQ=37)                                                                                                                                         aCCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                       <  6:24922/89‑1 (MQ=37)                                                                                                                                         aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc           >  2:323977/1‑149 (MQ=255)                                                                                                                                         aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc           <  6:189002/149‑1 (MQ=255)                                                                                                                                          ccTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc           >  4:28314/1‑148 (MQ=38)                                                                                                                                              gCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAAt       <  5:185505/148‑1 (MQ=255)                                                                                                                                                cAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACtt                      <  8:227750/131‑1 (MQ=255)                                                                                                                                                cAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGt    >  6:394271/1‑149 (MQ=255)                                                                                                                                                 aGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGtt   >  2:7158/1‑149 (MQ=255)                                                                                                                                                  gCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTc  <  7:265601/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTC  >  NC_002947/4605359‑4605652

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