Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,915 A→G 6.1% H26H (CAT→CAC) ‡ cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,9150AG 6.1% 396.4 / 11.2 125H26H (CAT→CAC) ‡cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base A (59/55);  new base G (4/4);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-02

Only 100 of 127 total aligned reads displayed.

TCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTT  >  NC_000913/2873644‑2874187
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                
   ggCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGGAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCg                                                                                                                                                                                                                                                                            <  2:471423/275‑1 (MQ=255)
            cgcgGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAACACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCCCCCCGGCGACCGTTTCTTTCCCGGACCCGGAGAGGCCCGTCAACCACAGCACAACCCCACGAtacccctagcgtat                                                                                                                                                                                                                                                                    >  4:170427/1‑261 (MQ=255)
            cgcgGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGGGTAg                                                                                                                                                                                                                                                                    >  4:86311/1‑274 (MQ=255)
               ggTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                                                                                                                                                <  4:290400/275‑1 (MQ=255)
                               tGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTg                                                                                                                                                                                                                                                 >  2:404632/1‑274 (MQ=255)
                                   gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGAACGTGGTGTAGCTCGCGTTGTTGCACAGGTACCgg                                                                                                                                                                                                                                            >  2:460481/1‑275 (MQ=255)
                                    tCAGCACCACCAGGTCGGCTTCAACCAGCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATGGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCgg                                                                                                                                                                                                                                            <  4:397381/274‑1 (MQ=255)
                                        caccacCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGAt                                                                                                                                                                                                                                                                                 <  2:210185/233‑1 (MQ=255)
                                                gTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCCTCGCTCAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTCTGTAACGCCTCCTCCAGCCCCCCGCCGACCGTTGATTTACCGGACCCGGTGCGACCGGTAAACCCCCGCACTACATCACCATGGCCGTGGGGTCACCCGGGTTGTCGCAAGCTGACCCGTTGGCCTTGCCAg                                                                                                                                                                                                                               >  4:173569/1‑275 (MQ=255)
                                                gTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAg                                                                                                                                                                                                                               >  2:277120/1‑275 (MQ=255)
                                                  ccGGCTCCAACCATTAACTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGACCGGCATCGCTAAAACCGATATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATCCGTACTGCCGCCGCCCTTATGTAACGCCTCCCCCACCGCACCGCCGCCCTTTGATCTACCTGACGCGGACCGGCCTGCACACACCACCACTACACCTCACTCACCGTGTGATCGCTCCGCTTGTTGCACAGTGACCGCATGGCTagcccacac                                                                                                                                                                                                                             >  4:95208/1‑267 (MQ=255)
                                                     gCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTGTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACAGTTGATTTACCGGACCCGGAGAGGCCGGTACACCAAAGCACTACACCACGCTTACCGAGGTTTAGCCCGCGTTGTTGGCAAGTGACCGGCAGGCCATACCCgacca                                                                                                                                                                                                                           >  1:130121/1‑272 (MQ=255)
                                                                    ttCCCCACGTCACCCACGCGACGGATATTCTCTTTACGATCGGCATCGCTACAACCGGGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAg                                                                                                                                                                                                           <  4:404631/275‑1 (MQ=255)
                                                                     tCGCCACTTCACCGACGCGACGGATATTCTCTGTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAg                                                                                                                                                                                                           <  2:320982/274‑1 (MQ=255)
                                                                      cGCCCCTGCACCTCCGCGACGTATATTCTCTTTACGATCGGCATCGCTAAAACCGCGATCGCTGCATAAGCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGATTTTATGTAACGCCTCCTCCAGCGCCCCTGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcg                                                                                                                                                                                                         <  1:104593/275‑1 (MQ=255)
                                                                      cGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGCTTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACAGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                          >  2:105028/1‑274 (MQ=255)
                                                                         cccTTCCCCGCCGCGACTGATATTCTCTTNNCGATCGGCATCGCTAAATCCGAGATCGCGGNAGAATCCGTGGGGAACATCGTCCACATCCAGCAGATACGTACTGAGGGCGAGTGTAGGTAACGCCGACTCCAGCNCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                                                                                                                                                      <  3:73371/273‑1 (MQ=255)
                                                                         ccattaaCCGACCCAACGGATATTCTCTTTCCGATCGGAATCGCTAAAACCGGGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                                                                                                                                                      <  4:277119/269‑1 (MQ=255)
                                                                         cACTTCACCGAGGCGAAGGATAATCTCTTTACGATCGGCATCGGTACAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                                                                                                                                                      <  1:75031/275‑1 (MQ=255)
                                                                         cACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGAAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                                                                                                                                                      <  2:173570/275‑1 (MQ=255)
                                                                                              atTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTTTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                                                                                                                                                 >  2:21789/1‑275 (MQ=255)
                                                                                                   ctTTACGATCGGCATCGCTTAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgcg                                                                                                                                                                            <  2:150286/275‑1 (MQ=255)
                                                                                                      tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGGCATGCAGCGCCATTCTTTATTCCCCAGCAAATCgcccccgc                                                                                                                                                                         >  1:149777/1‑275 (MQ=255)
                                                                                                       aCGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGcc                                                                                                                                                                        <  2:81818/275‑1 (MQ=255)
                                                                                                                    cTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGttgttg                                                                                                                                                                                                                                                        >  2:307053/1‑182 (MQ=255)
                                                                                                                        aaCCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                                                                                                                                                    <  4:307052/182‑1 (MQ=255)
                                                                                                                         aCCGGGATGCCTGCAGAACCCCGGTCGAAAATGGTCGCCATCCCGCAGATACGTACGGACGCCGAGTTTATGTAGCGCCCCCCCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcg                                                                                                                                                      <  4:460480/275‑1 (MQ=255)
                                                                                                                                           tCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAACGTGGCGACGAACCAGTGCATTCa                                                                                                                                     >  2:388401/1‑274 (MQ=255)
                                                                                                                                             cGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTTGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCACCAAATCGCGCCCGCCCCAGCGCGGAAAGTGGCGACGAACAAGAACATTCAAt                                                                                                                                   >  4:349050/1‑274 (MQ=255)
                                                                                                                                                ggCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCCGTGCGGAAAGTGGCGACGAACCAAAGCATTAAATTccc                                                                                                                               >  1:3211/1‑274 (MQ=255)
                                                                                                                                                ggCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGCAAGTGGCGACGAACCAGAGCATTCAATTCCa                                                                                                                               >  2:124757/1‑275 (MQ=255)
                                                                                                                                                ggCGAACATTGTCGCCACCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCACCAGCGCCCCGACGACCGTTGATATACCGGACCCGGAGAGGCCGCTCAACCACCGCACTACACCCCGATGACCGGGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGTCTATTCTAGACCACGTTTTCGTCGTGCAGCCCCATTATTTATCCCCCAGAACATCCCCCGCTCCCCCGTGCTGCAAGTGTCCACGAACCACCACattccattcc                                                                                                                                >  4:438730/1‑274 (MQ=255)
                                                                                                                                                            cGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCAc                                                                                                                   <  2:350452/275‑1 (MQ=255)
                                                                                                                                                                aTCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGa                                                                                                                >  2:6753/1‑274 (MQ=255)
                                                                                                                                                                  ccAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttg                                                                                                                                                                                                                                                           >  8:581704/1‑133 (MQ=255)
                                                                                                                                                                      cagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGc                                                                                                                                                                                                                                                                   <  6:105342/121‑1 (MQ=255)
                                                                                                                                                                      cagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGc                                                                                                                                                                                                                                                                   >  5:105342/1‑121 (MQ=255)
                                                                                                                                                                        gATACGTACGGACGCCGCGTTTATGTAAAGCCTCCTCCGGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCAAAGCACTACCCCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGAt                                                                                                        <  4:112584/274‑1 (MQ=255)
                                                                                                                                                                               aCTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                         >  5:284979/1‑106 (MQ=255)
                                                                                                                                                                               aCTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                         <  6:284979/106‑1 (MQ=255)
                                                                                                                                                                                  gACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATgg                                                                                                                                                                                                                                        <  8:1091748/136‑1 (MQ=255)
                                                                                                                                                                                  gACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATgg                                                                                                                                                                                                                                        >  7:1091748/1‑136 (MQ=255)
                                                                                                                                                                                      ccGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                          >  5:493567/1‑114 (MQ=255)
                                                                                                                                                                                      ccGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                          <  6:493567/114‑1 (MQ=255)
                                                                                                                                                                                                gTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTa                                                                                                                                                                                                                                     >  6:1069118/1‑125 (MQ=255)
                                                                                                                                                                                                gTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTa                                                                                                                                                                                                                                     <  5:1069118/125‑1 (MQ=255)
                                                                                                                                                                                                   aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGaa                                                                                                                        >  2:317910/1‑231 (MQ=255)
                                                                                                                                                                                                    cGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCgtggt                                                                                                                                                                                                                                                                        >  5:18378/1‑86 (MQ=255)
                                                                                                                                                                                                    cGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCgtggt                                                                                                                                                                                                                                                                        <  6:18378/86‑1 (MQ=255)
                                                                                                                                                                                                    cGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGCTGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTGCGAATGTCCTGAATTCAGATGGCCCAGCCGTTGCCGGGCTAACTggac                                                                            >  2:394452/1‑272 (MQ=255)
                                                                                                                                                                                                    cGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTc                                                                           >  4:205976/1‑275 (MQ=255)
                                                                                                                                                                                                    cGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAAACACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTTCCTGGCTAACTTGccc                                                                           >  4:50384/1‑273 (MQ=255)
                                                                                                                                                                                                      cctcctCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                          >  5:709726/1‑98 (MQ=255)
                                                                                                                                                                                                      cctcctCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                          <  6:709726/98‑1 (MQ=255)
                                                                                                                                                                                                       ctcctcCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCa                                                                                                                    <  4:317909/231‑1 (MQ=255)
                                                                                                                                                                                                        tcctccGGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgt                                                                                                                                                                                                                                                             >  7:563235/1‑93 (MQ=255)
                                                                                                                                                                                                        tcctccGGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgt                                                                                                                                                                                                                                                             <  8:563235/93‑1 (MQ=255)
                                                                                                                                                                                                        ccctccAGGGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTg                                                                        <  2:531998/273‑1 (MQ=255)
                                                                                                                                                                                                         cctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGtcc                                                                       >  2:514000/1‑272 (MQ=255)
                                                                                                                                                                                                         cctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACATTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGaa                                                                                                                        >  4:167065/1‑225 (MQ=255)
                                                                                                                                                                                                          ctccAGCGCCCCGGCGACCGTTGAGTTACCGGACCCGGAGAGGCCGGTAAACAACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGAATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCAc                                                                     <  1:94949/275‑1 (MQ=255)
                                                                                                                                                                                                             cAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACATTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCa                                                                                                                    <  2:167066/225‑1 (MQ=255)
                                                                                                                                                                                                                cgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                         <  8:217345/73‑1 (MQ=255)
                                                                                                                                                                                                                cgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATAc                                                                >  1:110292/1‑274 (MQ=255)
                                                                                                                                                                                                                 gcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTTGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATTCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGCAAGTGGCGACGAACCACACCCTTCAATGCCAGTTCGAAAGCATTGAATTCAAATGGCCGCGCAGTTCCCTGTCTAACTGGCCCGTGCACAATACCg                                                              >  4:127624/1‑275 (MQ=255)
                                                                                                                                                                                                                   cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                         <  7:173646/86‑1 (MQ=255)
                                                                                                                                                                                                                   cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                         >  8:173646/1‑86 (MQ=255)
                                                                                                                                                                                                                     cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAg                                                                                                                                                                                                                               >  8:1004809/1‑110 (MQ=255)
                                                                                                                                                                                                                     cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                                                                                                                                              <  7:1139964/143‑1 (MQ=255)
                                                                                                                                                                                                                     cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                                                                                                                                              >  8:1139964/1‑143 (MQ=255)
                                                                                                                                                                                                                       gCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGCGCATTCAATTCCAGTTCGAATGCCCTGAATTCAGATGGCCCAGCAGTTGCCTGGCTAACTGGCTCGTTCCCCCAACCGGcccc                                                         >  4:249820/1‑274 (MQ=255)
                                                                                                                                                                                                                                      cctgACCCGGGCGGGCCGGTGAACCAAAGGCATACAGCACGGTGACCGTGGTGTAGCTCGCGTTGTTGCACGGTGACCGGATGGCTATGGCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCt                                          <  4:124756/271‑1 (MQ=255)
                                                                                                                                                                                                                                      ccGGACCCGGGGAGGCCGGTAAACCACAGCACTACCCCCCGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTc                                         <  4:504445/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGCACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACCTTGCTCAGGCGATCGATaaa                           >  2:450897/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     ccGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                                                                                                                                                  <  6:909456/91‑1 (MQ=255)
                                                                                                                                                                                                                                                     ccGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                                                                                                                                                  >  5:909456/1‑91 (MQ=255)
                                                                                                                                                                                                                                                        gTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAAt                        >  2:144096/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               acaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtttt                                                                                                                                                                                                                     <  5:571073/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCag                                                                                                                                             <  6:377196/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                       acacCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt                                                                                                                                       >  6:1044822/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       acacCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt                                                                                                                                       >  8:708961/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                       acacCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                 >  5:613930/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                       acacCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGAGACGaa                                                                                                                                                 <  6:613930/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                        caccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGGCAGCGATAAAAATCAGCCCACCCgcacccg       >  4:494839/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                        aacccggATGACCGGGGTGTGGGTCTCGTTGTTGCACAGGGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCAAGGGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTGGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCGCGGTCACATTGGTCAGGCGATCGATAAAAATCAGCCCACCCGTCAcc        <  4:272469/268‑1 (MQ=255)
                                                                                                                                                                                                                                                                         accacGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCa                                                                                                                                                                                                                >  7:896780/1‑73 (MQ=38)
                                                                                                                                                                                                                                                                         accacGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCa                                                                                                                                     <  7:708961/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                   <  5:215148/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa                                                                                                                                    >  5:12759/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                           cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAAt                                                                                                                                                                                  <  6:168366/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGc                                                                                                                                          >  6:1018074/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                           cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGc                                                                                                                                          <  5:1018074/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                            acGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTc                                                                                                                                                                                                                 >  7:1099372/1‑69 (MQ=255)
                                                                                                                                                                                                                                                                            acGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt                                                                                                                                  >  5:1093244/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             cGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGAttt  >  4:68755/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                              gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtt                                                                                                                                                                                                                       <  7:549886/61‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                                                                                                                                                     <  6:651178/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTcc                                                                                                                                >  5:629804/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               aTGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATg                                                                                                                                                                                                              <  8:400557/69‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aTGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATg                                                                                                                                                                                                              >  7:400557/1‑69 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                
TCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTT  >  NC_000913/2873644‑2874187

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