Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,990,014 T→C 9.7% E251G (GAG→GGG)  hemC ← hydroxymethylbilane synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,990,0140TC 9.7% 240.8 / 2.1 99E251G (GAG→GGG) hemChydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base T (45/38);  new base C (7/5);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 101 total aligned reads displayed.

GCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCTAAACGTTTTAGTCCGGC  >  NC_000913/3989777‑3990252
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
gCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GTACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCACCCAGATTTCGCCTTCAATAAGATc                                                                                                                                                                                                                                                  >  2:875814/1‑239 (MQ=255)
 cTCACTAACTCTTCCCCAGGGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                                                                                                                                                                                                                                  <  2:646919/238‑1 (MQ=255)
  tCACTTACTCTTCACCAGCGGGAGACGGGCGGGTTAAAAGTATACCCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCCCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCg                                                                                                                                                                                                                                                 <  2:37446/238‑1 (MQ=255)
    aatAACTCTTCCCCAGCGGGAGACGGGCGGGTTACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCg                                                                                                                                                                                                                                              <  2:826649/237‑1 (MQ=255)
       aaCTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                                                                                                                                                                                                                            >  2:202606/1‑238 (MQ=255)
              ctcCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGCCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                                                                                                                                                                                                                                    >  2:561459/1‑239 (MQ=255)
               tcCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCCGCGTAGCTACCAAt                                                                                                                                                                                                                                   >  2:81323/1‑239 (MQ=255)
                  ggggggAGGCGGGCGGGTGAAAAGGATGCTCATGCCGGGTCGTCTCCGTTATAGCCTTCAGCGAGGATCT‑CCGGCGCTCCGTTATTCAGTAGCTCCTCTGCCAGCGAAACTCCCATGTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                                                                                                                                                                                                                                <  1:181764/236‑1 (MQ=255)
                    cGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGGGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                                                                                                                                                                                                                               >  2:345420/1‑238 (MQ=255)
                    cGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCCCGGCGTAGCTACCAATTGGc                                                                                                                                                                                                                               >  2:13876/1‑238 (MQ=255)
                    cGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑AGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCGGCGAAATCCCCCTTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGATTAATCTGCGAACCGTCCGGCGCGCCGCCCAGCGCACCCGGCCGGATTTCGCCCTCAATAAGATCGGCGTGGCCCCCAAATGGc                                                                                                                                                                                                                               >  1:813312/1‑238 (MQ=255)
                    cGGGAGACGGGCGAGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATACCCATTTGTTCGGAATCTTGCGGC‑GCACCGCGGCGTTCACCCCGAACAATCTGCCAACCTCCCGGCGCGCCGACCAGCGCACGCCGCCAGATTTCGCCATCAATAAGCCCGGAGTAGCtccacatgggc                                                                                                                                                                                                                               >  1:812959/1‑228 (MQ=255)
                        ggaCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGCTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCt                                                                                                                                                                                                                           <  2:623095/237‑1 (MQ=255)
                                      aaaGGATTCTCAGGCCGGGGCGTCTCCTGTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCga                                                                                                                                                                                                            <  2:677966/238‑1 (MQ=255)
                                        aGGATACTCATGCCGGGGCGTTTCCGTTATAGGCTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                                                                                                                                                                                                          <  2:850470/239‑1 (MQ=255)
                                          gtttcccctGCCGGGGCGTTCCCGTTATAGCCTTCAGCGAGGATCC‑CGCGGGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag                                                                                                                                                                                                           <  1:248681/228‑1 (MQ=255)
                                           attcTCATGCCGGGGCGTCTCCGTTATAGACTTCAGGGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGGCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCTCCTGACATCCGCCTTCGAGACg                                                                                                                                                                                                        <  1:34189/235‑1 (MQ=255)
                                            tACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                                                                                                                                                                                                                                   >  2:718436/1‑194 (MQ=255)
                                                 aTGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTc                                                                                                                                                                                                 >  2:151926/1‑239 (MQ=255)
                                                    ccGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCCCGAATAATCTGCGAACCGTCCGGCGCGCCGCCCCACCCCCGCAGCCGCATTTCGCCATCAAAACGCCCGCCGTAGCACCCCATTGGCACCCCACATCCCCCCTCCAAACGGGGATTAATg                                                                                                                                                                                              >  1:682795/1‑239 (MQ=255)
                                                        ggCGTCCCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgcg                                                                                                                                                                                           <  2:649247/238‑1 (MQ=255)
                                                          ggtttcCTTTATAGCCTCCAG‑GGGGGTCTCCGCGCCCTCCGTTATTCAGTAGCCCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGGCCGGCGCGGCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGcgcg                                                                                                                                                                                         <  1:202606/234‑1 (MQ=255)
                                                          cGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGcgcg                                                                                                                                                                                         <  2:682795/238‑1 (MQ=255)
                                                                        aCTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGCCGTAGCTACCAATTGGCACCTGCCATCCGCCTTCGAGCCGGGTATTCATGCCGCGTCCTGCGGTAAccc                                                                                                                                                                           >  2:319665/1‑236 (MQ=255)
                                                                         cTTCAGCGAGGATCG‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAAcgcg                                                                                                                                                                          <  2:181551/238‑1 (MQ=255)
                                                                         acctcgcgggggaccccgcgcgcgCCTGTTTTCAGTAGCTCCTCCGCCAGGCAAAGCCCCATTTGTTCGGCATCTTGGGGG‑GCCCCGGGGGGTTAACCGCGAATAATCTGCGAACCGTCGGGCGCGCCGAACAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAAcgcg                                                                                                                                                                          <  1:875814/223‑1 (MQ=255)
                                                                          ttCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGCACCGTCCGGCGCGCCGACCAGCGCACGCACCCAGATTTCGCCATCACTAAGCCCGCCTCAACTCCCAATTGGCCCCTGACAACCGCCTTCGAGACGGGGACTCATGGCCCGttacccggggacgcccc                                                                                                                                                                        >  1:471027/1‑224 (MQ=255)
                                                                          ttCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACCTCCGCCTTCGAGACGGGTATTCATGGCCCGTTCTGCgggacaccgca                                                                                                                                                                        >  1:817718/1‑230 (MQ=255)
                                                                          ttCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                                                                                                                                                                                                          <  1:262759/205‑1 (MQ=255)
                                                                          ttCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                                                                                                                                                                                                          >  2:262759/1‑205 (MQ=255)
                                                                          ttCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACATGACATCCGCCTTCGAGACGGGGATTCATGGCGCGTTCTGCGGTAACgcgcg                                                                                                                                                                        >  2:450028/1‑238 (MQ=255)
                                                                          ttCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGCCACCTGACATCCGCCTTCGAGACGGGACTTCATGGCGCGTTCTGCGGTAACgcgcg                                                                                                                                                                        >  2:308750/1‑238 (MQ=255)
                                                                           tCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGGATTCATGGCGCGTTCTGCGGTAACGCgcag                                                                                                                                                                       >  2:780870/1‑239 (MQ=255)
                                                                            cAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTTCGGTAACGCgcagt                                                                                                                                                                      >  2:628966/1‑239 (MQ=255)
                                                                            cAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCgcagt                                                                                                                                                                      >  2:141383/1‑239 (MQ=255)
                                                                            cAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCgcagg                                                                                                                                                                      >  2:612994/1‑238 (MQ=255)
                                                                            cAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGACGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTGCTGCGGTAACGCgcag                                                                                                                                                                       >  2:190393/1‑238 (MQ=255)
                                                                            cAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGCCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCCCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCCTGGCGCGTTCTGCGGTACCGCgccgt                                                                                                                                                                      >  2:867113/1‑239 (MQ=255)
                                                                               cGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAACTGGCCCCTGACATCCGCCTTCGCGCCGGGTATTCATGGCGCGTTCTGCGGTAAAGCgcacggcg                                                                                                                                                                   >  1:413973/1‑234 (MQ=255)
                                                                                          gcgcgcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGAATTCATGGCGCGTTCTCCGGTAACGCGCATTCC‑AGTTTCgttgt                                                                                                                                                        >  2:175323/1‑239 (MQ=255)
                                                                                              gcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCCATAAGCACGGCGTAGCTACCACTTGCCACCTGCCATCGGCCTTCGGCAAGGGTATTCATGGCGCGTTCCCCGGAAAcgccccctgcaattccgaggagaat                                                                                                                                                     >  1:36034/1‑217 (MQ=255)
                                                                                               cgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGCTTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGGATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTc                                                                                                                                                   >  2:764105/1‑239 (MQ=255)
                                                                                               cgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCAACGCGCGTTCTCCGGTAACGCGCAGTGC‑AGTTTCGTGGGGAtt                                                                                                                                                    >  2:247617/1‑238 (MQ=255)
                                                                                                            tAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGCGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAagca                                                                                                                                      >  2:149324/1‑239 (MQ=255)
                                                                                                            tAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAagca                                                                                                                                      >  2:57810/1‑239 (MQ=255)
                                                                                                              gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGCCATCCGCCTTCGAGACGGGGATTCATTGCGCGTTCTGCGGTAACCCGCCGTGC‑CGTTTCGGGGTGATTCAGGCCGGCAagcag                                                                                                                                     >  1:192301/1‑238 (MQ=255)
                                                                                                               ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGAGCGGCAagcagc                                                                                                                                    >  2:47356/1‑238 (MQ=255)
                                                                                                               ctctTCTGCCAGCCAACTCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCCGCGTTCACCGCGAATAATCTGCGACCCGTCCGGCGCGCCGACCAGCGCACGCAGCAAGATTTGGCCATCACAAAGTCCGGCGTACCTACCAACTGCAACCTGCCAGCCGCCTTCGAGTCGGGTATTCATGTCGCGTCCTGCGGTACCGCTCACTGC‑AGTTTCTCGGCGATTCAGCGCGGCAagcagc                                                                                                                                    >  1:604567/1‑238 (MQ=255)
                                                                                                                 cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGGGGGAACGCCAAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGACGCTc                                                                                                                                  >  1:164787/1‑238 (MQ=255)
                                                                                                                 cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCCCGGCAACAAGCTc                                                                                                                                  >  2:274623/1‑238 (MQ=255)
                                                                                                                 cttctGCCAGCGAAATCACCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCCCCGTCCAGCCCGGGATTCAGGAGCCATTCTGCGGGAAC‑CGCAGTGCGAGTTTCCGGGGGTTTAAGCGCGGCaaaaacgaga                                                                                                                                 >  1:482410/1‑231 (MQ=255)
                                                                                                                         aGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCCGGGC‑AGTTTCGTGGGGTTTCAGCGCGGCAAGCCGCTCGCGGGTggg                                                                                                                         >  1:738182/1‑237 (MQ=255)
                                                                                                                            gcaacccctagttGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCCCCGGGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAACTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGa                                                                                                                      <  2:462889/228‑1 (MQ=255)
                                                                                                                            gAAATCCCCTTTTGTTCGGCATCTTGCGGG‑GCACCGCGGCGTTCCCCGGGAATAATCTGCGAACCGTCCGGCGCGGCGACCAGCGCACGCAGCCAGATTGCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGa                                                                                                                      <  1:248434/239‑1 (MQ=255)
                                                                                                                            gAAATCCCCATGTGTTCGGCATCTTGCGGC‑CCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTg                                                                                                                       <  2:15015/238‑1 (MQ=255)
                                                                                                                             aaaTCCCCATGTGTTCGGCATCTCGCGGC‑GCACCGCGGCGGTCACCGCGAATACTCTGCGAACCGTCCGGCTCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGa                                                                                                                      <  2:469201/238‑1 (MQ=255)
                                                                                                                                        ttttCGGCATTTTGCGGCTGCCCCG‑GGCGTTCCCCGCGTAGAAGCTCCGAACCTGCCGGCGCGCGTGCCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCGCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAggcg                                                                                                          <  1:55602/237‑1 (MQ=255)
                                                                                                                                        tttgcGGCCTTTTGCGGC‑GCCCCGCGGCGTTCCCCGCGAATAAGCTGCGAACCGGCCGGCGCGCCGACCAGCGCACGGAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAggcg                                                                                                          <  1:151926/235‑1 (MQ=255)
                                                                                                                                        tgttcgggtcgttGCGGC‑GCCCCGCGGCGTTCACCGCGAATAAACTGCGAACCGGCCGGGGCGCCGGCCAGCGCACGCAGCAAGATTTCGCCATCAATAAGCTCGGCGTAGCTAGCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAggcg                                                                                                          <  1:81323/228‑1 (MQ=255)
                                                                                                                                        tggtCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCCAACCGGCCGGCGCGGCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAggcg                                                                                                          <  2:164787/236‑1 (MQ=255)
                                                                                                                                            cGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTGGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGCATTCATGGCGCGTTCTCCGGTAACGCGCCGGGC‑GGTTTCGTGGTGATTAAGCGCGGGCAGCAGCCGGCGGGTGCGGGAATCCTCACgggggact                                                                                                      >  1:662677/1‑236 (MQ=255)
                                                                                                                                            cGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCa                                                                                                                                          <  2:577966/203‑1 (MQ=255)
                                                                                                                                            cGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCa                                                                                                                                          >  1:577966/1‑203 (MQ=255)
                                                                                                                                                 tCTTGCGGC‑CCCCCGCGGCGCTCCCCGGCAATAATCCGCCAACCGTCCGGCGCGCCCACCAGCCCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAAt                                                                                                 <  1:458890/239‑1 (MQ=255)
                                                                                                                                                  cTTGCGGC‑GCCCCGGGGCGTCCCCCGCGAATAAACTGGGAACCGTCCGGGCCGCCGACCAGCGCCCGCAGCCCGCTTTCGCCATCAATAAGCGCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATa                                                                                                <  1:686633/239‑1 (MQ=255)
                                                                                                                                                    tGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATAc                                                                                               >  2:700984/1‑238 (MQ=255)
                                                                                                                                                       cgc‑gcCCGGCGGCGTTCACCCCGAATAATCCGCGACCCGGCCGGCGCGCCGACCAGCGCCACCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCAc                                                                                           <  1:175323/238‑1 (MQ=255)
                                                                                                                                                         c‑gcCCCGCGGGTTTCCCCGCGAATAATCCGCGAACCGTCCGGCGCGGCGACCAGCGCTCTCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGTTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCAcc                                                                                          <  1:97722/238‑1 (MQ=255)
                                                                                                                                                           gcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCg                                                                                         <  2:354509/238‑1 (MQ=255)
                                                                                                                                                             aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGCGGGCGGGAATCATCAAGGCGGCATTCAATACCCACCGCa                                                                                       >  2:800041/1‑238 (MQ=255)
                                                                                                                                                                 cggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCCCCCCtt                                                                                   >  2:170385/1‑238 (MQ=255)
                                                                                                                                                                      ttccccggaaatattCTGCGACCCGCCCGGCGCGCCGCCCAGCGCACGCGGCCAGATTTCGCCATACATAAGCGCGGCGTAGTTCCCGTTGGGCCGCGGCCATCCGCCTGCGGGCGGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAACCACCAAGGCGGCATTCAATCCCCACCGCACCTGGCCCTa                                                                             <  1:13876/230‑1 (MQ=255)
                                                                                                                                                                       tCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGGGTGCGTGAATCATCAAGGCGGCATTCAATACCCCCCGCACCTTGTCCTa                                                                             >  2:19850/1‑238 (MQ=255)
                                                                                                                                                                         aCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTg                                                                                                                                                                                    >  1:155590/1‑134 (MQ=255)
                                                                                                                                                                         aCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTg                                                                                                                                                                                    <  2:155590/134‑1 (MQ=255)
                                                                                                                                                                            gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAag                                                                                                                                        >  2:568906/1‑174 (MQ=255)
                                                                                                                                                                            gccaataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAag                                                                                                                                        <  1:568906/171‑1 (MQ=255)
                                                                                                                                                                                 aaaccTGGGACCCGCCCGGCGCGCCGACCCGCGCACGCAGCCAGATTTCGCCATCAATAAGCCCGGCGGAGCTACCAATTGGCCCCTGACCTCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑Aaga                                                                 <  1:57810/235‑1 (MQ=255)
                                                                                                                                                                                     ctcccacccGCCCGGCGCGCCCACCAGCGCCCGCCGCCCAATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCCCCTGACACCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGTTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGaaa                                                             <  1:319665/232‑1 (MQ=255)
                                                                                                                                                                                               tCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCACTTGGCACCTGACATCCGCCTTCGAGACGGGCATTCATGGCCCGTTCTGCGGTAACCCGCACCGC‑AGTTTCGTGGGGCTTCAGCGCGGCAAGCAGCTCCGGAGTGCGCGTATCATCCAGGCGGCATTCCAAACCCCCCCCACCTTGTCCCACCGCCGG‑ACAAGaaaaccggggggg                                                   >  1:550340/1‑229 (MQ=255)
                                                                                                                                                                                               tCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGGGAATCATCAAGGCGGCATTCAATACCCACCGCCCCATGTCCTACCGCCGG‑CAGAGAAATCCCggcgg                                                    >  2:861881/1‑235 (MQ=255)
                                                                                                                                                                                                  gggccgccctcccgcgcACGAACACGGATTTCGCAACCAAAAAGCTCGGCTGAGCTACCAATGGGACGCTGACAGCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGGGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCAGCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCaa                                                <  1:867113/226‑1 (MQ=255)
                                                                                                                                                                                                     gcgcCGACCAGCGCACGCAGCCAGATTGCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAAcg                                              <  2:738182/238‑1 (MQ=255)
                                                                                                                                                                                                        ccGACCAGGGCACGCAGCCAGATTTCCCCATCAATAAGCTCGGCTTAGCTACCAATTGGCCCCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACgcggc                                          <  2:337038/239‑1 (MQ=255)
                                                                                                                                                                                                         cGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGACTCATCAAGGCGGCATTCACTACCCACCCCCCCTTCTCCCACCCCCGGAAACAGAAACCGCGGG‑GGCAAccccgcc                                         >  1:470736/1‑233 (MQ=255)
                                                                                                                                                                                                                  cgccggccccgttttCGCCATCAATAGGCTCTGCGTACCGACCACTTGGCCGCTCGCCTCACCCTTCGGGCGGTGTATGCATGGCCCTTTCTGCGGCAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGAAAGCAGCGCGCGAGGGCGTGAATCATCAAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑Gt                               <  1:183333/227‑1 (MQ=255)
                                                                                                                                                                                                                  cccgcgcccAGATTTCCCCTTCAATAAGCTCGCGGTAGCTCCCAATGGGCACCGGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCATCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑Gt                               <  1:99194/232‑1 (MQ=255)
                                                                                                                                                                                                                  cccgcaTCCAGCTGCCGCCATCAATAAGCTCGGCGTGGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCCGCGGGAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑Gt                               <  2:849005/237‑1 (MQ=255)
                                                                                                                                                                                                                     gcaGCCAGATTTCGCCATCAATAAGCTCGAAGAGGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGAATTCATGGCGCGTGCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATGAAGCGCGGCAAGCAGCTCGCAAGTGCGTGAGTCACAAAGGGGGCGTTCAATAACCACCGCACCTTCTCCTACCGCCGG‑ACGAGAGACCTCGGCGGGCAAcgcgcccgcccacaggga                              >  1:505153/1‑224 (MQ=255)
                                                                                                                                                                                                                        gctaGATTTCGCCATAAATAAGCTCGGGGTAGCTCCCAATTGGTACCTGACACCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTcc                         <  1:700984/236‑1 (MQ=255)
                                                                                                                                                                                                                        gCCAGATTTCGCCAACAATAAGCCCGGCGTAGCTACCAATTGGCACCTGACCACCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTcc                         <  2:602771/239‑1 (MQ=255)
                                                                                                                                                                                                                        ccctattttCCCCTTCAATAAGTTGTGCGTTGTGCCCAATTGGCACCTTGCATCCGCCTTCGAGGCGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTcc                         <  1:428398/233‑1 (MQ=255)
                                                                                                                                                                                                                              tttCCCCATCAAAAAGGCCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCt                   <  2:817718/239‑1 (MQ=255)
                                                                                                                                                                                                                                tccacATCAAAAAGCCTGGGGGTGCTACCAATTGGCCCCTGACACCCGCCTTCGGAGCGGGTATTCATGGCGGGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCTaa                 <  1:170385/235‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGCCTCCAGACCTaaa                >  2:259816/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAAAACCACACGAACCTTGACCTACCGCCGG‑AAAGGAAATCTCGGGTGGCAACCCAGCGC‑AACACAGTGACTCCAAACCTaaa                >  1:28108/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cGCCATCAATAAGCCCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGAGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGGGAAATCTCGGGTGGCAACGCGGCGCGAATCC‑GTGACTCCAgcctcaa                 >  2:687834/1‑232 (MQ=255)
                                                                                                                                                                                                                                    cATCAATAAGCCCGGCGTAGCTACCAATTGGCACCGGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCTAAACg              <  2:9193/238‑1 (MQ=255)
                                                                                                                                                                                                                                              cTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGGGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCTAAACGTTTTAGTcag    >  2:615276/1‑236 (MQ=255)
                                                                                                                                                                                                                                               tCGGCGTCGCGCCCAATTGGCACCCGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCTAAACGTTTTAGTCCGGc  <  2:482410/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
GCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCT‑CGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGC‑GCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGC‑AGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGG‑AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC‑GTGACTCCAGACCTAAACGTTTTAGTCCGGC  >  NC_000913/3989777‑3990252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑