Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,294	T→G	7.8%	intergenic (+55/‑114)	eno → / → ftsB	enolase/cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,810,294	0	T	G	7.8%	194.0 / 7.7	64	intergenic (+55/‑114)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base G (2/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGC  >  NC_002947/1810146‑1810441                                                                                                                                                       |                                                                                                                                                        cGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACt                                                                                                                                                         >  2:177386/1‑149 (MQ=255)  gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                                        <  1:177386/149‑1 (MQ=255)  gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                                        >  6:121638/1‑149 (MQ=255)  gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                                        <  1:14215/149‑1 (MQ=255)  gACCGCGTTTCCAAGTACAAACAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGAAGCCGTAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                                        >  2:24180/1‑149 (MQ=255)    ccGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTc                                                                                                                                                      <  2:287888/149‑1 (MQ=255)      gcgTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCt                                                                                                                                                    <  1:191949/149‑1 (MQ=255)            ccAAGTACAACCCGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGGGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTa                                                                                                                                               <  6:204847/148‑1 (MQ=255)               aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTT‑CAGTAC‑GAACGTTCCTAcaactccg                                                                                                                                           <  6:246555/149‑8 (MQ=38)               aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTT‑CAGTAC‑GAACGTTCCTAcaactccg                                                                                                                                           <  6:252081/149‑8 (MQ=38)                      ccAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGt‑tt                                                                                                                                   >  2:60207/1‑149 (MQ=255)                      ccAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGt‑tt                                                                                                                                   <  2:332863/149‑1 (MQ=255)                           tgctgcGCATCGAAGAGGAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGGTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCC‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGt                                                                                                                               <  1:60207/149‑1 (MQ=255)                                       aaGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTgg                                                                                                                >  1:23778/1‑149 (MQ=255)                                          agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCt                                                                                                             >  5:379566/1‑149 (MQ=255)                                              aCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGc                                                                                                         <  5:93290/149‑1 (MQ=255)                                                tGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACgg                                                                                                                                        >  7:241744/1‑119 (MQ=255)                                                tGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACgg                                                                                                                                        <  8:241744/119‑1 (MQ=255)                                                  ggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGGTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAgg                                                                                                     >  1:334457/1‑149 (MQ=255)                                                    cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCGAGGCCAGGCt                                                                                                   <  5:121638/149‑1 (MQ=255)                                                     gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCACGGCTAAGCAAGTGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAgg                                                                                                     <  7:180687/146‑1 (MQ=255)                                                     gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCACGGCTAAGCAAGTGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAgg                                                                                                     >  8:180687/1‑146 (MQ=255)                                                       cAAAGCGGTTTACCGCTGTCGTGCCTAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCGAGGCCAGGCTTCg                                                                                                <  5:273526/149‑1 (MQ=255)                                                        aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGt                                                                                               <  6:379566/149‑1 (MQ=255)                                                               tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                                        >  2:37959/1‑149 (MQ=255)                                                               tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                                        >  2:43506/1‑149 (MQ=255)                                                               tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                                        >  5:150940/1‑149 (MQ=255)                                                                 tACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCg                                                                                       <  5:343471/148‑1 (MQ=255)                                                                   ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑Ct                                                                                                                  <  2:41648/119‑1 (MQ=255)                                                                   ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑Ct                                                                                                                  >  1:41648/1‑119 (MQ=255)                                                                    cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGATATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGTCTTCGTg                                                                                              <  7:185919/138‑1 (MQ=255)                                                                    cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCCGTCGACGAAG‑CTCTGGCCTCGGCCAGGCTTCGTGCTATTCGAgg                                                                                    <  1:37959/149‑1 (MQ=255)                                                                    cgcgGTCGTGCCGAGGTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGG‑TTCTGTCGACGAAGCC‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAgg                                                                                    <  1:43506/149‑1 (MQ=255)                                                                           gTGCCGAGTTTCGCGGCTAAGCAAGAGATCGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGAT‑A‑GTACTCAACTTTCCTACGTACAAACGCGT‑TTCTGTCGACGAAG‑C‑CTCGCCACGGCCAGGCATCGTGCTAATCGAGGCCCCGaac                                                                            >  6:364314/1‑148 (MQ=255)                                                                            tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                           <  8:182581/149‑1 (MQ=255)                                                                            tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                           <  2:334457/149‑1 (MQ=255)                                                                              ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTa                                                                          >  7:116276/1‑148 (MQ=255)                                                                                         ggCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                              <  8:116276/149‑1 (MQ=255)                                                                                           cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                              <  4:98909/147‑1 (MQ=255)                                                                                           cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                              >  3:98909/1‑147 (MQ=255)                                                                                               gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑t                                                       >  5:305889/1‑149 (MQ=255)                                                                                                      aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAAc                                                >  6:353495/1‑149 (MQ=255)                                                                                                      aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TGCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAAc                                                <  6:211271/149‑1 (MQ=255)                                                                                                         gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTTCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCt                                              <  2:312350/148‑1 (MQ=255)                                                                                                         gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGc                                                                                                         <  4:270199/90‑1 (MQ=255)                                                                                                         gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGc                                                                                                         >  3:270199/1‑90 (MQ=255)                                                                                                             aaaGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACGAACGGGT‑TTCTGTCGACGACG‑C‑CTGGGCGCGGCCAGGGTGCGTGGTATTCGAGGCCCCGAAGTAACAAGTACGCGGGCGCC‑GTTTTAACCTGGAt                                          >  5:195073/1‑148 (MQ=255)                                                                                                              aaGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATaa                                        <  6:150940/149‑1 (MQ=255)                                                                                                               aGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAAc                                       <  2:23778/149‑1 (MQ=255)                                                                                                                   agcagcCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGa                                   >  1:202343/1‑149 (MQ=255)                                                                                                                        ccGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCt                                              <  7:318790/133‑1 (MQ=255)                                                                                                                        ccGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTAGTGCTATTAGAGGCCCCGAAATAACAAGTAAGCGGCAGCC‑TTTTTAAACt                                              >  8:318790/1‑133 (MQ=255)                                                                                                                          ggAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTg                                             <  7:137692/132‑1 (MQ=255)                                                                                                                          ggAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CAGGCCACGGCCAGGCATAGTGGGATTAGAGGCCCCGAAATAAAAAATAACCGGCAACCATTTTTAAcag                                              >  8:137692/1‑130 (MQ=255)                                                                                                                             tacgTAGGAAAGTT‑CAGTAC‑GAACGTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTcc                         >  7:81272/4‑149 (MQ=255)                                                                                                                                 tAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAATTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                     >  6:282951/1‑149 (MQ=255)                                                                                                                                 tAGGAAAGTT‑CAGTAC‑GAACGTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                     >  5:110353/1‑149 (MQ=255)                                                                                                                                    gAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg                  <  6:305889/149‑1 (MQ=255)                                                                                                                                    gAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg                  >  7:209220/1‑149 (MQ=255)                                                                                                                                    gAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGACGCGCAGTCCCTATTgg                  >  2:244930/1‑149 (MQ=255)                                                                                                                                            c‑GTACTGAACTTTCATACTTACTAACTGGT‑TTCGGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTc            <  6:195073/148‑1 (MQ=255)                                                                                                                                               tACTGAACTTTCCTACGTACTAACGGGTGTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGAGCAGTCCCTATTGGTTGTTCCt          >  1:22048/1‑149 (MQ=255)                                                                                                                                                 cTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGTGCAGTCCCTATTGGTTGTTCCTCGt       >  7:37827/1‑149 (MQ=255)                                                                                                                                                  tGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGtt      >  3:140459/1‑149 (MQ=255)                                                                                                                                                  tGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGACGCGCAGTCCCTATTGGTTGTTCCTCGtt      <  1:244930/149‑1 (MQ=255)                                                                                                                                                  cgAACGTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGtt      <  6:110353/148‑1 (MQ=255)                                                                                                                                                    aaCTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg                  <  3:149620/135‑1 (MQ=255)                                                                                                                                                    aaCTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCt    >  3:169381/1‑149 (MQ=255)                                                                                                                                                      cTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGc  <  5:353495/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        CGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGT‑TTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGC  >  NC_002947/1810146‑1810441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑