Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,291,467	G→T	10.1%	L86L (CTG→CTT)	estC →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,291,467	0	G	T	10.1%	205.3 / 11.8	69	L86L (CTG→CTT)	estC	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base T (3/4);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACT  >  NC_002947/1291323‑1291609                                                                                                                                                  |                                                                                                                                               cGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGccc                                                                                                                                            <  4:66200/149‑1 (MQ=255) cGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGccc                                                                                                                                            >  3:66200/1‑149 (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)                                                      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                                                                        >  3:323009/1‑148 (MQ=255)                                                                      cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGt                                                                        >  2:319646/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                                                                            >  7:332602/1‑125 (MQ=255)                                                                                         tGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCtt                                                                            <  8:332602/125‑1 (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               <  5:66308/142‑1 (MQ=255)                                                                                                                                    atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa               >  6:66308/3‑144 (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)                                                                                                                                               gCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTAt      <  4:102544/142‑1 (MQ=255)                                                                                                                                               gCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTAt      >  3:102544/1‑142 (MQ=255)                                                                                                                                                  |                                                                                                                                               CGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTACCG‑CCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACT  >  NC_002947/1291323‑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: ‑