Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,424,498	C→G	5.9%	L33V (CTG→GTG)	PP_4764 →	histone deacetylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,424,498	0	C	G	5.9%	209.0 / 4.1	68	L33V (CTG→GTG)	PP_4764	histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

CAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCT  >  NC_002947/5424351‑5424646                                                                                                                                                        |                                                                                                                                                       cAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCt                                                                                                                                                       >  7:203384/1‑149 (MQ=255) cAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCt                                                                                                                                                       >  7:176903/1‑149 (MQ=255)       cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCg‑g                                                                                                                                                <  8:176903/149‑1 (MQ=255)       cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCg‑g                                                                                                                                                <  8:203384/149‑1 (MQ=255)       cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCg‑g                                                                                                                                                <  4:192756/149‑1 (MQ=255)              tGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCggg                                                                                                                                         >  2:33969/1‑149 (MQ=255)                gCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCAGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCt                                                                                                                                       >  6:319198/1‑149 (MQ=255)                  aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGa                                                                                                                                     >  5:153944/1‑149 (MQ=255)                  aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGa                                                                                                                                     >  5:148437/1‑149 (MQ=255)                  aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGC‑CGACCAGGTGGTCgtgc                                                                                                                                               <  8:61397/139‑4 (MQ=255)                  aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGC‑CGACCAGGTGGTCgtgc                                                                                                                                               >  7:61397/1‑136 (MQ=255)                   aGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGac                                                                                                                                    >  7:28700/1‑149 (MQ=255)                         gCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGa                                                                                                                              >  2:178136/1‑149 (MQ=255)                            cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGTTCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGAc                                                                                                                             <  1:90684/147‑1 (MQ=255)                            cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGTTCTACCACGAAGATTACAGCCCG‑GAGTTACCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGAc                                                                                                                             >  2:90684/1‑147 (MQ=255)                            cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCa                                                                                                                           <  8:28700/149‑1 (MQ=255)                               gCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGc                                                                                                                        >  6:334697/1‑149 (MQ=255)                               gCCCTGATTGCGACCCCATCCATCCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GACTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGc                                                                                                                        <  4:311036/149‑1 (MQ=255)                                  cttaTTGCTACCCCATCCATGCCGTTGCCGCTCATCTACCACGAACATTCCCGCCCGTCCGTTCCCCGCTGACTCCAACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAgg                                                                                                                         <  1:112127/146‑1 (MQ=255)                                     aTTGCGACCCCATCCATGCCGTTGCAGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGACG‑GCAGCGCGCTGACCACCGCCCAGGCCctgct                                                                                                                  >  4:196308/1‑149 (MQ=255)                                      ttGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCctgct                                                                                                                  <  3:196308/148‑1 (MQ=255)                                        gCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCtgctgc                                                                                                                <  5:162245/148‑1 (MQ=255)                                           aCCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGccc                                                                                                            <  5:319198/149‑1 (MQ=255)                                            ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACAAGGCCCTGCTCCGCCCg                                                                                                           >  5:302023/1‑149 (MQ=255)                                              ccatccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGcc                                                                                                             >  5:352734/1‑145 (MQ=255)                                              ccatccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGa                                                                                                         >  4:255149/1‑149 (MQ=255)                                              ccatccatGCCGTTGCCGATGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGcc                                                                                                             <  6:352734/145‑1 (MQ=255)                                                 tccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACa                                                                                                       >  6:245810/1‑148 (MQ=255)                                                     tGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGc                                                                                                 <  5:334697/149‑1 (MQ=255)                                                     tGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGc                                                                                                 <  1:33969/149‑1 (MQ=255)                                                          ttGCCGCTGATCTACCACGAAGATTACAGTCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGaccacc                                                                                                                                >  1:210474/1‑114 (MQ=255)                                                          ttGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGaccacc                                                                                                                                <  2:210474/114‑1 (MQ=255)                                                             ccGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACg                                                                                          <  5:245810/148‑1 (MQ=255)                                                                     cTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCg                                                                                 >  7:323200/1‑149 (MQ=255)                                                                        ccACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGccc                                                                              >  7:215530/1‑149 (MQ=255)                                                                        ccACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGccc                                                                              >  2:89654/1‑149 (MQ=255)                                                                          aCGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGcc                                                                                                >  7:102807/1‑129 (MQ=255)                                                                          aCGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGcc                                                                                                <  8:102807/129‑1 (MQ=255)                                                                          aCGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGc                                                                                <  6:334432/145‑1 (MQ=255)                                                                          aCGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGc                                                                                >  5:334432/1‑145 (MQ=255)                                                                              aGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGccc                                                                        <  2:298543/149‑1 (MQ=255)                                                                               gATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCa                                                                       >  8:242317/1‑149 (MQ=255)                                                                                aTTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                                      <  8:323200/149‑1 (MQ=255)                                                                                    cAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGAcc                                                                  >  4:353898/1‑149 (MQ=255)                                                                                                 ccccGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGAc                                                                                        >  8:10414/1‑115 (MQ=255)                                                                                                 ccccGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGAc                                                                                        <  7:10414/115‑1 (MQ=255)                                                                                                 ccccGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATc                                                       >  2:159960/1‑148 (MQ=255)                                                                                                  cccGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCCGCTACATCgc                                                     >  3:183192/1‑148 (MQ=255)                                                                                                    cGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                   <  8:215530/149‑1 (MQ=255)                                                                                                    cGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                   <  7:242317/149‑1 (MQ=255)                                                                                                     gCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGGTGCGCCCGGACATTCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                                   <  1:89654/149‑1 (MQ=255)                                                                                                               acCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACAt                                            >  3:223276/1‑148 (MQ=255)                                                                                                                 cGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTgcgc                                                                                                              <  5:272290/81‑1 (MQ=255)                                                                                                                 cGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTgcgc                                                                                                              >  6:272290/1‑81 (MQ=255)                                                                                                                       ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                   <  1:178136/149‑1 (MQ=255)                                                                                                                              aCAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcg                            >  8:50790/1‑149 (MQ=255)                                                                                                                               cAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgc                           <  1:159960/149‑1 (MQ=255)                                                                                                                                aaGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                                      <  6:138347/105‑1 (MQ=255)                                                                                                                                aaGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcg                          >  3:80192/1‑149 (MQ=255)                                                                                                                                       gCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc                   <  4:80192/149‑1 (MQ=255)                                                                                                                                            gctgccGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                   <  7:7613/124‑1 (MQ=255)                                                                                                                                            gctgccGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                                   >  8:7613/6‑129 (MQ=255)                                                                                                                                            ctgcACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTTCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgc               <  6:153944/148‑1 (MQ=255)                                                                                                                                            ctgcACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAgc               <  6:148437/148‑1 (MQ=255)                                                                                                                                               cACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACgg                                      <  6:185845/122‑1 (MQ=255)                                                                                                                                               cACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACgg                                      >  5:185845/1‑122 (MQ=255)                                                                                                                                                 cGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCgccgcc         >  2:2601/1‑149 (MQ=255)                                                                                                                                                 cGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCgccgcc         >  1:182944/1‑149 (MQ=255)                                                                                                                                                       ccTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGcc   >  1:304681/1‑149 (MQ=255)                                                                                                                                                        cTGGTCG‑GCAGCGTGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCTAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCt  <  6:302023/149‑1 (MQ=255)                                                                                                                                                        cTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCt  >  8:113464/1‑149 (MQ=255)                                                                                                                                                        cTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCt  >  3:148641/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                       CAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCG‑GAGTTCCCCGCTGA‑‑AC‑ACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACA‑TCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCCGCCTCGGCCT  >  NC_002947/5424351‑5424646

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