Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,872,965	C→T	6.6%	G219S (GGC→AGC)	aqpZ ←	aquaporin Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,872,965	0	C	T	6.6%	227.0 / 4.4	76	G219S (GGC→AGC)	aqpZ	aquaporin Z
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base T (2/3);  total (36/40)
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.93e-01

GTCCAACAGCATCGAAGAAACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGG  >  NC_002947/4872828‑4873107                                                                                                                                           |                                                                                                                                                gTCCAACAGCATCGAAGAAACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCAc                                                                                                                                      >  6:119678/1‑149 (MQ=255) gTCCAACAGCATCGAAGAAACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCAc                                                                                                                                      <  4:126236/149‑1 (MQ=255)          cATCGAAGAAACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAgcg                                                                                                                             >  2:276081/1‑149 (MQ=255)            tCGAAGAAACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAgcggc                                                                                                                           >  6:263485/1‑149 (MQ=255)                agaaACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCTGATCGGTGCTGCTATCGGCGgcgcgc                                                                                                                       <  4:153310/149‑3 (MQ=255)                 gaaACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCAc                                                                                                                      >  1:294752/1‑149 (MQ=255)                  aaaCCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCAcc                                                                                                                     <  1:18463/149‑1 (MQ=255)                      cTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAg                                                                                                                  >  7:160914/1‑148 (MQ=255)                      cTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAg                                                                                                                  >  1:57036/1‑148 (MQ=255)                      cTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGa                                                                                                                 <  6:67730/149‑1 (MQ=255)                       ttGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCTGATCGGTGCTGCTATCGGCGgcgc                                                                                                                         <  3:285033/140‑1 (MQ=255)                       ttGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCTGATCGGTGCTGCTATCGGCGgcgc                                                                                                                         >  4:285033/1‑140 (MQ=255)                       ttGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGaa                                                                                                                >  3:287893/1‑149 (MQ=255)                            cGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGc                                                                                                            <  5:263485/148‑1 (MQ=255)                             gTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCca                                                                                                          <  7:50260/149‑1 (MQ=255)                                 tCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAg                                                                                                                               >  2:293906/1‑124 (MQ=255)                                   accaccCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAAc                                                                                                               >  4:79023/1‑138 (MQ=255)                                   accaccCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAAc                                                                                                               <  3:79023/138‑1 (MQ=255)                                   accaccCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAgctg                                                                                                    >  4:206718/1‑149 (MQ=255)                                    ccaccCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAgctgc                                                                                                   <  4:287893/149‑1 (MQ=255)                                       cccTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGctgctg                                                                                                 >  1:282185/1‑148 (MQ=255)                                        ccTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGctgctg                                                                                                 >  4:258584/1‑147 (MQ=255)                                         cTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAg                                                                                              <  2:282185/149‑1 (MQ=255)                                         cTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAg                                                                                              <  5:119678/149‑1 (MQ=255)                                               gATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGc                                                                                                            >  3:65139/1‑129 (MQ=255)                                               gATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGc                                                                                                            <  4:65139/129‑1 (MQ=255)                                               gATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCa                                                                                        <  8:159455/149‑1 (MQ=255)                                                     cgtaccgtacGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCAcc                                                                                  >  1:76770/1‑149 (MQ=255)                                                         ccgtacGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAg                                                                                                       >  8:223254/1‑124 (MQ=255)                                                         ccgtacGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAg                                                                                                       <  7:223254/124‑1 (MQ=255)                                                         ccgtacGTGAAACCTACGCCGCCGGGCGCCGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCTGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACcgac                                                                               >  5:123367/1‑148 (MQ=255)                                                            tacGTGAAACCTATGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCAcc                                                                                                                     >  5:180103/1‑107 (MQ=255)                                                            tacGTGAAACCTATGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCAcc                                                                                                                     <  6:180103/107‑1 (MQ=255)                                                            tacGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAAc                                                                           <  8:160914/149‑1 (MQ=255)                                                               gTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAgg                                                                        >  1:50811/1‑149 (MQ=255)                                                                    aCCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGGTAAGGTTCTTTCGCCAGGCCCCGATACCGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGggccg‑g                                                                  <  8:78596/149‑1 (MQ=255)                                                                        aCGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑Gcccc                                                              >  1:148945/1‑149 (MQ=255)                                                                          gccgccGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACcgacg                                                                              <  1:44264/132‑1 (MQ=255)                                                                          gccgccGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACcgacg                                                                              >  2:44264/1‑132 (MQ=255)                                                                          gccgccGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACcga                                                                                <  2:52/130‑1 (MQ=255)                                                                          gccgccGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACcga                                                                                >  1:52/1‑130 (MQ=255)                                                                            cgccgGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGc                                                          >  4:283912/1‑149 (MQ=255)                                                                               cgGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCGGtagaac                                                       >  3:126059/1‑144 (MQ=255)                                                                                 ggCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTAcgc                                                     >  4:32983/1‑149 (MQ=255)                                                                                 ggCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTAcgc                                                     <  2:57036/149‑1 (MQ=255)                                                                                  gcgcTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGGCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAgg                                                                                             >  1:88534/1‑109 (MQ=255)                                                                                  gcgcTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAgg                                                                                             <  2:88534/109‑1 (MQ=255)                                                                                  aacgctggTGCACGCAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTAcgc                                                     >  8:44631/8‑149 (MQ=255)                                                                                     cTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCgg                                                 >  5:54730/1‑149 (MQ=255)                                                                                     cTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCgg                                                 <  1:276081/149‑1 (MQ=255)                                                                                            cAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGATCAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGAc                                          <  4:37299/149‑1 (MQ=255)                                                                                            cAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGa                                           <  6:54730/148‑1 (MQ=255)                                                                                            cAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGa                                           >  7:215428/1‑148 (MQ=255)                                                                                            cAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGa                                           <  2:76770/148‑1 (MQ=255)                                                                                             acgcAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCATCGGCGCCACCCACAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGAGCCCCGTGCTACGCGCGGGGTTGa                                           <  7:44631/146‑1 (MQ=255)                                                                                                cAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACaga                                       <  1:299927/149‑1 (MQ=255)                                                                                                 aGCGTTAAGGTTCTTTCGACAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACGGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACagag                                      <  8:215428/149‑1 (MQ=255)                                                                                                        aGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCTGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTg                                <  6:123367/148‑1 (MQ=255)                                                                                                          gTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGATGACAGAGGTATTGGTc                             >  6:100162/1‑149 (MQ=255)                                                                                                           ttCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCCTGCAACCCGCCGGGTTTCAAAAGGTATTTGta                             >  2:311982/1‑147 (MQ=255)                                                                                                                  gCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑Gcccc                                                              >  4:158994/1‑108 (MQ=255)                                                                                                                  gCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑Gcccc                                                              <  3:158994/108‑1 (MQ=255)                                                                                                                  gCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTg                                <  6:302662/138‑1 (MQ=255)                                                                                                                  gCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTg                                >  5:302662/1‑138 (MQ=255)                                                                                                                   ccAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGa                      >  6:16927/1‑147 (MQ=255)                                                                                                                   ccAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGa                      <  5:16927/147‑1 (MQ=255)                                                                                                                   ccAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATc                    <  4:124139/149‑1 (MQ=255)                                                                                                                                cAGCGCGCCGCCGATAGCAGCACCGTTCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGa       <  1:311982/149‑1 (MQ=255)                                                                                                                                cAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTgg        <  3:258584/148‑1 (MQ=255)                                                                                                                                 aGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAg            >  5:255018/1‑143 (MQ=255)                                                                                                                                 aGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAg            <  4:15263/143‑1 (MQ=255)                                                                                                                                 aGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAg            <  6:255018/143‑1 (MQ=255)                                                                                                                                 aGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAg            >  3:15263/1‑143 (MQ=255)                                                                                                                                   cgcgcCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtca    <  2:294752/149‑1 (MQ=255)                                                                                                                                   cgcgcCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtca    <  2:148945/149‑1 (MQ=255)                                                                                                                                     cgccgcCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcagg  <  3:283912/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                GTCCAACAGCATCGAAGAAACCTTGTTCGTGCTCACCACCCTGGGCGATGACCGTACCGTACGTGAAACCTACGCCGCCGGGCGCTGCGTGCAC‑CAGCGTTAAGGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCG‑GCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGG  >  NC_002947/4872828‑4873107

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