Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,280	C→A	6.3%	intergenic (+41/‑128)	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,280	0	C	A	6.3%	203.7 / 3.6	63	intergenic (+41/‑128)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base A (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

CCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGG  >  NC_002947/1810140‑1810425                                                                                                                                             |                                                                                                                                                       ccGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTAc                                                                                                                                              >  5:204847/1‑149 (MQ=255)       cGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACt                                                                                                                                        >  2:177386/1‑149 (MQ=255)        gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                       <  1:177386/149‑1 (MQ=255)        gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                       <  1:14215/149‑1 (MQ=255)        gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACtt                                                                                                                                       >  6:121638/1‑149 (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:252081/149‑8 (MQ=38)                     aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTT‑CAGTAC‑GAACGTTCCTAcaactccg                                                                                                                          <  6:246555/149‑8 (MQ=38)                            ccAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                   <  2:332863/149‑1 (MQ=255)                            ccAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                   >  2:60207/1‑149 (MQ=255)                                 tgctgcGCATCGAAGAGGAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGGTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGt                                                                                                               <  1:60207/149‑1 (MQ=255)                                             aaGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTgg                                                                                                >  1:23778/1‑149 (MQ=255)                                                agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCt                                                                                             >  5:379566/1‑149 (MQ=255)                                                    aCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑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‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAgg                                                                                     >  1:334457/1‑149 (MQ=255)                                                         gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTAc                                                                                                                                              <  8:347412/93‑1 (MQ=255)                                                         gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTAc                                                                                                                                              >  7:347412/1‑93 (MQ=255)                                                          cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCGAGGCCAGGCt                                                                                   <  5:121638/149‑1 (MQ=255)                                                           gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCACGGCTAAGCAAGTGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAgg                                                                                     <  7:180687/146‑1 (MQ=255)                                                           gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCACGGCTAAGCAAGTGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAgg                                                                                     >  8:180687/1‑146 (MQ=255)                                                             cAAAGCGGTTTACCGCTGTCGTGCCTAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCGAGGCCAGGCTTCg                                                                                <  5:273526/149‑1 (MQ=255)                                                              aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGt                                                                               <  6:379566/149‑1 (MQ=255)                                                                     tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                        >  5:150940/1‑149 (MQ=255)                                                                     tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                        >  2:43506/1‑149 (MQ=255)                                                                     tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                        >  2:37959/1‑149 (MQ=255)                                                                       tACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCg                                                                       <  5:343471/148‑1 (MQ=255)                                                                         ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑Ct                                                                                                  >  1:41648/1‑119 (MQ=255)                                                                         ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑Ct                                                                                                  <  2:41648/119‑1 (MQ=255)                                                                          cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGATATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGTCTTCGTg                                                                              <  7:185919/138‑1 (MQ=255)                                                                          cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCCGTCGACGAAG‑CTCTGGCCTCGGCCAGGCTTCGTGCTATTCGAgg                                                                    <  1:37959/149‑1 (MQ=255)                                                                          cgcgGTCGTGCCGAGGTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGGTTCTGTCGACGAAGCC‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAgg                                                                    <  1:43506/149‑1 (MQ=255)                                                                                 gTGCCGAGTTTCGCGGCTAAGCAAGAGATCGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGAT‑A‑GTACTCAACTTTCCTACGTACAAACGCGTTTCTGTCGACGAAG‑C‑CTCGCCACGGCCAGGCATCGTGCTAATCGAGGCCCCGaac                                                            >  6:364314/1‑148 (MQ=255)                                                                                  tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                           <  8:182581/149‑1 (MQ=255)                                                                                  tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                           <  2:334457/149‑1 (MQ=255)                                                                                    ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTa                                                          >  7:116276/1‑148 (MQ=255)                                                                                               ggCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                              <  8:116276/149‑1 (MQ=255)                                                                                                 cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                              >  3:98909/1‑147 (MQ=255)                                                                                                 cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                              <  4:98909/147‑1 (MQ=255)                                                                                                     gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑t                                       >  5:305889/1‑149 (MQ=255)                                                                                                            aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAAc                                >  6:353495/1‑149 (MQ=255)                                                                                                            aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTGCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAAc                                <  6:211271/149‑1 (MQ=255)                                                                                                               gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTTCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCt                              <  2:312350/148‑1 (MQ=255)                                                                                                               gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGc                                                                                         >  3:270199/1‑90 (MQ=255)                                                                                                               gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGc                                                                                         <  4:270199/90‑1 (MQ=255)                                                                                                                   aaaGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACGAACGGGTTTCTGTCGACGACG‑C‑CTGGGCGCGGCCAGGGTGCGTGGTATTCGAGGCCCCGAAGTAACAAGTACGCGGGCGCC‑GTTTTAACCTGGAt                          >  5:195073/1‑148 (MQ=255)                                                                                                                    aaGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATaa                        <  6:150940/149‑1 (MQ=255)                                                                                                                     aGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAAc                       <  2:23778/149‑1 (MQ=255)                                                                                                                         agcagcCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGa                   >  1:202343/1‑149 (MQ=255)                                                                                                                              ccGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCt                              <  7:318790/133‑1 (MQ=255)                                                                                                                              ccGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTAGTGCTATTAGAGGCCCCGAAATAACAAGTAAGCGGCAGCC‑TTTTTAAACt                              >  8:318790/1‑133 (MQ=255)                                                                                                                                ggAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTg                             <  7:137692/132‑1 (MQ=255)                                                                                                                                ggAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CAGGCCACGGCCAGGCATAGTGGGATTAGAGGCCCCGAAATAAAAAATAACCGGCAACCATTTTTAAcag                              >  8:137692/1‑130 (MQ=255)                                                                                                                                   tacgTAGGAAAGTT‑CAGTAC‑GAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTcc         >  7:81272/4‑149 (MQ=255)                                                                                                                                       tAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAATTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt     >  6:282951/1‑149 (MQ=255)                                                                                                                                       tAGGAAAGTT‑CAGTAC‑GAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt     >  5:110353/1‑149 (MQ=255)                                                                                                                                          gAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg  >  7:209220/1‑149 (MQ=255)                                                                                                                                          gAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg  <  6:305889/149‑1 (MQ=255)                                                                                                                                          gAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGACGCGCAGTCCCTATTgg  >  2:244930/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                       CCGTTCCGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTT‑C‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAG‑C‑CTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCC‑TTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGG  >  NC_002947/1810140‑1810425

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