Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,318,911	G→C	7.0%	A42A (GCC→GCG)	treSA ←	trehalose synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,318,911	0	G	C	7.0%	174.8 / 4.2	57	A42A (GCC→GCG)	treSA	trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (23/30);  new base C (2/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

CAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGC  >  NC_002947/3318764‑3319052                                                                                                                                                    |                                                                                                                                               cAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTggg                                                                                                                                               >  4:74878/1‑149 (MQ=255) cAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTggg                                                                                                                                               >  2:81430/1‑149 (MQ=255) cAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTggg                                                                                                                                               <  5:351065/149‑1 (MQ=255)  aGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGc                                                                                                                                              >  2:203838/1‑149 (MQ=255)     cccTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTc                                                                                                                                           <  1:81430/149‑1 (MQ=255)       cTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCgg                                                                                                                                         <  4:163022/149‑1 (MQ=255)             gCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAgg                                                                                                                                   >  1:382594/1‑149 (MQ=255)              ccGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAggg                                                                                                                                  >  6:362542/1‑149 (MQ=255)                  tctcCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTtgc                                                                                                                              >  2:82994/1‑149 (MQ=255)                      cGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGt                                                                                                                                 <  5:93493/142‑1 (MQ=255)                        aCAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCcac                                                                                                                       <  3:74878/149‑1 (MQ=255)                              gcTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGgcgc                                                                                                                 >  1:354249/1‑149 (MQ=255)                                ttCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGa                                                                                                               <  5:238421/149‑1 (MQ=255)                                ttCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGa                                                                                                               >  6:366862/1‑149 (MQ=255)                                  ccACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGa                                                                                                               <  5:244063/147‑1 (MQ=255)                                     caACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGc                                                                                                          >  6:183307/1‑149 (MQ=255)                                      aaCGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGcc                                                                                                         >  8:288367/1‑149 (MQ=255)                                       aCGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCt                                                                                                        <  2:382594/149‑1 (MQ=255)                                        cGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTg                                                                                                       <  2:347302/149‑1 (MQ=255)                                         gCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGa                                                                                                      <  3:288253/149‑1 (MQ=255)                                         gCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGa                                                                                                      <  2:239222/149‑1 (MQ=255)                                         gCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGa                                                                                                      <  7:205179/149‑1 (MQ=255)                                          cTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAg                                                                                                     <  5:362542/149‑1 (MQ=255)                                            tCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTa                                                                                                   >  5:67498/1‑149 (MQ=255)                                                   ccAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTg                                                                                            >  3:272400/1‑149 (MQ=255)                                                       ggCACCGAGCACCGAGCAACCCTCGGGCGCGATGAGGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGccc                                                                                        <  7:373742/149‑1 (MQ=255)                                                          accgagcaccgagcAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGt                                                                                     >  8:66429/1‑149 (MQ=255)                                                           ccgagcaccgagcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGTTTAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTc                                                                                    <  1:101022/149‑1 (MQ=255)                                                           ccgagcaccgagcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCCGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTc                                                                                    >  1:409134/1‑149 (MQ=255)                                                              agcaccgagcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTg                                                                                 <  5:366862/149‑1 (MQ=255)                                                              agcaccgagcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTCCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTg                                                                                 <  1:82994/149‑1 (MQ=255)                                                               gcaccgagcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTgg                                                                                <  6:317597/149‑1 (MQ=255)                                                                    gagcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGCCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGCAGGc                                                                           >  5:355171/1‑149 (MQ=255)                                                                     agcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGcc                                                                          <  4:272400/149‑1 (MQ=255)                                                                      gcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCt                                                                                  <  3:36916/140‑1 (MQ=255)                                                                      gcaACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCt                                                                                  >  4:36916/1‑140 (MQ=255)                                                                           ccTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACAGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAg                                                                     >  1:76609/1‑148 (MQ=255)                                                                               gggCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTg                                                                                            <  8:279687/121‑1 (MQ=255)                                                                               gggCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTg                                                                                            >  7:279687/1‑121 (MQ=255)                                                                               gggCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATg                                                                 <  5:183307/148‑1 (MQ=255)                                                                                     gatgatGGCGTCGGGGTAGACCGTGAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGCGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGt                                                          >  8:115021/1‑149 (MQ=255)                                                                                     gatgatGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTg                                                                                                       <  3:21085/104‑1 (MQ=255)                                                                                     gatgatGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTg                                                                                                       >  4:21085/1‑104 (MQ=255)                                                                                     gatgatGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGt                                                          >  3:235887/1‑149 (MQ=255)                                                                                     gatgatGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGCGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGt                                                          >  8:120692/1‑149 (MQ=255)                                                                                                     tAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGAc                                          <  2:354249/149‑1 (MQ=255)                                                                                                         ccGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTAt                                      <  1:203838/149‑1 (MQ=255)                                                                                                           gTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCTCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGa                                    <  4:235887/149‑1 (MQ=255)                                                                                                             cAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCTAGCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGAc                                   >  1:115398/1‑149 (MQ=255)                                                                                                              aGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGGT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACg                                  <  1:253784/148‑1 (MQ=255)                                                                                                              aGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACg                                  <  7:66429/148‑1 (MQ=255)                                                                                                                ccACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGt                               >  8:308768/1‑149 (MQ=255)                                                                                                                 cacaCCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTc                              >  5:245479/1‑149 (MQ=255)                                                                                                                   cacCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCgg                            <  6:245479/149‑1 (MQ=255)                                                                                                                     ccGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGc                          >  5:304124/1‑149 (MQ=255)                                                                                                                        aGGCGATTTCGGTGGCGCGGCGGGGCTGCGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTgg                       <  7:115021/149‑1 (MQ=255)                                                                                                                        aGGCGATTTCGGTGGCGCGGCGGGGCTGCGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTgg                       <  7:120692/149‑1 (MQ=255)                                                                                                                                 cGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTcg              >  3:99687/1‑149 (MQ=255)                                                                                                                                            cggGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTg   <  2:115398/149‑1 (MQ=255)                                                                                                                                             gtggCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGTCTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGGGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGc  <  6:355171/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               CAGGCCCTGTACGCCGATCTCCGACAGGCGCTTCCACAACGCTTCGTGGGCCAGGGCACCGAGCACCGAGCAACCCTCGGGCGCGATGATGGCGTCGGGGTAGACCGTCAGCCACACCGAGGCGATTTCGGTGGCGCGGCGGGGCTGGGCCTCGGCATAGGGTTGCT‑GCCACAGGCGCGACTGGCCTGAGTACAGGCTGGCCCGGTCCTGGGAGGCCTTGAGCATGGCGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGC  >  NC_002947/3318764‑3319052

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