Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,887 T→C 6.5% intergenic (+34/‑115) garD → / → prlF D‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,276,8870TC 6.5% 312.3 / 8.0 92intergenic (+34/‑115)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (38/48);  new base C (4/2);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

TGGCGGTACCCGTCATTAAAATGGCAACCCGCACCGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAGCCAGGCCAGGACAGCATTCATTACGAAATTCTGCCTGGTGGGC  >  NC_000913/3276621‑3277152
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                         
tGGCGGTACCCGTCATTAAAATGGCAACCCGCACCGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                                                                                                                                                   <  2:222218/275‑1 (MQ=255)
 ggCGGTACCCGGCATTAAAATGGCAACCCGCACCGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                                                                                                                                                   <  2:107585/274‑1 (MQ=255)
         ccGTCCTTAAAAGGGCAACCCCAACCGAGCTGGCGAACCGCGGGTTTGATGTAATGGATATTAATGCGGGCACCATCTATACCGGCGAAGAACCTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAGACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGACCTTAGTGCCATGTTTTTa                                                                                                                                                                                                                                                           <  4:193992/274‑1 (MQ=255)
            tCATTAAAATGGCAACCCGCACCGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGt                                                                                                                                                                                                                                                       >  1:67564/1‑275 (MQ=255)
                   aaTGGCAACCCGCACGGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAAGGCGGGAACCATCCCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAg                                                                                                                                                                                                                                                 <  1:145861/274‑1 (MQ=255)
                      ggCAACCCGCACCGAGCTCGCGAACCGCTGTTTTTATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGcccc                                                                                                                                                                                                                                             <  4:451237/275‑1 (MQ=255)
                           cccGCACCGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAATGCTGGCACCATCGCTACCGGCTAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACg                                                                                                                                                                                                                                         <  2:30293/274‑1 (MQ=255)
                                aCCGGGCTGGCGAACCGCTTGTTTGATTTAATGGATATTAATGCGGGCCCCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                                                                                                                                                   <  4:492208/275‑1 (MQ=255)
                                         gCGAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGGGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgca                                                                                                                                                                                                                          <  2:371910/275‑1 (MQ=255)
                                           gAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAACTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                                                                                                                                        <  2:189881/275‑1 (MQ=255)
                                               cgcgGGTTGGATTGAATGGATATTAATGCGGGAACCATCGTCACCTGCGAAGTAACTATTGAAGGTGTGGGCTGTAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGTGAAGAAGAAAACCTGTTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgca                                                                                                                                                                                                                    <  4:440416/275‑1 (MQ=255)
                                                cttgtttttttttAATGGGTATTAATGGGNNCCCCATCGCTGCCGGCGAAGAAACTATTGGNGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGNGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgcat                                                                                                                                                                                                                   <  3:67564/265‑1 (MQ=255)
                                                            tAAGGGATATTAATGCGGGCCACAGCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                                                                                                                                       <  4:62339/275‑1 (MQ=255)
                                                                        aTGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                                                                                                                                        >  2:299950/1‑246 (MQ=255)
                                                                           cGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTTTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAtt                                                                                                                                                                                            >  2:536278/1‑271 (MQ=255)
                                                                           cGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGat                                                                                                                                                                                        >  1:143162/1‑274 (MQ=255)
                                                                            gggAACCATCGCTACCGGCGAAGAAGCTATTGTAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGTAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgca                                                                                                                                                                                                                    <  4:299949/246‑1 (MQ=255)
                                                                             ggCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGtt                                                                                                                                                                                                                                                               >  2:86761/1‑202 (MQ=255)
                                                                             ggCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                                                                                                                                                      >  2:475980/1‑275 (MQ=255)
                                                                             ggCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGCCGACAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCCGCATAACCAGCTGGCGGTGTTTAACNCGGCACCGGCGACCTGCTTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGGGCCATGTTTTTATTTTTTAAAGCCCCCACGCCACTTAATAAttcattggttntn                                                                                                                                                                                                                 >  3:110759/1‑241 (MQ=255)
                                                                                 ccATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTa                                                                                                                                                                                                                                                           <  4:86760/202‑1 (MQ=255)
                                                                                  cagcGCTACCGGGAAAGAAACTAGTGAATNNNTGGGCTGGGAGTTTGTTCACTTTATTCTCGACGTCGCCAGCGGGAAGAATAAAACCTTCTCGGAGCAATGGGGGCTGCATAACCAGCTGGCTGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaaa                                                                                                                                                                                 <  3:143162/272‑1 (MQ=255)
                                                                                  cATCGCTACCGGCGAAGAAACTATTGGAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaaa                                                                                                                                                                                 <  1:110759/275‑1 (MQ=255)
                                                                                               gaagaaACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                                                                                                                                                    <  2:134302/275‑1 (MQ=255)
                                                                                               gaagaaACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGGGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                                                                                                                                                    <  4:245350/275‑1 (MQ=255)
                                                                                                              gagtggggCTGGAAGTTGTTCCACTTTAGTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCACGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATg                                                                                                                                                      <  2:397711/269‑1 (MQ=255)
                                                                                                                   gggCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTTGTTCAAGCTTATAATTTGAGCAGACAACAGGTCTTAAACAGACAGGTAATCATTCCCGc                                                                                                                                                 >  4:382340/1‑274 (MQ=255)
                                                                                                                   gggCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAACAGGAAAGTAAATCATGCCCgat                                                                                                                                                >  4:417504/1‑273 (MQ=255)
                                                                                                                       tGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCGCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCt                                                                                                                                                <  2:417505/271‑1 (MQ=255)
                                                                                                                              tgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAATCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTTCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAATGCCCCTAAGTCCATTAATAATGCATTTGCATTTCCCTAGTTCAAGCTTACAGTTTGAGCAGAAATCAGGTTGTAAAAGGACCGTGAACCATCCCCgtaaatccccggc                                                                                                                                      >  4:305784/1‑262 (MQ=255)
                                                                                                                               gttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgctc                                                                                                                                     <  2:382341/274‑1 (MQ=255)
                                                                                                                                     ctgtcttttcaaCGTCGCCAGCGGAAAGTATAAAACCCTCACGGTGCAATGGGGGCTGCAAAACCATCTGGCGGTGTTTAACCCGGAACCGTGTACCTGATTTCTTTTGGAACCTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCAACGTCCATTAATAATGCAGTTGCATTACCGTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCt                                                                                                                              <  4:118686/264‑1 (MQ=255)
                                                                                                                                               gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCa                                                                                                                    >  4:10429/1‑275 (MQ=255)
                                                                                                                                               gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCa                                                                                                                    >  4:130971/1‑275 (MQ=255)
                                                                                                                                                    cGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAg                                                                                                                                                                                                                                                                        >  6:182428/1‑122 (MQ=255)
                                                                                                                                                    cGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAg                                                                                                                                                                                                                                                                        <  5:182428/122‑1 (MQ=255)
                                                                                                                                                    cGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGaa                                                                                                               >  4:240602/1‑275 (MQ=255)
                                                                                                                                                        aGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCt                                                                                                                                                                                                                                                                           <  4:542284/115‑1 (MQ=255)
                                                                                                                                                               agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                                                                                                                                                   >  7:40583/1‑116 (MQ=255)
                                                                                                                                                               agaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                                                                                                                                                   <  8:40583/116‑1 (MQ=255)
                                                                                                                                                                          ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAt                                                                                                                                                                                                                                                          >  6:500217/1‑114 (MQ=255)
                                                                                                                                                                          ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAt                                                                                                                                                                                                                                                          <  5:500217/114‑1 (MQ=255)
                                                                                                                                                                              cGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                                                                                                                                                >  7:838316/1‑136 (MQ=255)
                                                                                                                                                                              cGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATAGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                                                                                                                                                <  8:838316/136‑1 (MQ=255)
                                                                                                                                                                                 aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTTCCTCGCAAAACATGGCACTAAGGCCATGt                                                                                                                                                                                                                                                                >  8:936490/1‑101 (MQ=255)
                                                                                                                                                                                     aTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                                                                                                                                             >  7:517460/1‑148 (MQ=255)
                                                                                                                                                                                        ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                                                                                                                                                   >  6:996114/1‑91 (MQ=255)
                                                                                                                                                                                        ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCa                                                                                                                                                                                                                                                                   <  5:996114/91‑1 (MQ=255)
                                                                                                                                                                                         gggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                                                                                                                                             <  6:438174/144‑1 (MQ=255)
                                                                                                                                                                                         gggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                                                                                                                                             >  5:438174/1‑144 (MQ=255)
                                                                                                                                                                                           gCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaa                                                                                                                                                                                  >  2:75320/1‑169 (MQ=255)
                                                                                                                                                                                              gCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTg                                                                                                                                                                                                                                                        <  7:5101/96‑1 (MQ=255)
                                                                                                                                                                                              gCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTg                                                                                                                                                                                                                                                        >  8:5101/1‑96 (MQ=255)
                                                                                                                                                                                               cATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                                                                                                                                              <  4:75319/169‑1 (MQ=255)
                                                                                                                                                                                                 tAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAt                                                                                                                                                                                                                                                          <  7:82942/91‑1 (MQ=255)
                                                                                                                                                                                                 tAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAt                                                                                                                                                                                                                                                          >  8:82942/1‑91 (MQ=255)
                                                                                                                                                                                                     cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                                                                                                                                        >  5:55513/1‑121 (MQ=255)
                                                                                                                                                                                                     cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                                                                                                                                        <  6:55513/121‑1 (MQ=255)
                                                                                                                                                                                                       gCTGGCGGGGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGa                                                             <  2:256604/274‑1 (MQ=255)
                                                                                                                                                                                                           gCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                                                                                                                                       <  8:517460/148‑1 (MQ=255)
                                                                                                                                                                                                           gCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCt                                                        >  4:341295/1‑275 (MQ=255)
                                                                                                                                                                                                            cGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCtt                                                       <  2:341296/275‑1 (MQ=255)
                                                                                                                                                                                                            cGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtt                                                                                                                                                                                                                                                               >  6:143593/1‑75 (MQ=38)
                                                                                                                                                                                                                    aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                                                                                                                                             <  7:841714/117‑1 (MQ=255)
                                                                                                                                                                                                                    aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                                                                                                                                             >  8:841714/1‑117 (MQ=255)
                                                                                                                                                                                                                     aCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                                                                                                                                       >  5:521799/1‑122 (MQ=255)
                                                                                                                                                                                                                     aCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                                                                                                                                       <  6:521799/122‑1 (MQ=255)
                                                                                                                                                                                                                        cGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGccc                                                                                                                                                                                                                                              >  8:445743/1‑80 (MQ=255)
                                                                                                                                                                                                                        cGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGccc                                                                                                                                                                                                                                              <  7:445743/80‑1 (MQ=255)
                                                                                                                                                                                                                           cACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                                                                                                                                                   >  3:136291/1‑88 (MQ=255)
                                                                                                                                                                                                                           cACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATTCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAgccagg                                        >  2:294794/1‑275 (MQ=255)
                                                                                                                                                                                                                             ccGGTTACCTGGTTTCTTTTGGAACGTCCCTCGCCAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAgccaggc                                       <  4:536277/274‑1 (MQ=255)
                                                                                                                                                                                                                             ccGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAgccaggc                                       <  2:130972/274‑1 (MQ=255)
                                                                                                                                                                                                                              cGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                                                                                                                                       >  5:506637/1‑113 (MQ=255)
                                                                                                                                                                                                                              cGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                                                                                                                                       <  6:506637/113‑1 (MQ=255)
                                                                                                                                                                                                                               ggTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaata                                                                                                                                                                                                                               <  1:136291/88‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTg                                                >  2:62768/1‑257 (MQ=255)
                                                                                                                                                                                                                                         tttCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAgc                                            <  4:62767/257‑1 (MQ=255)
                                                                                                                                                                                                                                              tttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCAAATGCATTattaa                                                                                                                                                                                                             >  8:951731/1‑87 (MQ=255)
                                                                                                                                                                                                                                              tttGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCAAATGCATTattaa                                                                                                                                                                                                             <  7:951731/91‑5 (MQ=255)
                                                                                                                                                                                                                                                            cGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCACCGCGCCAGTGCGTGAGGCCTTAAAACTTAAGCCCGGCCAGGCCAGCATTCATTACGAAATTCTGCCtgg       >  4:365452/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 aaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                                                                                                                                            <  8:917918/73‑1 (MQ=37)
                                                                                                                                                                                                                                                                 aaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAg                                                                                                                                                               >  6:468980/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                 aaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAg                                                                                                                                                               <  5:468980/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc                                                                                                                                 >  5:109466/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  aCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAGCCAGGCCAGGACAGCATTCATTACGAAATTCTGCCTGGTGGGc  <  2:240603/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aTGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATaa                                                                                                                                                                                              <  7:1003650/84‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aTGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATaa                                                                                                                                                                                              >  8:1003650/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                    aTGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct                                                                                                                                    >  4:10769/1‑142 (MQ=255)
                                                                                                                                                                                                                                                                        ccTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg                                                                                                                                <  2:10769/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCt                                                                                                                                                                           >  6:803838/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                         cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCt                                                                                                                                                                           <  5:803838/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACt                                                                                                                         >  5:456737/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          ttAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                                                                                                                                            >  5:790620/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                          ttAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                                                                                                                                            <  6:790620/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                         
TGGCGGTACCCGTCATTAAAATGGCAACCCGCACCGAGCTGGCGAACCGCTGGTTTGATTTAATGGATATTAATGCGGGCACCATCGCTACCGGCGAAGAAACTATTGAAGAGGTGGGCTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTGACCACTGAATCAAAGGTCACGATACGCGGACAAACAACTATCCCCGCGCCAGTGCGTGAGGCCTTAAAACTGAAGCCAGGCCAGGACAGCATTCATTACGAAATTCTGCCTGGTGGGC  >  NC_000913/3276621‑3277152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: