Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,004,478 G→A 13.6% L74L (CTG→TTG)  yigI ← 4HBT thioesterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,4780GA13.6% 186.2 / 2.8 81L74L (CTG→TTG) yigI4HBT thioesterase family protein
Reads supporting (aligned to +/- strand):  ref base G (21/45);  new base A (8/2);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 5.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCATTGGTTAATGAAATGTTGATTTTATGCTTCTTTGTTG  >  NC_000913/4004216‑4004751
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                           
cTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGATACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCAGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAGGGTTATTCCCACGATAAACAGACAGGCGGCGagacc                                                                                                                                                                                                                                                                             >  1:77523/1‑274 (MQ=255)
             ttCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCCCACCACACCCGCGGCGACATCCCCCGCCGa                                                                                                                                                                                                                                                                  >  2:373087/1‑274 (MQ=255)
                                ggcggtggcACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCGGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTAATTCCCACGCACAACAGACCGGCGGCGACATCCAGCGCCCACGCAATGACCCCCCCGTGccaca                                                                                                                                                                                                                                          >  1:45454/1‑274 (MQ=255)
                                ggcggtggcACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAAAGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTACTTCATCTTCACTGATGGTTTCGTGGCGGTTTAATGTATTTCCCACCCACACCACACCTCCGGCGACATCCCGCGCCCACGCCATGACCCCGCCGTCCaaac                                                                                                                                                                                                                                          >  1:139747/1‑277 (MQ=255)
                                ggcggtggcACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCACGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGCTAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACACTCAACACCGCCGCAATGACCCCTCCGCGCaaaa                                                                                                                                                                                                                                          >  1:89150/1‑278 (MQ=255)
                                     tggcACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTACAGGACTTCCCACGCCCAACAGAACGGGCGCGACACCCAGCTCCCCACCAAAAACCCCCCCGTGCaaaaa                                                                                                                                                                                                                                         >  4:330321/1‑273 (MQ=255)
                                             gCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTgcg                                                                                                                                                                                                                               <  2:330321/274‑1 (MQ=255)
                                                  ataAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAg                                                                                                                                                                                                                          <  2:318427/274‑1 (MQ=255)
                                                   tacaGCTGGTCGTCATTGTGTAATGCAACGCGGGAGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGcc                                                                                                                                                                                                                     <  1:129994/275‑1 (MQ=255)
                                                   taAAGCTTTTCTTCATTGTGTAATTCAACGCGGGCGCCGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGcc                                                                                                                                                                                                                     <  1:46494/278‑1 (MQ=255)
                                                      agttgtgCTGCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAGTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCGTGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCGGCTGGCGTAGTTCATGTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGcc                                                                                                                                                                                                                     <  4:164800/272‑1 (MQ=255)
                                                       gCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCCAACCGGCGGCGACATCCAGCGCCGACGCCATGACCCCGCCGTGCACAATGCTTTTCGCCCCGTTGccc                                                                                                                                                                                                                    >  1:139287/1‑275 (MQ=255)
                                                        cTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCcca                                                                                                                                                                                                                   <  2:160631/275‑1 (MQ=255)
                                                           ttcttcATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTTAGGTACTTCCCCCCCACACCAGACCGGCCGCGGCACCCAGCCCCGAAGCAATGACCCCGCCGGTTAAAATGCTTTTGCGCCAGTTGCccacca                                                                                                                                                                                                               >  2:86814/1‑275 (MQ=255)
                                                           ttcttcATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTTCCACGCACACCAGACTGGCGGCGGCATCCAGCGCCGACGCAATGACCCCCCCGCTCAAAATGCCTGTCGCCCAGTTGCccacca                                                                                                                                                                                                               >  2:92659/1‑275 (MQ=255)
                                                           ttcttcATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGGTCGTTCCACGCACACCACACCGGGGGCGAGCTCCACCGCCCACGCCATGACCCCCCCCGCTCAAAACTCTTGCCCCCAGTTTCccaccc                                                                                                                                                                                                                >  4:186594/1‑274 (MQ=255)
                                                                   ttgtacaTTCAAGCGCGGGGAACGGCGGCTTTGTTGCCTGCAGGCCACATGCTACTAGTGGCAGTCAAACGATCGCCCCTTCCTGGGCACAGATAAGTAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGTTACTTCCCGCGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTgg                                                                                                                                                                                                        <  4:53956/269‑1 (MQ=255)
                                                                       tAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGGACTTTCCACGCACACCAGACCGGCGGCGACATACAGCGCCGACGCAATGACCCCGCCCTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCCTTGGCTGa                                                                                                                                                                                                    >  2:496835/1‑274 (MQ=255)
                                                                            cAACGCGGGGGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTa                                                                                                                                                                                              <  2:186594/275‑1 (MQ=255)
                                                                                cgcgGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACTCGAAGAACAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTCAGGGTCTTTCCCCCCACCCCAAAACAGCGTGCGCATCCAGGCGCGCAGCCATTGCCCCCCCCGGCAAAAATCCTTGCCCCCCGTTGCACACCCTCATTCACTGAATTTTaaaa                                                                                                                                                                                          >  4:360557/1‑274 (MQ=255)
                                                                                cgcgGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTATGTCCCACGCACACCCGACCGGCGGCGACATCAGCCGCCCACCCAATTACCCCCCCCTGCAAAATGTGTTTCGCCCACTTGCCCACCATCATTCGCTAGTTTTTcaac                                                                                                                                                                                          >  4:37913/1‑274 (MQ=255)
                                                                                cgcgGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTAGTTCCCACGCACACAAGACCGGCGGCGGCCTCCAGCGCCGACGCAATGACCCCGCCGTGTCCAATGCTTTGCGCCCAGTTGCCCACCAACATTggagttgttttcaag                                                                                                                                                                                           >  2:137353/1‑262 (MQ=255)
                                                                                cgcgGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTTCCACGCCCACCAGACCGGCCGCGAAATCCCGCGCCCACGCCACTACCCCCCCCTGCACAAAGCCTTGCGCCCACTTTCCCCCCCTCATTCTCTCTTTTTTAAAg                                                                                                                                                                                          >  4:216365/1‑275 (MQ=255)
                                                                                cgcgGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGGGCCGACGCCATGGCCCCGCCCTGCAAAATGCCTTGCGCCCAGTTGCCCACCCTCATTGGCTGATTTTTAAAg                                                                                                                                                                                          >  2:368527/1‑275 (MQ=255)
                                                                                cgcgGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAAACCGGCGGCGAAAACCCGCGCCGACGCAATGAACCCGCCGTGCAAAATGCTTTTCGCCCAGTTGCCCACCAACATTGGCTGATTTTTAAAg                                                                                                                                                                                          >  2:545848/1‑275 (MQ=255)
                                                                                   ggggngcgacgacTTGTTTGCCGCGACACAACAGGCCCCGAGAAGCNAAAAAACTCACGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCACCATCCGCNGTAGCCGCTGGCGTAGGTCATCGTCACTGGTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGAACACCAGACCGGCGGCGACATCCAGCGCAGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCcca                                                                                                                                                                                                                   <  3:93368/240‑1 (MQ=255)
                                                                                    nggaacggcgacTTTAGTGCCGCCACGCAACAAGCTACCGGAGGCGGTAAAACCNTGGCNCCTCCCTGGGCGCGGGTAAGCAANGCGACGATCAATGGTCCCCGTCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGc                                                                                                                                                                                                                                                               <  4:594684/203‑1 (MQ=255)
                                                                                    ggcgacggcgacTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAt                                                                                                                                                                                                   <  2:216365/263‑1 (MQ=255)
                                                                                    ggcgacggcgacTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAAGTACTTCCCACGCACACCAGAACGGCGGCCGCATCCACCCCCGCCCCCATGACCCCCCCGTGCACAATTCTTTGCGCCCAAGCTCCCCCCATCTTTTGCCGATTTTTAAAAGcc                                                                                                                                                                                       >  4:105392/1‑274 (MQ=255)
                                                                                      cgacggcgacTTTGTTGCCTGAACGCAACAGGCAACTAGTAGCAGTAAACCGCTCTCCCTGGCCTGGGCGCAGAGAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGc                                                                                                                                                                                    <  4:99587/275‑1 (MQ=255)
                                                                                            tgacTTAATTGCCTGTACGGAACAGGCTCCTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATGAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGCTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCg                                                                                                                                                                              <  4:496835/274‑1 (MQ=255)
                                                                                                         tGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCAGGTTAAGGTACTTCCCACGCACACCACACCGGCCGCGCCCTCCAGCGCCGACGCAATGACCCCGCCGTTCCAAATGCTTTGCGCCCAGTTGCCCCACATCCTTGGCTGATTTTTAAATGCCAGCCGTTCGAAATCTTTTTct                                                                                                                                                                >  4:259096/1‑273 (MQ=255)
                                                                                                          gCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGCAGATCAACGGTCCCCATCCGCGATAGCCGCTGGCGAAATTCATCTTCACTGATGGTTTCGTGCCGGGTCATGGTACATCCTACGCACCCCAAGCCCGCGGCCACATCCAGCGCCGACCCACTCAACCCGCAGTGCCACATGCTTTTCGCCCACTCTCCCACCATCAATTGCCCATTTTTAACAGGCCAGTGGTGACATCTCtttgctg                                                                                                                                                                 >  4:81931/1‑270 (MQ=255)
                                                                                                            aaggaaCAGGCTACTATTGGCAGTAAACCGCTCGCCCCGGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCCTCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTaa                                                                                                                                                             <  4:368527/270‑1 (MQ=255)
                                                                                                                  cAGGCTACTGGTAGCAGTAAAACGCTCGCCCCGGCCTGGGCGCCGATAATCAACGCGAAGATAAATGGTCCCCACCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCa                                                                                                                                                      <  4:545848/275‑1 (MQ=255)
                                                                                                                     gCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTGCTTCCCCCGCACACCCAACCGGGGGCGACCTCCCAGGCCGACGCAACGGACCCCCCGCTCCAAAAGCTTTTCGCCCAGTTTCCCCACCCCATTGGCCGAGTTTTAAAGGCCAGCTCTGCGCAACTTTTTTCCTTACCCTCCAGt                                                                                                                                                    >  4:197797/1‑274 (MQ=255)
                                                                                                                     gCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGaa                                                                                                                                                                           <  1:78896/252‑1 (MQ=255)
                                                                                                                     gCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAAGATCTTCCCACCCACCCCAGACCGGCGGGGAAACCCAGTGCCCGCCGCATTACCCCCCCTCGCAACAAGCTTTGCCTCCCGTTTGCCCCCCACCTTGGCAGTGTTTTAAAGGTCAGCTGTTACAAGTCTTTTCTCGTactccccagt                                                                                                                                                     >  4:235407/1‑265 (MQ=255)
                                                                                                                      ctactaGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTc                                                                                                                                                  <  2:171613/275‑1 (MQ=255)
                                                                                                                       tactaGTAGCAGTAAAACGCGCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTcc                                                                                                                                                 <  2:600643/275‑1 (MQ=255)
                                                                                                                            gtagCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACCCAATGACCCCGCCGTGCAAAATGGTTTGCGCCCAGTTGCCCACCAACCTTGGCTGATTTTTAAAGGCCAGCTGTGGGAACTCTTTTTCGGAACGCTCCCGGTCCATtcc                                                                                                                                            >  2:339075/1‑275 (MQ=255)
                                                                                                                               gCAGTAAAACGCTCGCCCTTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCCTCACTGATGGTTTCGTGGCGGGTTTAGGTAATTCCCATCCACACCCAAACGGACGCGACATCACGAGCCGACACAATTACCCCGCCGTTCAAAATGCTTTTCGCCCAGTTTCCCACACACACTGGCCGATTTTTTAAGGCGACGTGGGCGCAATCTTTTTTTTAAAGCTCCACtttcctgccaca                                                                                                                                         >  4:149293/1‑265 (MQ=255)
                                                                                                                                cggTAAAACCCTCGCCCCTGCCTGGGCGCCGATAATCAACCCGAAGCTCCATGGTCCCCATCCGCGATAGCCGCTGGCGTNGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATgcg                                                                                                                                     <  3:139287/276‑1 (MQ=255)
                                                                                                                                    aaaaCGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATgcg                                                                                                                                     <  2:584410/274‑1 (MQ=255)
                                                                                                                                       aCGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATgcg                                                                                                                                     <  2:572713/271‑1 (MQ=255)
                                                                                                                                             gCCCCTGCCTGGGCGCAGATAATCAAAGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCAATGATGGTTTCGTTGCCGGTTGACGCGCTTCTCCCGCACACCGATCCGGGGTCTTCAACCAGCGACGGAGCAATCAAACCGACATGCCACATCATTTGCGGCCCTTTGACCACCGACTATTGGTGATTTTTAAAGGCCAGTTGGCGCAATGTTTTTTCGTGACCACTCAGTTCCATTCCCAAttggcgcgtaaaag                                                                                                                           >  4:36144/1‑262 (MQ=255)
                                                                                                                                               cccTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGc                                                                                                                         <  2:197797/275‑1 (MQ=255)
                                                                                                                                                  tgcctgGGCGCGAATAAGTAACGCGAAGATCAATGTTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTTCTGGGTTCGTGGCGGGTTAAGGTACTTCCCACGCAAACCAGACCGGCGGCGACAGCCAGCGCCGACGCAATGACACCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGc                                                                                                                         <  4:198832/272‑1 (MQ=255)
                                                                                                                                                         gcgcAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCACACCCGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTTCGCCCAGTTGCCCCCCATCATTGGCTGATTTTTAAAGGCCAGCGGGGCGAACTCTTTTTCGTac                                                                                                                                                             >  1:647/1‑228 (MQ=255)
                                                                                                                                                                    cAACGCGAAGATCAATGGTCCCCATCCGCGGTAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAATTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCa                                                                                                    <  2:105392/275‑1 (MQ=255)
                                                                                                                                                                      taagagacagcAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCAc                                                                                                   <  2:565298/264‑1 (MQ=255)
                                                                                                                                                                       cgcgAAGAGCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACGGATGGTTTCGTGGCGGGGTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACt                                                                                              <  1:137787/278‑1 (MQ=255)
                                                                                                                                                                           aaGACCAATGGTCCCCAATCGCGGTAGCCGCTGGTGTATTTCATCTTCACTGATGGTTGCGTGGGGGGTTAAGGTACTACCCACGCGCACCAGGACGGCGGCGACAACCAGCGCCGACGCAATGACCCCGCCGTGCATAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTGAAGGCCAGCTGTGCGAACTCTTTTTCATAACGCTCCAGTTCCAACCCCGATGCGCGGTTAAATGGCATATGATAAACAAACATCTCTCCCACt                                                                                              <  4:92659/274‑1 (MQ=255)
                                                                                                                                                                                          cATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGtt                                                                              <  2:360557/275‑1 (MQ=255)
                                                                                                                                                                                            tCCGCGATCGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACGTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCa                                                                            <  2:286330/275‑1 (MQ=255)
                                                                                                                                                                                               gcgctgggctccggCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCTGCGGCGACATCCAGCGCCGACGCCATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCAACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCg                                                                         <  4:339075/263‑1 (MQ=255)
                                                                                                                                                                                                cgATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCAACCCAGACCGGCGGCGACATCCAGCGCCGGCGCAATGTCCCCGCCGTGCAACATTCTTTGCGCCCAATTGCCCACCATCCATTGCCGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCTTTACCCTCCAGCTTCATCCCCAATGCCCGGTTAAATGGCAAGTATTTAACAACAATCTCACCCCCCACTTTCCAGGCCTTTTCatca                                                                         >  4:388396/1‑271 (MQ=255)
                                                                                                                                                                                                 gttggCCGCTTGCGTATGTCTTCTTCCCTGATGGTTTCGTGGCGGGTTAAGGTACTTCCACCGCAAACCAGACCGTCGCCGACATCCATCGCCGACGCAATGACCCCGCCGTGCAAAATGCTGTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCCGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGt                                                                       <  4:306120/271‑1 (MQ=255)
                                                                                                                                                                                                    aGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTGCTTCCCACCCCCACCAGACCGGCGGCGCCATCCCACGCCCACCCAATTACCCCCCCCTGCCAAAAGCCTTTCGCCCCGTTTCCCACCAACACTGGGTCATTTTTAAAAGGCAACTGCTCGGACACTTTTTCGTAAAGCTCCCTTTCCATCCCCCATGCGCGGTTTAATGGGCTTTGATATAACAACACATCCCCCACTAATTTCAGGGCTTGTTCAACGGTCa                                                                     >  4:56170/1‑274 (MQ=255)
                                                                                                                                                                                                      ccGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGTTCTTCCCACCCACACCAAACCGGCGGCGACATTCAGCGCCGCCGCAATTAACCCCCCGTGCAAAATGCTTTGCGCCCACTTGCCCACCAACATTTGCCGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCACTTCCATCCCCAATTCGCGGTTCAATGGCACGGTTTACACAAAACTCTCCCCCACAAAATTCCGGGCTTTTTTCGCTGTCACTa                                                                  >  4:239466/1‑275 (MQ=255)
                                                                                                                                                                                                         cTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGGACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACgg                                                               <  2:37913/275‑1 (MQ=255)
                                                                                                                                                                                                           ggCGTAGTGCATCTTCACGGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCAAACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCa                                                             <  2:61288/275‑1 (MQ=255)
                                                                                                                                                                                                           ggCGTAGGTCAGCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCa                                                             <  2:115236/275‑1 (MQ=255)
                                                                                                                                                                                                           cgCGTATTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCa                                                             <  2:289924/274‑1 (MQ=255)
                                                                                                                                                                                                            gggTAGTTAATCTTCACGGATGGTTTAGGGGAGGGTTAAGGTATTTCCCCCGGACACCAAACCGGCGGCGACGTCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTGTTTCGTAACGCTCCAGTTCCATCCCCAAGGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCa                                                             <  4:298825/272‑1 (MQ=255)
                                                                                                                                                                                                              gTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGaa                                                                                                                                                                           <  2:306465/163‑1 (MQ=255)
                                                                                                                                                                                                               tAGTTCATTGTCACTGATTGTTTCGTGGCGGGGTAAGGCACTTCCCACGCACACCAGACCTGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACatat                                                      <  1:153512/278‑1 (MQ=255)
                                                                                                                                                                                                                  ttCATCTTCACTGATGGTTTCTTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCCGCGCCGACGCAATTACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACatat                                                      <  2:439878/275‑1 (MQ=255)
                                                                                                                                                                                                                      tCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCa                                                                                                                                                                                                                                                              <  2:547075/78‑1 (MQ=255)
                                                                                                                                                                                                                           acggATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTa                                              <  2:112746/271‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTTTCGTGGCGGGTTAAGGTCCTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGGGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCa                                       <  2:527705/275‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCa                                       <  2:25411/275‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTTTCGTGGCGGGTTAAGGTACTACCCACGCACACCAAACCGGCGGCTACATCCAGCACCGACGAAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCa                                       <  2:170465/275‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATcc                                                 >  4:546487/1‑242 (MQ=255)
                                                                                                                                                                                                                                                         tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCNAGTTCCATCCCCAATGCCCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTTTTCNGCGGTCAGTACGGCAGACATCTCATCCTTANNATTCCTTTGTTAATGAAATTtttattttattnt           >  3:73173/1‑276 (MQ=255)
                                                                                                                                                                                                                                                         tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCATTGGTTAATGAAATGTTGATTTTCTGct           >  1:66755/1‑279 (MQ=255)
                                                                                                                                                                                                                                                            cACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTac                                             <  2:546487/242‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATa                                                                                                                 >  2:444698/1‑172 (MQ=255)
                                                                                                                                                                                                                                                              cGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTaa                                                                                            >  4:345252/1‑193 (MQ=255)
                                                                                                                                                                                                                                                                  cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTc                                                                                        <  2:345252/193‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCATTGGTTAATGAAATGTTGATTTTATGCTTCTttgttg  <  2:235407/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                           
CTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCAGACATATCATCCTTACACTTCATTGGTTAATGAAATGTTGATTTTATGCTTCTTTGTTG  >  NC_000913/4004216‑4004751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: