Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,825,791	A→G	6.2%	F1800S (TTC→TCC)	pvdL ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,825,791	0	A	G	6.2%	189.2 / 2.3	64	F1800S (TTC→TCC)	pvdL	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GTCAGCACGAAGAAGTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTGCAGGGTGCC  >  NC_002947/4825644‑4825939                                                                                                                                                     |                                                                                                                                                     gTCAGCACGAAGAAGTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGaa                                                                                                                                                     <  5:10234/149‑1 (MQ=255)    aGCACGAAGAAGTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTc                                                                                                                                                  <  8:296784/149‑1 (MQ=255)       aCGAAGAAGTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTg                                                                                                                                               >  7:140779/1‑149 (MQ=255)         gaagaaGTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGGCGGAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGc                                                                                                                                              >  5:262871/1‑148 (MQ=255)           agaagTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCCGACAAAGC‑GCTGAAGTCTTGCCAg                                                                                                                                           >  5:299606/1‑149 (MQ=255)                  cTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGAt                                                                                                                                    <  2:266619/149‑1 (MQ=255)                    cgcgCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATg                                                                                                                                   >  8:66402/1‑148 (MQ=255)                     gcgcTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATg                                                                                                                                   >  1:1011/1‑147 (MQ=255)                           gTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGc                                                                                                                           >  4:289839/1‑149 (MQ=255)                            ttGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCt                                                                                                                          <  5:335833/149‑1 (MQ=255)                            ttGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCt                                                                                                                          <  2:199083/149‑1 (MQ=255)                                ccTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCt                                                                                                                      <  2:1011/149‑1 (MQ=255)                                     aCCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCgg                                                                                                                  <  1:167906/148‑1 (MQ=255)                                       caggcaggCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCa                                                                                                               <  8:140779/149‑1 (MQ=255)                                        aggcaggcGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCAc                                                                                                              >  5:13727/1‑149 (MQ=255)                                              gcgcgcAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTg                                                                                                        >  7:11830/1‑149 (MQ=255)                                               cgcgcAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTgg                                                                                                       <  6:347947/149‑1 (MQ=255)                                                 cgcAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTa                                                                                                     <  8:11830/149‑1 (MQ=255)                                                    agcagcGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACgg                                                                                                  >  6:221877/1‑149 (MQ=255)                                                      cagcGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGaa                                                                                                                                                     >  3:320635/1‑96 (MQ=255)                                                      cagcGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGaa                                                                                                                                                     <  4:320635/96‑1 (MQ=255)                                                            cccccGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGt                                                                                          >  2:19425/1‑149 (MQ=255)                                                                     cAGGTCGGAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCCGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTg                                                                                  <  6:262871/148‑1 (MQ=255)                                                                     cAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCCGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTg                                                                                  <  6:299606/148‑1 (MQ=255)                                                                      aGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTggg                                                                                <  1:253004/149‑1 (MQ=255)                                                                           gAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACg                                                                           >  4:334873/1‑149 (MQ=255)                                                                                 cTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTg                                                                                  >  6:32246/1‑136 (MQ=255)                                                                                 cTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTg                                                                                  <  5:32246/136‑1 (MQ=255)                                                                                 cTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgc                                                                     <  5:249705/149‑1 (MQ=255)                                                                                 cTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgc                                                                     <  5:144654/149‑1 (MQ=255)                                                                                    gTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgcgc                                                                   >  5:172067/1‑148 (MQ=255)                                                                                    gTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgcgc                                                                   >  3:303247/1‑148 (MQ=255)                                                                                    gTGCGCCTGGTCGTAGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgcgc                                                                   >  7:310356/1‑148 (MQ=255)                                                                                     tGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgc                                                                     >  3:106862/1‑145 (MQ=255)                                                                                     tGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTgc                                                                     <  4:106862/145‑1 (MQ=255)                                                                                     tGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCaa                                                                 <  5:221877/149‑1 (MQ=255)                                                                                         ccTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGtt                                                             >  4:222048/1‑149 (MQ=255)                                                                                              tcgtcgGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGCGGATGAAGTATTGCCAGTCCAGATGCAGGCCGATTTCTTCGGCCACGCACTGTTACGGAACGCCGTCCATGCTGGGGAACTTGTTGCGCAAGGTTTctt                                                         >  1:88614/1‑147 (MQ=255)                                                                                                   ggCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCa                                                   >  2:333389/1‑149 (MQ=255)                                                                                                    gCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCAc                                                  >  1:307373/1‑149 (MQ=255)                                                                                                    gCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCAc                                                  >  2:121850/1‑149 (MQ=255)                                                                                                       aGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGAt                                               <  7:66402/149‑1 (MQ=255)                                                                                                         ttgttgCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGAt                                               >  7:8856/1‑147 (MQ=255)                                                                                                           gttgCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAgc                                           >  7:279142/1‑149 (MQ=255)                                                                                                              gCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTg                                                                                  <  8:198666/107‑1 (MQ=255)                                                                                                              gCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTg                                                                                  >  7:198666/1‑107 (MQ=255)                                                                                                                aaaCGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTg                                      <  4:303247/149‑1 (MQ=255)                                                                                                                       gCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAgcgc                                <  2:88614/148‑1 (MQ=255)                                                                                                                          ggcaggcaTCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGc                            <  3:289839/149‑1 (MQ=255)                                                                                                                               gcaTCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGaa                       <  1:121850/149‑1 (MQ=255)                                                                                                                                caTCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGt                                                                                                      <  7:46515/69‑1 (MQ=255)                                                                                                                                caTCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGt                                                                                                      >  8:46515/1‑69 (MQ=255)                                                                                                                                caTCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGGACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAg                      <  2:307373/149‑1 (MQ=255)                                                                                                                                caTCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGGCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAg                      >  8:128291/1‑149 (MQ=255)                                                                                                                                   cGTCTGACAAAGC‑GCTGGAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGc                                     >  1:150417/1‑131 (MQ=255)                                                                                                                                   cGTCTGACAAAGC‑GCTGGAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGc                                     <  6:258666/131‑1 (MQ=255)                                                                                                                                   cGTCTGACAAAGC‑GCTGGAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGc                                     >  5:258666/1‑131 (MQ=255)                                                                                                                                   cGTCTGACAAAGC‑GCTGGAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGc                                     <  2:150417/131‑1 (MQ=255)                                                                                                                                    gTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGt                  <  6:172067/149‑1 (MQ=255)                                                                                                                                    gTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGt                  <  8:310356/149‑1 (MQ=255)                                                                                                                                     tCTGCCAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTc                 >  2:236793/1‑149 (MQ=255)                                                                                                                                        gACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTc                 >  6:297470/1‑146 (MQ=255)                                                                                                                                        gACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTc                 <  5:297470/146‑1 (MQ=255)                                                                                                                                        gACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCAc              >  5:12124/1‑149 (MQ=255)                                                                                                                                           aaaGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTg           >  8:171652/1‑149 (MQ=255)                                                                                                                                              gc‑gcTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTGCAg        <  7:128291/149‑1 (MQ=255)                                                                                                                                               c‑gcTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTGCAgg       <  8:279142/149‑1 (MQ=255)                                                                                                                                                     aaGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTGCAGGGTGcc  >  7:9082/1‑149 (MQ=255)                                                                                                                                                     aaGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTGCAGGGTGc   <  7:171652/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GTCAGCACGAAGAAGTGCTCGCGCTCGTTGGCCTTGACCAGGCAGGCGCGCAGCAGCGGCCCCCGCTCCAGGTCGAAGGGCTGGTGCGCCTGGTCGTCGGCCAGTTGTTGCAAACGTTGCTGGCAGGCATCGTCTGACAAAGC‑GCTGAAGTCTTGCCAGTCCAGATGCAGGCCGCTGTCTTCGGCCACGCACTGGTACGGCACGCCGTCGATGCTGGGGAACGTGGTGCGCAAGGTTTCGTGGCGCACGATCAGCGCCTGCAGCGCACGCTCGAAGGCGTCCACGTGCAGGGTGCC  >  NC_002947/4825644‑4825939

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