Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,623,395	G→T	7.3%	R115R (CGG→AGG)	PP_4941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,623,395	0	G	T	7.3%	300.2 / 9.0	94	R115R (CGG→AGG)	PP_4941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/41);  new base T (4/3);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

TGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCGTAGGCA  >  NC_002947/5623250‑5623543                                                                                                                                                   |                                                                                                                                                       cgcggCCGGCATTGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAgg                                                                                                                                                     >  6:396647/5‑149 (MQ=255)    aGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTc                                                                                                                                                  >  8:24765/1‑149 (MQ=255)     gCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTc                                                                                                                                                  >  2:382685/1‑148 (MQ=255)      ccGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga                                                                                                                                                <  5:378133/149‑1 (MQ=255)       cGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGCCGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTc                                                                                                                                                  >  1:304301/1‑147 (MQ=255)          cAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCgaga                                                                                                                                              >  5:38592/1‑147 (MQ=255)          cAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGGCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGGCGAGaaa                                                                                                                                            <  6:115636/149‑1 (MQ=255)               ccggccgCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGcc                                                                                                                                       >  2:212402/1‑149 (MQ=255)                cggccgCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGTAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGccc                                                                                                                                      >  2:5300/1‑149 (MQ=255)                cggccgCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGccc                                                                                                                                      >  8:76065/1‑149 (MQ=255)                              cAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCttt                                                                                                                       <  3:244823/149‑1 (MQ=255)                              cAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCttt                                                                                                                       <  8:199300/149‑1 (MQ=255)                                aCGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCttttt                                                                                                                     <  6:38592/149‑1 (MQ=255)                                 cGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTa                                                                                                                    <  1:365761/149‑1 (MQ=255)                                 cGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTa                                                                                                                    <  5:390458/149‑1 (MQ=255)                                   tACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGc                                                                                                                  >  1:327687/1‑149 (MQ=255)                                   tACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGc                                                                                                                  <  2:327687/149‑1 (MQ=255)                                          cTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGtt                                                                                                           >  5:359658/1‑149 (MQ=255)                                             cAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑aa                                                                                                       <  5:5002/149‑1 (MQ=255)                                              aGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAg                                                                                                                              >  1:276206/1‑126 (MQ=255)                                              aGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAg                                                                                                                              <  2:276206/126‑1 (MQ=255)                                                tGGCGCAGC‑AAGACCCGTCAGGTTAAGGTAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAACCCAGCTGCTTTTTAGCCGCAGTTC‑AAGc                                                                                                     <  1:5300/149‑1 (MQ=255)                                                  gcgcAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCAc                                                                                                   <  7:76065/148‑1 (MQ=255)                                                    gcagc‑aAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAg                                                                                                >  2:71837/1‑149 (MQ=255)                                                    gcagc‑aAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAg                                                                                                >  2:243522/1‑149 (MQ=255)                                                        c‑aAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGtt                                                                                            >  5:247143/1‑149 (MQ=255)                                                            gACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTtcg                                                                                         >  5:156601/1‑149 (MQ=255)                                                               ccGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTtcgtcg                                                                                      >  2:299226/1‑149 (MQ=255)                                                               ccGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTtcgtcg                                                                                      >  5:160451/1‑149 (MQ=255)                                                                cGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGc                                                                                     <  5:396647/149‑1 (MQ=255)                                                                 gTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGcc                                                                                    >  3:186040/1‑149 (MQ=255)                                                                 gTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGcc                                                                                    >  1:253984/1‑149 (MQ=255)                                                                 gTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGcc                                                                                    >  5:384409/1‑149 (MQ=255)                                                                  tCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCt                                                                                   <  1:165504/149‑1 (MQ=255)                                                                  tCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCt                                                                                   <  3:269757/149‑1 (MQ=255)                                                                     ggTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCtt                                                                                  <  8:42715/147‑1 (MQ=255)                                                                     ggTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCtt                                                                                  >  7:42715/1‑147 (MQ=255)                                                                     ggTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTg                                                                                 >  6:163554/1‑148 (MQ=255)                                                                         aaGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAAc                                                                             <  5:163554/148‑1 (MQ=255)                                                                         aaGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCGGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAAc                                                                             <  6:332940/148‑1 (MQ=255)                                                                          aGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGt                                                                           <  2:98001/149‑1 (MQ=255)                                                                           ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTTGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTtc                                                                                          <  1:106049/133‑1 (MQ=255)                                                                           ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCgaga                                                                                                                                              >  3:240530/1‑83 (MQ=255)                                                                           ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCgaga                                                                                                                                              <  4:240530/83‑1 (MQ=255)                                                                           ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTtc                                                                                          >  2:106049/1‑133 (MQ=255)                                                                           ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTg                                                                          <  4:312714/149‑1 (MQ=255)                                                                           ggAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTGCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTGC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTg                                                                          >  5:129179/1‑149 (MQ=255)                                                                                 gcACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAAt                                                                      >  1:47081/1‑147 (MQ=255)                                                                                  cACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAg                                                                   >  3:113141/1‑149 (MQ=255)                                                                                     gCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCt                                                                <  6:359658/149‑1 (MQ=255)                                                                                       cATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgagct                                                                                                                                             <  6:99408/72‑3 (MQ=39)                                                                                       cATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgagct                                                                                                                                             >  5:99408/1‑70 (MQ=255)                                                                                         tCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGag                                                            <  2:253984/149‑1 (MQ=255)                                                                                         tCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGACCAAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCACAGTGCAATCAGCCTTgcg                                                            >  6:338909/1‑147 (MQ=255)                                                                                            aCGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagcca                                                         <  6:384409/149‑1 (MQ=255)                                                                                             cGCCAAGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCATCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTCGAAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagcca                                                         <  1:299226/149‑1 (MQ=255)                                                                                              gCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagccagc                                                       >  3:379152/1‑149 (MQ=255)                                                                                                 aCGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCATCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCATTGTTTCTTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCAt                                                    <  5:338909/149‑1 (MQ=255)                                                                                                       gAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTg                                                >  8:23953/1‑147 (MQ=255)                                                                                                       gAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTg                                                >  4:368779/1‑147 (MQ=255)                                                                                                       gAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTg                                                >  4:84739/1‑147 (MQ=255)                                                                                                       gAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCATCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGc                                               <  6:247143/148‑1 (MQ=255)                                                                                                           gCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCt                                          >  5:341272/1‑149 (MQ=255)                                                                                                           gCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCt                                          >  1:38793/1‑149 (MQ=255)                                                                                                             ttGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTgg                                        >  8:257980/1‑149 (MQ=255)                                                                                                             ttGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTgg                                        >  2:216232/1‑149 (MQ=255)                                                                                                             ttGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTgg                                        >  2:220969/1‑149 (MQ=255)                                                                                                               gCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTgg                                        <  1:212402/147‑1 (MQ=255)                                                                                                                ccGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCa                                      <  1:383209/148‑1 (MQ=255)                                                                                                                 cGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTA‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAg                                     >  3:372648/1‑148 (MQ=255)                                                                                                                    tCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAAc                                 <  1:71837/149‑1 (MQ=255)                                                                                                                      ggCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCCCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACct                               <  2:304301/149‑1 (MQ=255)                                                                                                                        ccTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACctgc                             <  2:132363/149‑1 (MQ=255)                                                                                                                         cTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACctgct                            >  6:94621/1‑149 (MQ=255)                                                                                                                              cGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCa                                      >  6:383813/1‑134 (MQ=255)                                                                                                                              cGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCa                                      <  5:383813/134‑1 (MQ=255)                                                                                                                              cGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCAc                        <  5:94621/148‑1 (MQ=255)                                                                                                                                 cagccagcAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGa                    <  4:186040/149‑1 (MQ=255)                                                                                                                                   gccagcAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATc                  <  1:216232/149‑1 (MQ=255)                                                                                                                                        tttcTCGACCCTCGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGcc               >  8:377948/4‑147 (MQ=255)                                                                                                                                        tttcTCGACCCTCGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGcc               <  7:377948/144‑1 (MQ=255)                                                                                                                                          gcTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAg           <  3:368779/149‑1 (MQ=255)                                                                                                                                           cTCGACCCTCGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCg                                                  <  5:13116/109‑1 (MQ=255)                                                                                                                                           cTCGACCCTCGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCg                                                  >  6:13116/1‑109 (MQ=255)                                                                                                                                           cTCGACCCTCGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGt          >  6:270754/1‑149 (MQ=255)                                                                                                                                           cTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagcca                                                         <  2:64950/102‑1 (MQ=255)                                                                                                                                           cTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGagcca                                                         >  1:64950/1‑102 (MQ=255)                                                                                                                                           cTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGt          >  3:362409/1‑149 (MQ=255)                                                                                                                                              gACCCGAGGTTCGAGTAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCTTATCCTTGCAAAGTTCAATCAGCATTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAATTCGt       <  2:283462/149‑1 (MQ=255)                                                                                                                                              gACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCg        >  1:5877/1‑148 (MQ=255)                                                                                                                                               aCCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCg        >  4:46129/1‑147 (MQ=255)                                                                                                                                               aCCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCGTa      <  7:24765/149‑1 (MQ=255)                                                                                                                                                cccGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCGTAg     >  8:192415/1‑149 (MQ=255)                                                                                                                                                 ccGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCGTAgg    >  1:13014/1‑149 (MQ=255)                                                                                                                                                   gAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCGTAGGCa  >  5:388771/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       TGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGC‑AAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCGCCCCAA‑CCAGCTGCTTTTTAGCCGCAGTTC‑AAGCACCAGGGTTTCGTCGCCTTGCAACGTGCAATCAGCCTTGAGCCAGCCATCGTGCAGGCTGGCAGGAACCTGCTGCACCTGATCGCCGAAGTCGTAGGCA  >  NC_002947/5623250‑5623543

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