Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,291,462	A→C	9.6%	K85Q (AAG→CAG)	estC →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,291,462	0	A	C	9.6%	208.4 / 10.6	71	K85Q (AAG→CAG)	estC	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base C (3/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGGTGGCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACT  >  NC_002947/1291315‑1291609                                                                                                                                                     |                                                                                                                                                    gtgggtggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCaa                                                                                                                                                    <  5:301929/149‑1 (MQ=255)    ggtggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGc                                                                                                                                                  <  4:178195/148‑1 (MQ=255)    ggtggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCt                                                                                                                                                 >  6:122420/1‑149 (MQ=255)         cGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGccc                                                                                                                                            >  3:66200/1‑149 (MQ=255)         cGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGccc                                                                                                                                            <  4:66200/149‑1 (MQ=255)         cGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGccc                                                                                                                                            <  8:249745/149‑1 (MQ=255)            aaCCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTc                                                                                                                                          >  6:241914/1‑148 (MQ=255)             aCCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCGTCCCTGCAGTTGGTCGGCGAGGGCCAGCTTGCCCTCg                                                                                                                                         <  1:246995/149‑1 (MQ=255)                      gACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGt                                                                                                                               >  2:27308/1‑149 (MQ=255)                           gtgGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGcc                                                                                                                           >  2:273374/1‑148 (MQ=255)                            tgGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGAGCAAGCTGGCCCTCGATGCCCCGGTGGcc                                                                                                                           <  1:27308/147‑1 (MQ=255)                                 cGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAactact                                                                                                                    >  5:169312/1‑149 (MQ=255)                                          cGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGt                                                                                                           <  1:109404/149‑1 (MQ=255)                                          cGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGt                                                                                                           >  2:78991/1‑149 (MQ=255)                                            aCAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTg                                                                                                                              >  3:268338/1‑128 (MQ=255)                                            aCAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTg                                                                                                                              <  4:268338/128‑1 (MQ=255)                                            aCAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTg                                                                                                                   >  8:279964/1‑139 (MQ=255)                                            aCAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTg                                                                                                                   <  7:279964/139‑1 (MQ=255)                                            aCAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTc                                                                                                         >  4:236681/1‑149 (MQ=255)                                             cAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCg                                                                                                        <  5:122420/149‑1 (MQ=255)                                              aaGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGc                                                                                                       <  3:236681/149‑1 (MQ=255)                                               aGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGcc                                                                                                      <  8:183356/149‑1 (MQ=255)                                                ggACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGccc                                                                                                     <  1:78991/149‑1 (MQ=255)                                                ggACGGCCAGCAGGCCGGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCCAGACCTTCGCGGCCGTTACCC‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGccc                                                                                                     <  4:199927/149‑1 (MQ=255)                                                      ccAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCgg                                                                                               <  8:84611/149‑1 (MQ=255)                                                        agcagGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGGCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGt                                                                                             <  7:176655/149‑1 (MQ=255)                                                          cagGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTaa                                                                                           <  1:219125/149‑1 (MQ=255)                                                           agGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAg                                                                                          <  1:273374/149‑1 (MQ=255)                                                             gCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCt                                                                                                                                                 >  6:321621/1‑92 (MQ=255)                                                             gCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCt                                                                                                                                                 <  5:321621/92‑1 (MQ=255)                                                              ccTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCa                                                                                        <  7:242838/148‑1 (MQ=255)                                                               cTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCTGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGt                                                                                      <  6:230212/149‑1 (MQ=255)                                                               cTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGt                                                                                                                               <  5:369206/108‑1 (MQ=255)                                                               cTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGt                                                                                                                               >  6:369206/1‑108 (MQ=255)                                                               cTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGt                                                                                      >  6:157294/1‑149 (MQ=255)                                                               cTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGt                                                                                      >  1:229040/1‑149 (MQ=255)                                                                    acaGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATc                                                                                 <  5:241914/149‑1 (MQ=255)                                                                              cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGt                                                                        >  2:319646/1‑148 (MQ=255)                                                                              cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGt                                                                        >  3:323009/1‑148 (MQ=255)                                                                                tCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAg                                                                     >  4:261294/1‑149 (MQ=255)                                                                                  gCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGCTGGTCGGCGAGGGCAATCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAgct                                                                   >  2:138021/1‑149 (MQ=255)                                                                                     aaCGTGTTCTCCTGCACCATCACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGTGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCACAGGCGGGTAAGCAGTCGATCACCTTGc                                                                          <  7:221068/139‑1 (MQ=255)                                                                                     aaCCAGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGACGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTAGCCCAGGCGGGTAAGCAGTCGATCACATTGc                                                                          >  8:221068/1‑139 (MQ=255)                                                                                          gTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAg                                                                                          <  3:168353/118‑1 (MQ=255)                                                                                          gTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAg                                                                                          >  4:168353/1‑118 (MQ=255)                                                                                          gTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGc                                                            >  5:139356/1‑148 (MQ=255)                                                                                                 tGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCtt                                                                            <  8:332602/125‑1 (MQ=255)                                                                                                 tGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCtt                                                                            >  7:332602/1‑125 (MQ=255)                                                                                                 tGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCCCCGTGc                                                    >  6:314557/1‑149 (MQ=255)                                                                                                      cAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGcc                                               >  4:68356/1‑148 (MQ=255)                                                                                                      cAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCCATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGcc                                               >  4:73332/1‑148 (MQ=255)                                                                                                           ccTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCgg                                                 >  1:161409/1‑142 (MQ=255)                                                                                                           ccTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCgg                                                 <  2:161409/142‑1 (MQ=255)                                                                                                            cTTAGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCCGTCGCTCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCgg                                         >  6:243155/1‑149 (MQ=255)                                                                                                                   gCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCg                                  >  5:212425/1‑149 (MQ=255)                                                                                                                   gCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCg                                  >  2:238463/1‑149 (MQ=255)                                                                                                                           cG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCAGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAgctgc                          >  4:296423/1‑149 (MQ=255)                                                                                                                                   cAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCgg                    <  6:169312/148‑1 (MQ=255)                                                                                                                                   cAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGc                   <  1:238463/149‑1 (MQ=255)                                                                                                                                      ttGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGa                <  4:323009/149‑1 (MQ=255)                                                                                                                                        ggTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa               <  3:261294/148‑1 (MQ=255)                                                                                                                                         gTCGGCGAGGGCAAGCTGGCCCTCTATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCCTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGc             <  3:68356/149‑1 (MQ=255)                                                                                                                                         gTCGGCGAGGGCAAGCTGGCCCTCCATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCCTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGc             <  3:73332/149‑1 (MQ=255)                                                                                                                                            ggCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTTCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt          >  8:6151/1‑149 (MQ=255)                                                                                                                                            atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa               >  6:66308/3‑144 (MQ=255)                                                                                                                                            atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa               <  5:66308/142‑1 (MQ=255)                                                                                                                                             tcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt          <  6:258663/147‑1 (MQ=255)                                                                                                                                             tcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt          >  5:258663/2‑148 (MQ=255)                                                                                                                                                aGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTg                                                                           >  2:86054/1‑80 (MQ=255)                                                                                                                                                aGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTg                                                                           <  1:86054/80‑1 (MQ=255)                                                                                                                                                 gggCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTAt      >  5:199248/1‑148 (MQ=255)                                                                                                                                                    cAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACt  <  1:319646/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTGGGTGGCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACT  >  NC_002947/1291315‑1291609

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