Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 78,830 T→A 5.6% intergenic (+31/+18) setA → / ← leuD broad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091378,8300TA 5.6% 318.1 / 5.7 90intergenic (+31/+18)setA/leuDbroad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base A (3/2);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

TGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAG  >  NC_000913/78565‑79096
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
tGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTCCTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCCTTATTTTTAACCGCAAAGGTTAAAGTCGTTTGATGCCGTGGACGATTGCGGAAAGCCCGGGCATTTGCCCggg                                                                                                                                                                                                                                                                   >  4:367113/1‑275 (MQ=255)
tGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCggg                                                                                                                                                                                                                                                                   >  2:228427/1‑275 (MQ=255)
  tATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCa                                                                                                                                                                                                                                                                 >  4:384808/1‑275 (MQ=255)
    tttATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTCAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAg                                                                                                                                                                                                                                                               <  2:94670/275‑1 (MQ=255)
    tttATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTAAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAg                                                                                                                                                                                                                                                               <  2:427770/275‑1 (MQ=255)
    tttATCGGCATTGTTGCGGGTATTGGGATGCTATAGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAg                                                                                                                                                                                                                                                               <  2:376249/275‑1 (MQ=255)
                  tgCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATa                                                                                                                                                                                                                                                 <  2:478868/275‑1 (MQ=255)
                           tGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAt                                                                                                                                                                                                                                                                            >  2:176717/1‑239 (MQ=255)
                           tGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgt                                                                                                                                                                                                                                        >  4:447469/1‑275 (MQ=255)
                               aTGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGa                                                                                                                                                                                                                                                                        <  4:176716/239‑1 (MQ=255)
                                tGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGttt                                                                                                                                                                                                                                   <  2:347854/275‑1 (MQ=255)
                                  cTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAATGGTTAACGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCCAGGGGATTAATTCATAAACGCAGGTTGTTTTg                                                                                                                                                                                                                                 >  2:122563/1‑275 (MQ=255)
                                    atcgTTTCAGGATCTAATGTCTGGAAGAGCGAGGGCAGCCACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATCCTGGCTGGCGCTATTCCGGGAGCAATTGCACAAAGTTGGTGGCACTTTGCTGTCTACTGGTTAATTTCGGTTATTGCAGTTGTCGCATTTTTTTTAACCGCATCGGTTACCGACGTTTTATGATGTGGACGATAGCGGACAGCCCGGTCCGTTGACCGGGCAAGGGGATTATTGCATATACGACGGTTNGTTTgat                                                                                                                                                                                                                               >  4:512177/4‑273 (MQ=255)
                                                aggggacggttgAAGAGGTGGGGTCGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGGAATTCTGGCCGGCGTTATTACGGGAGCAATTGCGCAAAGTGTGGAGCACTTTGCGGTCTACTGGGTAATTGCGGTTATTTCGGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGc                                                                                                                                                                                                                    <  4:311707/263‑1 (MQ=255)
                                                                   ggggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTTCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCATGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAg                                                                                                                                                                                                                        >  2:311708/1‑251 (MQ=255)
                                                                       cAGCTACCACCTTATTTCATAACATGATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:233992/275‑1 (MQ=255)
                                                                       cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:228426/275‑1 (MQ=255)
                                                                       cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  2:447470/275‑1 (MQ=255)
                                                                       cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  2:343667/275‑1 (MQ=255)
                                                                       cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  2:367114/275‑1 (MQ=255)
                                                                       cAGCTAACTCCTTTTTTGCTAACATAATTTCTTCCGGGGGAATTCTGGCTGGGGTGTTTCAGGGAGAAATTGCAAAAAGTTGTGGGAAATTTGCTGTGTAATGGGTAATGGCGGTTATTTCGGTTGTCGAATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:483246/275‑1 (MQ=255)
                                                                        atgtACCACCTTTTTTAATAAAAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:446065/271‑1 (MQ=255)
                                                                                  ttATTTACGAACATTATTTCTACCGGGGTAATTCTTGGTGGGGATATTAAGGGGGTAAATGAAAAAAGTAGGGGCAACTTTGCTGTCTACAGGGTAATTGCGGTTATTTCTGTTTTCGCATTATTTTTAACCGCAAAGTTTAAAGACGTTTGATGACGTTGACGATAGCGGAAGGCCCGGTCATTTGACCGGGCAAGGAGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATTGCGTCGTCGTGCTGCAAGGTAAGCCCAATa                                                                                                                                                                                  <  4:507599/274‑1 (MQ=255)
                                                                                  gtATTTACTATCACTATTTTTGCGCGGGTAATTTTTGCTGGCGTTTTTCAGGGCGCAATTGCACAAAGTTGGGGGCATGTTGGTTTTTCCTACGGAATTGCGGTTATTGCTGTTGTCGCATTATGTTTAGCCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATa                                                                                                                                                                                  <  4:122562/273‑1 (MQ=255)
                                                                                       tACTAACATTAGGTCTACCGGGGTAAGTCTGGCTGTCGTTATTCAGGTAGCAATTGAACAAAGTTGGGGGCAAGTTGCTGTCTACTGGGTAATTGCTGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTTATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAATCGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTcc                                                                                                                                                                            <  2:96556/275‑1 (MQ=255)
                                                                                          tAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGTTGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCGTAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCa                                                                                                                                                                                                >  2:248909/1‑252 (MQ=255)
                                                                                          tAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTGTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGa                                                                                                                                                                         <  4:339863/275‑1 (MQ=255)
                                                                                          tAACAGTAGTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGa                                                                                                                                                                         <  2:368856/275‑1 (MQ=255)
                                                                                           aaCAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGAc                                                                                                                                                                        >  2:455571/1‑275 (MQ=255)
                                                                                              aGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGTTGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCGTAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:248908/252‑1 (MQ=255)
                                                                                                        ccGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGAATTATTTTTAACCGCAAAGGTTAAAGAAGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAACGGGAAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGc                                                                                                                                                            >  1:47640/1‑274 (MQ=255)
                                                                                                            ggTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTgg                                                                                                                                                       <  2:384809/275‑1 (MQ=255)
                                                                                                                         ggcGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                                                                                                                                                         <  7:5041/148‑1 (MQ=255)
                                                                                                                                ttCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGt                                                                                                                                   <  4:455570/275‑1 (MQ=255)
                                                                                                                                     ggAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCgg                                                                                                                                                                                                                                                                    >  4:539441/1‑141 (MQ=255)
                                                                                                                                         cAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCaa                                                                                                                                                                                                                                                                <  2:539442/141‑1 (MQ=255)
                                                                                                                                            ttGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAt                                                                                                                                                                                                                                                      >  6:170521/1‑148 (MQ=255)
                                                                                                                                                 caAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAt                                                                                                                                                                                                                                                                            >  8:755459/1‑121 (MQ=255)
                                                                                                                                                 caAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAt                                                                                                                                                                                                                                                                            <  7:755459/121‑1 (MQ=255)
                                                                                                                                                    aGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCtt                                                                                                                                                                                                                              >  2:257774/1‑164 (MQ=255)
                                                                                                                                                        gggggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATa                                                                                                                                                                                                                          <  4:257773/164‑1 (MQ=255)
                                                                                                                                                         ggggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCCAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGCCGCCGTGCTGCAATGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCCGCCCTGGTAACCCCGTAGGTTTtctcccacgt                                                                                                           >  4:337484/1‑268 (MQ=255)
                                                                                                                                                              aCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATTACGTGGACGATAGCGGAAAGCCCGGTCATTTTACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCCATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCATGGGCCGCGGCATGACACGATGTTAAAGCGATAGGTTTTCTCCCCCCCTTTCAc                                                                                                     >  4:171276/1‑275 (MQ=255)
                                                                                                                                                              aCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTGCTCTCCCGCTTTCAc                                                                                                     >  1:21500/1‑275 (MQ=255)
                                                                                                                                                                 ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                                                                                                                                                         <  5:300387/108‑1 (MQ=255)
                                                                                                                                                                 ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                                                                                                                                                         >  6:300387/1‑108 (MQ=255)
                                                                                                                                                                          aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAc                                                                                                                                                                                                                                                                       >  5:922467/1‑101 (MQ=255)
                                                                                                                                                                          aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAc                                                                                                                                                                                                                                                                       <  6:922467/101‑1 (MQ=255)
                                                                                                                                                                          aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCCTAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCa                                                                                                                                                                           >  2:268093/1‑193 (MQ=255)
                                                                                                                                                                           cggggTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGTGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTc                                                                                         <  2:263514/272‑1 (MQ=255)
                                                                                                                                                                              ggTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGAcc                                                                                                                                                                       <  4:268092/193‑1 (MQ=255)
                                                                                                                                                                                 aaTTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATc                                                                                   >  4:386976/1‑274 (MQ=255)
                                                                                                                                                                                   ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATCGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGGGATAGGTTTTCTCTACCGCTTTCACCTCTTGCGCTTCCAGATCCa                                                                                 >  2:89587/1‑274 (MQ=255)
                                                                                                                                                                                   ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATgg                                                                                                                                                                                                               <  5:170521/148‑1 (MQ=255)
                                                                                                                                                                                   ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCa                                                                                                                                                                                                >  1:76801/1‑163 (MQ=255)
                                                                                                                                                                                   ttGCGGTTATTTCTGTTGTCGCACTATTTTTAAACGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGCTTAATTCATAAACGTAGGTTGTTTTGCTTCATAAGTGGTATTGGCGTCGTCGTGCTGCTAGGTAAGCCACATTCTGTCCAGACCGGTCATCATGCCGTTGCGGCGGAAGGCATCACAGGTAAACGGATAGGTTATCACCCCCGCATTCCACTCTTGCGCCTcaggaaccc                                                                                 >  4:27390/1‑266 (MQ=255)
                                                                                                                                                                                       ggTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  3:76801/163‑1 (MQ=255)
                                                                                                                                                                                              tCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATggg                                                                     >  2:169774/1‑274 (MQ=255)
                                                                                                                                                                                                   tgtCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGAAAGGCATCGATGGTAAAGAGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCTACGTCGCAATGAATccc                                                                 >  1:147590/1‑274 (MQ=255)
                                                                                                                                                                                                      cGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                                                                                                                                           >  6:272988/1‑69 (MQ=255)
                                                                                                                                                                                                      cGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                                                                                                                                           <  5:272988/69‑1 (MQ=255)
                                                                                                                                                                                                       gCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGAt                                                                                                                                      >  2:211636/1‑201 (MQ=255)
                                                                                                                                                                                                         attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCa                                                                                                                                                                                                                                           >  6:971239/1‑98 (MQ=255)
                                                                                                                                                                                                         attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCa                                                                                                                                                                                                                                           <  5:971239/98‑1 (MQ=255)
                                                                                                                                                                                                           tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTa                                                                                                                                  <  4:211635/201‑1 (MQ=255)
                                                                                                                                                                                                           tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCgg                                                                                                                                                                                                                                                                    >  5:618390/1‑71 (MQ=255)
                                                                                                                                                                                                           tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCgg                                                                                                                                                                                                                                                                    <  6:618390/71‑1 (MQ=255)
                                                                                                                                                                                                                       cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                <  8:884902/79‑1 (MQ=255)
                                                                                                                                                                                                                       cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                <  8:695662/79‑1 (MQ=255)
                                                                                                                                                                                                                       cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                >  7:884902/1‑79 (MQ=255)
                                                                                                                                                                                                                       cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                >  7:695662/1‑79 (MQ=255)
                                                                                                                                                                                                                         aaGGTTAAAGACGTTTGATGACTTGGACGATAGAGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCAAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCaa                                          <  2:474285/275‑1 (MQ=255)
                                                                                                                                                                                                                              tAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGt                                     <  2:337485/275‑1 (MQ=255)
                                                                                                                                                                                                                                   aCGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTc                                <  2:386977/275‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTCTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCaaa                                         >  2:291645/1‑265 (MQ=255)
                                                                                                                                                                                                                                      tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTa                                                                                                                                                                                           >  4:370490/1‑117 (MQ=255)
                                                                                                                                                                                                                                      tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCCCCTCTTGCGCTTCCAGATCCAAGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCTCAAACAGTTCGTCCAc                             >  1:113639/1‑275 (MQ=255)
                                                                                                                                                                                                                                      tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCAc                             >  1:33419/1‑275 (MQ=255)
                                                                                                                                                                                                                                      tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACAAGCGCAAACAGTTCGTCCAc                             >  2:319287/1‑275 (MQ=255)
                                                                                                                                                                                                                                      gttGATGTCGTGGATGATAGCTGAAGGCCCGGTCATTTGACCGGGCCAGGGGATTAATTCATAAACGCAGGTTTGTTTGCTTCATAAGAGGCAATGGCTTCGTCGGGCTGCGAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCAc                             <  4:89586/274‑1 (MQ=255)
                                                                                                                                                                                                                                          atgaCGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGcc                                                                                                                                                                                       <  2:370491/117‑1 (MQ=255)
                                                                                                                                                                                                                                          atgaCGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                                                                                                                                                     >  7:1128162/1‑39 (MQ=37)
                                                                                                                                                                                                                                          atgaCGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                                                                                                                                                     <  8:1128162/39‑1 (MQ=37)
                                                                                                                                                                                                                                            gaCGTGGACGTTAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCAGAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGc                       <  2:512178/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  gACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCg                   <  2:155097/273‑1 (MQ=255)
                                                                                                                                                                                                                                                               aGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGgc    >  2:150498/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               aGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGCTGTTTTGCTTCATAAGCGCTAATGGCGTCGTAGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGTCATCGCTGGTAAAGCGATAGATTTTCTCCCGCGCTTTCACTTCTTCCGCTTGCAGACCCAAGTCGAAATGTATCCCTGGTTACGCTTGCACCACCGCAAGCGTTTCTCCCAATACTGCACCGCTTAATTTCCCCgta    >  4:136568/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                 cccGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGgcag  >  2:355335/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                                                                                                                                               >  6:792201/1‑52 (MQ=255)
                                                                                                                                                                                                                                                                   cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                                                                                                                                               <  5:792201/52‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAGGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                                                                                                                                                        <  7:854996/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAGGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                                                                                                                                                        >  8:854996/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                   cGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAg                                                                                                                               >  7:527078/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg                                                                                                                                    >  6:38876/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                     gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg                                                                                                                                    <  5:38876/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
TGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAG  >  NC_000913/78565‑79096

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