Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	396,848	Δ1 bp	6.8%	coding (98/924 nt)	PP_0329 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	396,848	0	C	.	6.8%	286.2 / 13.4	88	coding (98/924 nt)	PP_0329	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base . (4/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TTTGTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCT  >  NC_002947/396700‑396991                                                                                                                                                      |                                                                                                                                                 tttGTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTc                                                                                                                                                  <  5:23744/149‑1 (MQ=255)     tCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAAt                                                                                                                                              >  6:157267/1‑149 (MQ=255)           tATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAg                                                                                                                                        <  7:173519/149‑1 (MQ=255)             tCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAgg                                                                                                                                       <  5:367319/148‑1 (MQ=255)               cAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                    <  1:229788/149‑1 (MQ=255)                aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                    >  5:425594/1‑148 (MQ=255)                aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                    >  8:392804/1‑148 (MQ=255)                aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                    >  3:111535/1‑148 (MQ=255)                     gggTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGc                                                                                                                             <  7:391373/149‑1 (MQ=255)                       gtgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTT‑CAATCCTGCGGATTGGAACGCtc                                                                                                                           <  1:176971/149‑2 (MQ=255)                        tgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTc                                                                                                                          <  6:269855/149‑1 (MQ=255)                                  aTTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGttt                                                                                                                 <  3:115837/148‑1 (MQ=255)                                    ttttttGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTcc                                                                                                               >  5:105844/1‑148 (MQ=255)                                      ttttGCTTTCCGTCAGGCGTTCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGTAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCaa                                                                                                             <  6:187997/148‑1 (MQ=255)                                       tttGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAAt                                                                                                            <  3:30405/148‑1 (MQ=255)                                         tGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGa                                                                                                          <  6:425594/148‑1 (MQ=255)                                         tGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGa                                                                                                          >  5:251835/1‑148 (MQ=255)                                          gCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGa                                                                                                          >  3:400383/1‑147 (MQ=255)                                          gCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCAAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGa                                                                                                          >  2:237210/1‑147 (MQ=255)                                           cTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGTCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGa                                                                                                          >  1:119751/1‑147 (MQ=255)                                           cTTTCCGTCAGGCGTGCTAACGACAAGCACTGAGTCCATCTCATCTGTTGCTAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCTTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACa                                                                                                        <  2:352560/148‑1 (MQ=255)                                            tttCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGt                                                                                                      >  2:438044/1‑149 (MQ=255)                                             ttCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTg                                                                                                     >  5:325310/1‑149 (MQ=255)                                              tCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGACCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTg                                                                                                     <  1:30587/148‑1 (MQ=255)                                                 gTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                 <  5:157267/149‑1 (MQ=255)                                                      gCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTg                                                                                            >  3:339231/1‑149 (MQ=255)                                                       cGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                           <  7:392804/149‑1 (MQ=255)                                                         tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCa                                                                                                                                         <  4:236207/102‑1 (MQ=255)                                                         tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCa                                                                                                                                         >  3:236207/1‑102 (MQ=255)                                                         tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                         <  2:119751/149‑1 (MQ=255)                                                            tCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTTGTCAGTTTGctggtc                                                                                      >  4:262728/1‑149 (MQ=255)                                                            tCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                      >  5:64892/1‑149 (MQ=255)                                                            tCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                      >  2:379175/1‑149 (MQ=255)                                                              aCGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctg                                                                                    <  6:251835/149‑1 (MQ=255)                                                                 aCAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctgg                                                                                   >  5:387002/1‑147 (MQ=255)                                                                       aCTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAata                                                                           >  1:435956/1‑149 (MQ=255)                                                                         tGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATAtt                                                                         >  1:153784/1‑149 (MQ=255)                                                                         tGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTCGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctggt                                                                                  <  5:176618/140‑1 (MQ=255)                                                                         tGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTCGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctggt                                                                                  >  6:176618/1‑140 (MQ=255)                                                                            gTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                      <  6:64892/149‑1 (MQ=255)                                                                               catctcatctGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                    <  1:379175/148‑1 (MQ=255)                                                                               catctcatctGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                    <  4:400383/148‑1 (MQ=255)                                                                                     atctGTTGCGAAAAG‑CTGGTAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCa                                                              <  2:322787/148‑1 (MQ=255)                                                                                     atctGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACa                                                                                                        >  8:249401/1‑106 (MQ=255)                                                                                     atctGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACa                                                                                                        <  7:249401/106‑1 (MQ=255)                                                                                     atctGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAg                                                             >  4:71318/1‑149 (MQ=255)                                                                                        tGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGAc                                                          >  8:144608/1‑149 (MQ=255)                                                                                         gTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACt                                                         >  1:407518/1‑149 (MQ=255)                                                                                          ttGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACt                                                         <  2:435956/148‑1 (MQ=255)                                                                                             cGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttt                                                     <  4:111535/149‑1 (MQ=255)                                                                                             cGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttt                                                     <  6:325310/149‑1 (MQ=255)                                                                                               aaaaG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttttt                                                   <  6:105844/149‑1 (MQ=255)                                                                                                aaaG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTc                                                  <  4:169878/149‑1 (MQ=255)                                                                                                      tGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGa                                             >  2:124149/1‑149 (MQ=255)                                                                                                       gggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                            >  6:223847/1‑149 (MQ=255)                                                                                                        ggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                            <  3:71318/148‑1 (MQ=255)                                                                                                         gAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAAACCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGt                                          >  1:102967/1‑149 (MQ=255)                                                                                                          aaCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTg                                         >  1:424483/1‑149 (MQ=255)                                                                                                             ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                      >  1:414529/1‑149 (MQ=255)                                                                                                             ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                      >  2:69636/1‑149 (MQ=255)                                                                                                               ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                      <  5:223847/147‑1 (MQ=255)                                                                                                                cTGCTTCTTTGCCAGCGGACTAAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCg                                   >  4:3466/1‑149 (MQ=255)                                                                                                                  gCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                 >  8:70548/1‑149 (MQ=255)                                                                                                                        ttGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                 >  1:268663/1‑79 (MQ=255)                                                                                                                        ttGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                 <  2:268663/79‑1 (MQ=255)                                                                                                                         tGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCt                             <  4:339231/146‑1 (MQ=255)                                                                                                                            ctgCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGtt                       >  6:430617/3‑149 (MQ=255)                                                                                                                             aGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTg                      <  6:387002/149‑1 (MQ=255)                                                                                                                              gCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                     <  1:438044/149‑1 (MQ=255)                                                                                                                                  aCTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                      <  6:195568/128‑1 (MQ=255)                                                                                                                                  aCTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                      >  5:195568/1‑128 (MQ=255)                                                                                                                                    tCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCgg               <  4:423302/149‑1 (MQ=255)                                                                                                                                      aGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCgg               <  3:262728/147‑1 (MQ=255)                                                                                                                                       gCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa            >  4:10894/1‑149 (MQ=255)                                                                                                                                        cAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGt                               >  5:368227/1‑129 (MQ=255)                                                                                                                                        cAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGt                               <  6:368227/129‑1 (MQ=255)                                                                                                                                        cAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAAt           <  2:407518/149‑1 (MQ=255)                                                                                                                                          agagGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAg         <  2:424483/149‑1 (MQ=255)                                                                                                                                           gagGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa            >  6:387164/1‑145 (MQ=255)                                                                                                                                             gggAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                 <  7:321488/122‑1 (MQ=255)                                                                                                                                             gggAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                 >  8:321488/1‑122 (MQ=255)                                                                                                                                             gggAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATg      <  2:342285/149‑1 (MQ=255)                                                                                                                                              gaaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg     >  7:64209/3‑149 (MQ=255)                                                                                                                                              gaaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg     >  3:324127/3‑149 (MQ=255)                                                                                                                                               aaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaa             >  4:86229/2‑140 (MQ=255)                                                                                                                                               aaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaa             <  3:86229/139‑1 (MQ=255)                                                                                                                                               aaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCGCCCGGTGTCTGCAGTTGGCCCCGGAAATAGATGGc    >  8:113789/2‑149 (MQ=255)                                                                                                                                                aGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGcc   >  2:26220/1‑149 (MQ=255)                                                                                                                                                 gCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCt  <  3:10894/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 TTTGTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAG‑CTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCT  >  NC_002947/396700‑396991

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