Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	388,722	T→A	11.3%	L198Q (CTG→CAG)	soxB →	sarcosine oxidase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	388,722	0	T	A	11.3%	175.5 / 10.1	69	L198Q (CTG→CAG)	soxB	sarcosine oxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base A (5/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

CACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTG  >  NC_002947/388576‑388870                                                                                                                                                    |                                                                                                                                                     cacTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTgg                                                                                                                                                    <  8:251744/149‑1 (MQ=255)  acTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCATGGGGCTTCGCCCGCGCTGCCGACGCCCTggg                                                                                                                                                   >  1:284993/1‑149 (MQ=255)  acTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCATCGGAGCCACCGTACAGCGCCGTGGTGGCGTTGCACGTCACGATGCCGT‑GGCCTGCGGCTTCGCCCGCGCTGCCGACGCCCTcgg                                                                                                                                                   >  2:132739/1‑149 (MQ=255)    tGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGTCGACGCCCAGGGCg                                                                                                                                                 >  5:58979/1‑149 (MQ=255)      cccAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGTCGACGCCCAGGGCGt                                                                                                                                                <  1:62237/148‑1 (MQ=255)      cccAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTc                                                                                                                                               >  4:29088/1‑149 (MQ=255)          ggTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGTCGACGCCCAGGGCGTCggca                                                                                                                                           <  6:58979/149‑4 (MQ=255)           gTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCt                                                                                                                                          >  3:135322/1‑149 (MQ=255)                ggCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑AGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATc                                                                                                                                      >  5:164847/1‑148 (MQ=255)                    gAAATCCCGTACCTGGACTGCTCGAAGACCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCCCGAGGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGc                                                                                                                                  >  7:63391/1‑148 (MQ=255)                    gAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGc                                                                                                                                  >  1:1781/1‑148 (MQ=255)                    gAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTCGCGTTGCCCGTCACGATGCCCT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGc                                                                                                                                  >  7:170634/1‑148 (MQ=255)                          ccGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaaa                                                                                                                               <  8:14187/145‑1 (MQ=255)                          ccGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaaa                                                                                                                               >  7:14187/1‑145 (MQ=255)                               ccTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGa                                                                                                                      <  3:311354/149‑1 (MQ=255)                                cTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGAt                                                                                                                     >  5:250767/1‑149 (MQ=255)                                   gACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATc                                                                                                                    <  2:80053/147‑1 (MQ=255)                                   gACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCg                                                                                                                   <  1:304043/148‑1 (MQ=255)                                   gACTGATCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATc                                                                                                                    <  2:74555/147‑1 (MQ=255)                                         tCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCg                                                                                                            >  1:222273/1‑149 (MQ=255)                                         tCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGACACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGTCATCGGCTTCCg                                                                                                            >  7:126116/1‑149 (MQ=255)                                            aagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGt                                                                                                                                                >  8:210127/1‑110 (MQ=255)                                            aagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGt                                                                                                                                                <  7:210127/110‑1 (MQ=255)                                            aagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCaa                                                                                                         <  3:116670/149‑1 (MQ=255)                                             agaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAg                                                                                                        <  5:356389/149‑1 (MQ=255)                                                  acCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGaaaa                                                                                                   >  3:54413/1‑149 (MQ=255)                                                    ccGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACg                                                                                                 >  8:128072/1‑149 (MQ=255)                                                     cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACgg                                                                                                >  6:82939/1‑149 (MQ=255)                                                      gCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgc                                                                                               <  8:103301/149‑1 (MQ=255)                                                       cTACCCGCTCCGCGGCGCCACCGTACAGCGCCGTGGTGGCGGTGCCCGTCACGCTGCCGG‑TGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTTATCCAGCAAACCGAACTGATCGGCTTCCGCAAGGAAAACGgcg                                                                                              <  6:164847/149‑1 (MQ=255)                                                           ccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGt                                                                                          >  3:44170/1‑149 (MQ=255)                                                            cGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAg                                                                                                        >  8:105672/1‑134 (MQ=255)                                                            cGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAg                                                                                                        <  7:105672/134‑1 (MQ=255)                                                                    ggCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                 <  3:29088/149‑1 (MQ=255)                                                                         cACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGG‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaaa                                                                            <  6:20637/149‑1 (MQ=255)                                                                           tacggtggcgcCGAGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaaa                                                                            <  4:361381/140‑1 (MQ=255)                                                                           ccGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAAc                                                                           >  2:178383/1‑148 (MQ=255)                                                                             gTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGAACCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACcaa                                                                        >  8:236515/1‑149 (MQ=255)                                                                                cAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACcaacaa                                                                     <  4:135322/149‑1 (MQ=255)                                                                                 aGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACaaa                                                                    <  7:128072/149‑1 (MQ=255)                                                                                     ccGTGGTGGCGTTGCCCGTCACGATGCCGTAGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGc                                                                 >  1:35399/1‑149 (MQ=255)                                                                                           tggCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCg                                                          <  2:35399/149‑1 (MQ=255)                                                                                            ggCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCgg                                                         >  6:159426/1‑149 (MQ=255)                                                                                             gCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGgc                                                        <  3:165924/149‑1 (MQ=255)                                                                                                   cccGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCaa                                                   <  2:1781/148‑1 (MQ=255)                                                                                                    ccGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAgc                                                 <  5:82939/149‑1 (MQ=255)                                                                                                    ccGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAgc                                                 <  7:236515/149‑1 (MQ=255)                                                                                                      gTCACGATGCCTT‑GGCCTGTTGCTTCGACCTCGCTGCCGACTCCCTTGGCGTCTACCTGCTCAATCAAACCTAAGTGATCGGCTTCCTCAAGGAAAAATTCGAGGTCATCGTTGTGTAAACCAACAAAGGCTTCATCGGCTCCAAgcgc                                               <  8:63391/149‑1 (MQ=255)                                                                                                             tGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCggtg                                        <  2:284993/149‑3 (MQ=255)                                                                                                               ccGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGt                                       <  1:178383/148‑1 (MQ=255)                                                                                                                 cggggCCTGGTGCCTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTgg                                     <  1:132739/146‑1 (MQ=255)                                                                                                                       cTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccg                               >  3:67765/1‑149 (MQ=255)                                                                                                                         ggggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgcc                             >  8:152810/1‑149 (MQ=255)                                                                                                                            gCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                          <  4:54413/149‑1 (MQ=255)                                                                                                                            gCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCATCCGCAAGGAAAACGGCGCGGTCATCGGGGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                          <  3:51618/149‑1 (MQ=255)                                                                                                                             cTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                          >  8:173780/1‑148 (MQ=255)                                                                                                                                 gCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAACCGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTc                     >  1:101409/1‑149 (MQ=255)                                                                                                                                 gCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTc                     <  4:44170/149‑1 (MQ=255)                                                                                                                                   ccGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCg                   >  7:91938/1‑149 (MQ=255)                                                                                                                                      tcaggtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg                  <  8:314971/141‑1 (MQ=255)                                                                                                                                      tcaggtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg                  >  7:314971/7‑147 (MQ=255)                                                                                                                                         tGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGcaca             >  3:318992/1‑149 (MQ=255)                                                                                                                                         ggtcGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGcaca             >  1:200882/4‑149 (MQ=255)                                                                                                                                          gtcGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACAt            >  1:162461/3‑149 (MQ=255)                                                                                                                                                gCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGcc        >  1:362670/1‑147 (MQ=255)                                                                                                                                                 cccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAg     >  7:242330/1‑149 (MQ=255)                                                                                                                                                 cccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAg     >  4:251359/1‑149 (MQ=255)                                                                                                                                                  ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                          <  6:340464/127‑1 (MQ=255)                                                                                                                                                  ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                          >  5:340464/1‑127 (MQ=255)                                                                                                                                                    tGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTg  >  7:228217/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGT‑GGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTG  >  NC_002947/388576‑388870

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