Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,796,713	C→T	9.4%	S427S (TCG→TCA)	vdh ←	vanillin dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,796,713	0	C	T	9.4%	179.0 / 6.9	64	S427S (TCG→TCA)	vdh	vanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/25);  new base T (3/3);  total (36/28)
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.74e-01

GGCCGTGCTGGAGGGTGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAGGAC  >  NC_002947/3796565‑3796860                                                                                                                                                     |                                                                                                                                                    ggCCGTGCTGGAGGGTGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATAccc                                                                                                                                                     <  3:85034/149‑1 (MQ=255) ggCCGTGCTGGAGGGTGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATAccc                                                                                                                                                     <  3:202661/149‑1 (MQ=255)           gAGGGTGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCAc                                                                                                                                             <  7:55855/147‑1 (MQ=255)                tGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTgg                                                                                                                                      >  4:91682/1‑149 (MQ=255)                tGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTgg                                                                                                                                      <  7:283593/149‑1 (MQ=255)                    ccAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGcc                                                                                                                                   >  5:294482/1‑148 (MQ=255)                                 gTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTTATATGGCAGATACCCGACTCTACCCGTTGGGCCAAGGCCAGGGcgcg                                                                                                                     <  1:196313/149‑1 (MQ=255)                                     aCTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTg                                                                                                                 <  3:91682/149‑1 (MQ=255)                                     aCTGATCGAGGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTg                                                                                                                 <  2:94004/149‑1 (MQ=255)                                       tGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGgt                                                                                                               >  6:298364/1‑149 (MQ=255)                                                gCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCt                                                                                                      >  6:45476/1‑149 (MQ=255)                                                            ccGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCg                                                                                           <  6:294482/148‑1 (MQ=255)                                                            ccGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGa                                                                                          >  2:242076/1‑149 (MQ=255)                                                             cGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTTATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGAt                                                                                         <  5:36134/149‑1 (MQ=255)                                                              gAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTcgc                                                                                                          >  2:367079/1‑131 (MQ=255)                                                              gAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTcgc                                                                                                          <  1:367079/131‑1 (MQ=255)                                                               aaGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGc                                                                                                                              >  2:65922/1‑110 (MQ=255)                                                               aaGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGc                                                                                                                              <  1:65922/110‑1 (MQ=255)                                                                aGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGAc                                                                                      <  5:45476/149‑1 (MQ=255)                                                                     ccATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGa                                                                                   <  7:160661/147‑1 (MQ=255)                                                                           ccGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTc                                                                           >  1:76557/1‑149 (MQ=255)                                                                             gCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTcc                                                                          >  5:17835/1‑148 (MQ=255)                                                                             gCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTcc                                                                          >  5:73547/1‑148 (MQ=255)                                                                                 gACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTc                                                                     >  7:155088/1‑149 (MQ=255)                                                                                  aCTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAg                                                                                                                                <  8:343738/89‑1 (MQ=255)                                                                                  aCTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAg                                                                                                                                >  7:343738/1‑89 (MQ=255)                                                                                       aCCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCTGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCgg                                                                                                        >  1:311716/1‑108 (MQ=255)                                                                                       aCCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCTGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCgg                                                                                                        <  2:311716/108‑1 (MQ=255)                                                                                       aCCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCTGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCgg                                                                                                        <  2:316047/108‑1 (MQ=255)                                                                                       aCCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCTGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCgg                                                                                                        >  1:316047/1‑108 (MQ=255)                                                                                         cccGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGc                                                               >  2:22147/1‑147 (MQ=255)                                                                                          ccGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAg                                                            >  5:323348/1‑149 (MQ=255)                                                                                            gCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGa                                                                                       >  2:251223/1‑120 (MQ=255)                                                                                            gCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGa                                                                                       <  1:251223/120‑1 (MQ=255)                                                                                             ccAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTg                                                         >  2:236999/1‑149 (MQ=255)                                                                                              cAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTg                                                         >  3:308161/1‑148 (MQ=255)                                                                                                 aCGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTgcag                                                      >  2:218363/1‑148 (MQ=255)                                                                                                    gCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAg                                                  <  5:298364/149‑1 (MQ=255)                                                                                                    gCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAg                                                  >  5:305018/1‑149 (MQ=255)                                                                                                        cTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGTGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGa                                                                        >  7:48615/1‑123 (MQ=255)                                                                                                        cTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGTGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGa                                                                        <  8:48615/123‑1 (MQ=255)                                                                                                         tGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTtc                                             >  3:242256/1‑149 (MQ=255)                                                                                                          gcgcTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTtcg                                            >  1:236836/1‑149 (MQ=255)                                                                                                          gcgcTTAATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTtcg                                            >  2:194569/1‑149 (MQ=255)                                                                                                            gcTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTtcgtc                                          <  2:76557/149‑1 (MQ=255)                                                                                                             ctgcATCGGGCAGGGTCGGGCCGTTGATATGGCAGATACCCGTCTCCACCCGTTGGGCGAAGGCCAGGGCGCGGGTGGTGTCGAGGCTGAATATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTTCTGCAAGGCTtcgtc                                          <  6:288117/145‑1 (MQ=255)                                                                                                              ttCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCTGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCt                                               <  5:84986/142‑1 (MQ=255)                                                                                                              ttCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCTGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCt                                               >  6:84986/1‑142 (MQ=255)                                                                                                                    gTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGACCTCCGAGTCGTTGGCCAGCTGCAGCAAg                                                  >  8:43910/1‑133 (MQ=255)                                                                                                                    gTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAg                                                  <  7:43910/133‑1 (MQ=255)                                                                                                                     tGAACGGTCGGGCTGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGgc                                 <  1:194569/149‑1 (MQ=255)                                                                                                                       aaCGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGgcgc                               >  6:295228/1‑149 (MQ=255)                                                                                                                               ggCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTacc                       <  2:236836/149‑1 (MQ=255)                                                                                                                                  cGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTaccact                    <  4:242256/149‑2 (MQ=255)                                                                                                                                  cGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTaccacc                    >  6:368061/1‑149 (MQ=255)                                                                                                                                  cGTTGATATGGCAGATACCCGACTACACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTaccacc                    >  3:332667/1‑149 (MQ=255)                                                                                                                                      gATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTAccaccgcca                >  3:292818/1‑149 (MQ=255)                                                                                                                                      gATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTAccaccgcc                 <  8:155088/148‑1 (MQ=255)                                                                                                                                        tatGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTAccaccgccacc              >  3:285191/1‑149 (MQ=255)                                                                                                                                           ggCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCggg           >  8:341449/1‑149 (MQ=255)                                                                                                                                            gCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGc          >  6:211853/1‑149 (MQ=255)                                                                                                                                            gCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGc          >  8:116755/1‑149 (MQ=255)                                                                                                                                               gATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGa       <  1:242076/149‑1 (MQ=255)                                                                                                                                                aTACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGaa      <  4:308161/149‑1 (MQ=255)                                                                                                                                                 tACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAg     <  1:22147/149‑1 (MQ=255)                                                                                                                                                  aCCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAg     >  7:188289/1‑148 (MQ=255)                                                                                                                                                    ccGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAGGa   >  4:314507/1‑148 (MQ=255)                                                                                                                                                    ccGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAGGAc  >  3:59863/1‑149 (MQ=255)                                                                                                                                                    ccGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAGGAc  <  5:368061/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGCCGTGCTGGAGGGTGACCCAGCGCAACTGGGTGAACTGATCGATGGCCGTGCGGCTGCCGAAGCTGCCATAGCCGCTGGACTTGACCCCGCCAAACGGCATCTGCGCTTCATCGTGAACGGTCGGGCCGTTGATATGGCAGATACCCGACTCCACCCGTTGGGCCAAGGCCAGGGCGCGGCTGGTGTCGCGGCTGAAAATGGCCGATGACAGACCGAACTCCGAGTCGTTGGCCAGCTGCAGCAAGGCTTCGTCGCCTTCGGCGCGCAGTACCACCGCCACCGGGCCGAAGGAC  >  NC_002947/3796565‑3796860

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