Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,895	T→C	5.6%	K47E (AAG→GAG) ‡	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,895	0	T	C	5.6%	267.6 / 3.9	89	K47E (AAG→GAG) ‡	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base C (3/2);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGG  >  NC_002947/807747‑808041                                                                                                                                                      |                                                                                                                                                    tCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCtt                                                                                                                                                     <  6:107287/149‑1 (MQ=255)  cgcgGCCCTGGTCAGCAAGGTGCTCGACCAGGATGATGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGAAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTg                                                                                                                                                    <  7:90370/149‑1 (MQ=255)       cccTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAg                                                                                                                                               >  6:170569/1‑149 (MQ=255)          tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCCACAGGCGCTGCGCGGTCTTGATGACCTGCGCGGCGGTCATGTCCCGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGATCAGG‑Gt                                                                                                                                           >  7:165315/1‑149 (MQ=255)               aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                       >  6:4389/1‑148 (MQ=255)                gCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTa                                                                                                                                     >  2:310399/1‑149 (MQ=255)                  aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGTACAGG‑GTGCGGTACt                                                                                                                                   <  7:349652/149‑1 (MQ=255)                  aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGCTCATGTCCGGCACGATGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACt                                                                                                                                   >  6:366928/1‑149 (MQ=255)                    ggTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTgg                                                                                                                                 >  3:319049/1‑149 (MQ=255)                       gCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTcc                                                                                                                              >  1:151453/1‑149 (MQ=255)                        cTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                             <  4:61754/149‑1 (MQ=255)                         tCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACAACAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTg                                                                                                                            >  5:183270/1‑149 (MQ=255)                           gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                             >  8:252142/1‑146 (MQ=255)                           gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                             <  7:252142/146‑1 (MQ=255)                            aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTgg                                                                                                                           >  3:32476/1‑147 (MQ=255)                            aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                                  >  4:223967/1‑138 (MQ=25)                            aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                                  <  3:223967/140‑3 (MQ=25)                               aGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGtt                                                                                                                      <  4:32476/149‑1 (MQ=255)                                  aTGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                   <  2:116847/149‑1 (MQ=255)                                   tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                   >  2:336574/1‑148 (MQ=255)                                   tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGCCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                   >  5:194476/1‑148 (MQ=255)                                           aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTggg                                                                                                           <  4:319049/148‑1 (MQ=255)                                           aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                          >  1:355153/1‑149 (MQ=255)                                           aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                          <  5:170569/149‑1 (MQ=255)                                            aCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                          <  1:171383/148‑1 (MQ=255)                                               gCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTAc                                                                                                      <  6:265899/149‑1 (MQ=255)                                                cGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTa                                                                                                       >  3:142164/1‑147 (MQ=255)                                                      gcgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGTCAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCa                                                                                                <  1:310399/149‑1 (MQ=255)                                                       cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                              <  8:76081/149‑1 (MQ=255)                                                       cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                              <  2:151453/149‑1 (MQ=255)                                                          tGCGCGGTCTCGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTc                                                                                           >  3:253965/1‑149 (MQ=255)                                                            cgcgGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgc                                                                                         <  2:84123/149‑1 (MQ=255)                                                             ggggTGTTGTTGACCTGCGCGGGGGTCATGTCCGGCACGAAGTATTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgcg                                                                                        <  6:194476/147‑1 (MQ=255)                                                               ggTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCAGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGcc                                                                                      >  2:143636/1‑149 (MQ=255)                                                                   ttgttgACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                  >  7:10279/1‑149 (MQ=255)                                                                   gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                  >  2:221444/5‑149 (MQ=255)                                                                   gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                  >  4:93474/5‑149 (MQ=255)                                                                         aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                            <  1:336574/149‑1 (MQ=255)                                                                         aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTAG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGCAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                            >  6:289396/1‑149 (MQ=255)                                                                           cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                           >  3:35609/1‑148 (MQ=255)                                                                           cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                           >  8:67872/1‑148 (MQ=255)                                                                              cgcgGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAg                                                                        <  5:66872/148‑1 (MQ=255)                                                                                cggcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCt                                                                     >  8:270216/1‑149 (MQ=255)                                                                                  gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   <  2:138920/149‑1 (MQ=255)                                                                                  gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   <  7:67872/149‑1 (MQ=255)                                                                                   cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                                  <  4:253965/149‑1 (MQ=255)                                                                                   cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                                  <  4:51058/149‑1 (MQ=255)                                                                                      tCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGaa                                                                >  2:60983/1‑148 (MQ=255)                                                                                         tGTCCGGAACTAAGTAGTAGAAACCCCGGGTATCGCACACCCGTTCG‑AAGCCCCGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                            <  8:28121/149‑1 (MQ=255)                                                                                          gTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                            >  4:97630/1‑148 (MQ=255)                                                                                           tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   <  5:87628/140‑1 (MQ=255)                                                                                           tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   >  6:87628/1‑140 (MQ=255)                                                                                             cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAgg                                                                                                                                              >  5:26893/1‑64 (MQ=255)                                                                                             cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTcgc                                                        <  1:60983/149‑1 (MQ=255)                                                                                              ggCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTTACGTGGAAGCCTTTGAACAGCTcgcg                                                       <  6:183270/149‑1 (MQ=255)                                                                                                cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                      >  2:286983/1‑148 (MQ=255)                                                                                                cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                      >  7:284163/1‑148 (MQ=255)                                                                                                cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGa                                                     <  5:4389/149‑1 (MQ=255)                                                                                                   gAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                      >  2:141098/1‑145 (MQ=255)                                                                                                        agtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCgg                                              <  4:35609/148‑1 (MQ=255)                                                                                                         gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcg                                            >  7:220357/1‑149 (MQ=255)                                                                                                         gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgc                                             <  7:270216/148‑1 (MQ=255)                                                                                                          tagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                           >  1:149557/1‑149 (MQ=255)                                                                                                            gAAACCACGGGTATCGCACACCGGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTGGCGGAAGATCGgcgc                                           >  7:59424/1‑147 (MQ=255)                                                                                                            gAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                           >  6:33613/1‑147 (MQ=255)                                                                                                               aCCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                      >  2:342658/1‑149 (MQ=255)                                                                                                                ccACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCTGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTg                                     <  1:286983/149‑1 (MQ=255)                                                                                                                 cACGGGCATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGGCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                                    <  8:158044/149‑1 (MQ=255)                                                                                                                    cagTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                  >  1:56034/3‑148 (MQ=25)                                                                                                                     agTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                  >  3:95580/2‑147 (MQ=25)                                                                                                                       tATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGt                              >  5:21947/1‑149 (MQ=255)                                                                                                                         tCGCACACCATAGCG‑TACCCCAGGGCCTTGCACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACACCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGtt                             <  8:165315/148‑1 (MQ=255)                                                                                                                         tCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGTCAGtt                             <  8:284163/148‑1 (MQ=255)                                                                                                                           gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                            <  1:82327/147‑1 (MQ=255)                                                                                                                           gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                            >  2:82327/1‑147 (MQ=255)                                                                                                                           gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                          <  1:141098/149‑1 (MQ=255)                                                                                                                            cacaCCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                         <  3:97630/149‑1 (MQ=255)                                                                                                                               ccctgTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaaa                       <  4:95580/145‑1 (MQ=37)                                                                                                                               acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGTAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAg                      <  8:220357/149‑1 (MQ=255)                                                                                                                               acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                         >  3:37701/1‑146 (MQ=255)                                                                                                                                ccAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGa                     <  2:355153/149‑1 (MQ=255)                                                                                                                                 cAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTCAACAGCTCGCGGAAGATCGGCGCGAGCTGGTACAGTTCGGAAAGa                     >  4:49877/1‑148 (MQ=255)                                                                                                                                    ttCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCaa                 <  5:33613/149‑1 (MQ=255)                                                                                                                                       g‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCAGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGac              >  3:128662/1‑148 (MQ=255)                                                                                                                                         aaGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATg             >  8:50674/1‑149 (MQ=255)                                                                                                                                             ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg                <  7:190377/142‑1 (MQ=255)                                                                                                                                             ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg                >  8:190377/1‑142 (MQ=255)                                                                                                                                                    cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGTGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg  >  2:28929/1‑149 (MQ=255)                                                                                                                                                    cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg  <  4:37701/149‑1 (MQ=255)                                                                                                                                                    cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg  >  4:101239/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGG  >  NC_002947/807747‑808041

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