Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,446,635	G→C	9.4%	K343N (AAG→AAC)	gpD →	tail formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,446,635	0	G	C	9.4%	180.5 / 9.2	64	K343N (AAG→AAC)	gpD	tail formation
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (3/3);  total (34/30)
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.68e-01

GCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAA  >  NC_002947/3446494‑3446781                                                                                                                                              |                                                                                                                                                    gCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGagc                                                                                                                                              <  3:219366/149‑1 (MQ=255) gCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGagc                                                                                                                                              <  6:160359/149‑1 (MQ=255)  cTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGa                                                                                                                                                <  3:99200/146‑1 (MQ=255)  cTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGa                                                                                                                                                >  4:99200/1‑146 (MQ=255)  cTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTCag                                                                                                                                               >  1:137214/1‑147 (MQ=255)     gATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCt                                                                                                                                          <  8:28675/149‑1 (MQ=255)       tGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCTGCTCAACCTTGAGCGGCttg                                                                                                                                        <  7:213699/149‑3 (MQ=255)                ccTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACg                                                                                                                               >  2:294376/1‑149 (MQ=255)                 cTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACgg                                                                                                                              >  3:105089/1‑149 (MQ=255)                 cTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACgg                                                                                                                              <  1:186629/149‑1 (MQ=255)                     tggACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGt                                                                                                                          <  2:137214/149‑1 (MQ=255)                          ggggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAACGACCCCAAGCCGCTCAAGCTTGAGAAGCTGTAACCCCCAGTCCGTacacc                                                                                                                     >  8:295267/1‑148 (MQ=255)                            ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGc                                                                                                                   >  8:200796/1‑149 (MQ=255)                            ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGc                                                                                                                   >  8:251278/1‑149 (MQ=255)                                gAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCAt                                                                                                                >  1:100041/1‑148 (MQ=255)                                  gCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCAc                                                                                                             >  7:194595/1‑149 (MQ=255)                                  gCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCAc                                                                                                             >  4:282951/1‑149 (MQ=255)                                  gCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCAc                                                                                                             >  4:19085/1‑149 (MQ=255)                                    aGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTg                                                                                                           <  3:38464/149‑1 (MQ=255)                                      gTATTTACCCAGGCCGGCTGGACTACGACCGTCGAGTCCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCCGGCAAGAAAAACAAAGACCCCACGCCGCGCACGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGga                                                                                                         <  8:357450/149‑1 (MQ=255)                                        aTTTACACAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACACCACGGCCGTACACGGCCATCACTGgag                                                                                                        >  2:249108/1‑148 (MQ=255)                                              ccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCaa                                                                                                 >  3:80784/1‑149 (MQ=255)                                                     ggCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGaa                                                                                                      >  2:93703/1‑137 (MQ=255)                                                     ggCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGaa                                                                                                      <  1:93703/137‑1 (MQ=255)                                                     ggCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTAt                                                                                          <  3:19085/149‑1 (MQ=255)                                                                  ccGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGAAGCTGTAACCAAACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCAATATCAGTTCAACAg                                                                             >  5:236580/1‑149 (MQ=255)                                                                     tCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTg                                                                          >  2:1131/1‑149 (MQ=255)                                                                       gAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCCCGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGc                                                                         >  6:112140/1‑148 (MQ=255)                                                                              aCGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGAtc                                                                 <  1:294376/149‑1 (MQ=255)                                                                                 gcggcAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGAtcctc                                                              <  3:284329/149‑1 (MQ=255)                                                                                      aagaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGCCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑aa                                                        >  1:333273/1‑149 (MQ=255)                                                                                       agaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AAc                                                       <  4:105089/149‑1 (MQ=255)                                                                                        gaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACg                                                      >  6:219340/1‑149 (MQ=255)                                                                                        gaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACg                                                      >  6:224426/1‑149 (MQ=255)                                                                                          agGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACgcc                                                    >  1:298634/1‑149 (MQ=255)                                                                                             gCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCAACACTGGAGAAACCAATGGCTATCTCAGTTCACCAGCTGCAACACATCCTCCCC‑AACgccggc                                                 >  6:306483/1‑149 (MQ=255)                                                                                              caaggccaaggcCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACgccggcc                                                <  2:100041/149‑1 (MQ=255)                                                                                                aggccaaggccTCGGGCAAGAAAAAGAAAGAGACCAAGCCTCTCAAGGTTGAGCAGCTGTAACCCCACGGGCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGc                                              <  8:340209/149‑1 (MQ=255)                                                                                                   ccaaggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGt                                                                                    >  4:52708/1‑109 (MQ=255)                                                                                                   ccaaggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGt                                                                                    <  3:52708/109‑1 (MQ=255)                                                                                                          cTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCg                                    >  6:33466/1‑149 (MQ=255)                                                                                                          cTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGACCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCg                                    >  2:304212/1‑149 (MQ=255)                                                                                                             gggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCTCAAAGCCGGCGttt                                 <  5:112140/149‑1 (MQ=255)                                                                                                              ggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTa                                <  3:282951/149‑1 (MQ=255)                                                                                                                     aaaagaaaGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTccc                         >  1:21080/1‑149 (MQ=255)                                                                                                                      aacgaaagaCACCAAGCCGCTCAAGGATGAGCAGCTGTAACCCCACCGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCg                        <  8:356075/146‑1 (MQ=255)                                                                                                                       aagaaagaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCg                        <  1:249108/149‑1 (MQ=255)                                                                                                                        agaaagaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGa                      <  1:1131/149‑1 (MQ=255)                                                                                                                            agaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCa                  <  8:194595/149‑1 (MQ=255)                                                                                                                                ttacaGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCgg                                      <  6:178450/121‑1 (MQ=255)                                                                                                                                ttacaGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCgg                                      >  5:178450/5‑125 (MQ=255)                                                                                                                                  cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGcta            >  1:231678/1‑149 (MQ=255)                                                                                                                                  cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGcta            >  1:228836/1‑149 (MQ=255)                                                                                                                                  cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGct             >  5:225925/1‑148 (MQ=255)                                                                                                                                  cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGct             >  5:179803/1‑148 (MQ=255)                                                                                                                                    aGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTccc                         <  1:354775/134‑1 (MQ=255)                                                                                                                                    aGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTccc                         >  2:354775/1‑134 (MQ=255)                                                                                                                                      ctgCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTAt        >  3:166118/3‑149 (MQ=255)                                                                                                                                       cGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATg       <  2:298634/149‑1 (MQ=255)                                                                                                                                        gCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATgg      <  6:179803/149‑1 (MQ=255)                                                                                                                                         cTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACATGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATggg     <  6:300690/149‑1 (MQ=255)                                                                                                                                         cTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATggg     >  1:157653/1‑149 (MQ=255)                                                                                                                                            aaGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCGACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGaa  <  7:200796/149‑1 (MQ=255)                                                                                                                                            aaGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGaa  <  7:251278/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    GCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCC‑AACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAA  >  NC_002947/3446494‑3446781

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