Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,958 G→T 6.5% intergenic (‑45/+28) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9580GT 6.5% 344.8 / 11.3 107intergenic (‑45/+28)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (47/53);  new base T (3/4);  total (50/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-02

Only 100 of 110 total aligned reads displayed.

AACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGT  >  NC_000913/4470684‑4471229
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
aaCACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAgg                                                                                                                                                                                                                                                                                  >  2:218660/1‑275 (MQ=255)
     ggTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACCGCGGTAAAAAGTCTGGCTCATTTTAGGGAGcc                                                                                                                                                                                                                                                                             >  4:102365/1‑275 (MQ=255)
             tAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGANCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGCTAGCTGCCGGTGCATTTTACGCCGAGATAGTTTTGCTCATGANTTCTCCTTTATTACACCGGTAAAAANTCTTGATCATTATATGGAgtcgacctttt                                                                                                                                                                                                                                                                     >  3:92508/1‑265 (MQ=255)
                 gATGTCGCCCACTTGCAGGCCAGGGGGTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                                                                                                                                                                                                 <  4:370243/275‑1 (MQ=255)
                   tGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                                                                                                                                                                                                               <  2:66211/275‑1 (MQ=255)
                    gTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCCCATTTGTGGGAGCCAGCCTTTTTTTACCa                                                                                                                                                                                                                                                              >  4:217496/1‑275 (MQ=255)
                                aaGGCCATCGGGTTCGACGATCGCTTTTACGTTAGCCAGCGACGGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCaattaa                                                                                                                                                                                                                                                  <  2:149912/273‑1 (MQ=255)
                                    ccAGCGGCTTCTACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGc                                                                                                                                                                                                                                             <  2:50908/275‑1 (MQ=255)
                                       ggggtctCGACGATAGTTTTTGCGTTATCCAACGACGTACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAATGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAg                                                                                                                                                                                                                                          <  4:469786/269‑1 (MQ=255)
                                       gCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTTCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCNACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATATGTTTTCTCATGCTTTCTNCTTTATTACAGCGGTAAACAGTCTTGCTACTTATATGGAGCCAGANTNTT‑TTATCAATCCGCaattaatttaccat                                                                                                                                                                                                                                          >  3:37675/1‑269 (MQ=255)
                                       gCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACATCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAg                                                                                                                                                                                                                                          >  2:95349/1‑275 (MQ=255)
                                          gCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                       <  2:46375/275‑1 (MQ=255)
                                          gCTTCGACGAGCGATTTTACGTTATCCCGCGACGGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                       <  1:17172/275‑1 (MQ=255)
                                                       cTTTTACCCTATCCAGCGACTGCCCTTGCTGTGCGTCTACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCTTGATCTGACCTGAAGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGTTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGaa                                                                                                                                                                                                                          <  1:157627/275‑1 (MQ=255)
                                                          ttACGTTATTCAGCGACTGAATGGCTCGGGCAGCCACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                                                                                                                                                                                                        <  2:294815/274‑1 (MQ=255)
                                                                  tCCGGCGACTGCAGGGGCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGGTTTCGGATTTACCGGGATCTGACCGGAGGGGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                                                                                                                                                                                                <  4:502622/274‑1 (MQ=255)
                                                                  tCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGCTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCGGTCGCGATAGTTTTGCGCAGGCTTTCTCCTTTCTTCCAGCGTTCAACAGTCTAGCTCCTTATCTGGAGCCAGACTTTTATTAACAA‑CCGAACTTAATTGGCCGACACCCCACTAGGGGCAAACATTTTTTc                                                                                                                                                                                                               >  4:89254/1‑275 (MQ=255)
                                                                  tCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAATTCTGGCTCATTATAGGGATCCAGACTTTTATTACCAA‑CCGCAATTAATTGGACAGGCCCACATTATGGGAAAACTCTTTTTc                                                                                                                                                                                                               >  4:72303/1‑275 (MQ=255)
                                                                   ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                              <  1:92508/275‑1 (MQ=255)
                                                                   ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGGTTTCGGATTTACCGGGATCTTACCGAAGGTGATGATCATATTGCCTAGATCAACGCCCTGTACGTAAGGACCGATATCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                              <  2:16371/275‑1 (MQ=255)
                                                                            gACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCa                                                                                                                                                                                                                                           >  4:337841/1‑237 (MQ=255)
                                                                                tGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                       <  2:337842/237‑1 (MQ=255)
                                                                                      tgCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGT‑AAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                                                                                                                                                                                             <  2:311046/272‑1 (MQ=255)
                                                                                        ctgCGACGTCGGCGGGTACGTCGCCCGGTTTCGGATTTACCGGGTTCTGACCGGAGGGGATGATCATATTGCCCAGAGCAACGCCCTGTACGTAAGGACCGATAGCGGCCGGTGCATTTTCCGTAGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                                                                                                                                          <  4:430710/272‑1 (MQ=255)
                                                                                                tCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTCCAGCGGTAAAAAGTCTGGCTCATTCACGGGGGTCTGACTTTTATTCCCAA‑CCGGACTTAATTGGCCAGCCCCAATATAAGGGAAAAATTTTTTTCCAAGTATTTGCATTTGCGCCTGATATcct                                                                                                                                                                                  >  4:486674/1‑272 (MQ=255)
                                                                                                                                  ggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                                                                                                                                                                                                               <  5:543853/148‑1 (MQ=255)
                                                                                                                                          cggggggTGATGTTCAGATTTCCCTGTTCACTGCCCTGTCCGTAAGAGCCGATGTGTGCCGGTTCATTTTCCGTCCGGATAGTTTTGCTCATGATTTCTCCTTTATTACAGAGGTAAAAAGTCTGGCTAATTATAGGGAGCCAGGCTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGc                                                                                                                                       <  4:455164/270‑1 (MQ=255)
                                                                                                                                            ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt                                                                                                                                     <  2:217497/275‑1 (MQ=255)
                                                                                                                                             gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTgg                                                                                                                                                                                                                                              >  4:126981/1‑169 (MQ=255)
                                                                                                                                             gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa                                                                                                                                      >  2:137958/1‑273 (MQ=255)
                                                                                                                                                 tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAg                                                                                                                                                                                                                                          <  2:126982/169‑1 (MQ=255)
                                                                                                                                                   atgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTACGTCGCGCTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCGTTATAGGGGGCCAGACTTTTATTACCAA‑ACGCCATTAATTCGCCAACACCACATTATGGTAAAACTCTTTTTCACAGTATTTGCCTTTGAACGGGATACCCTTGGGCCGTTTTCACACGGACAAGCTGGATAAAACCGGTTCGGCCTGGCCGac                                                                                                                              >  4:447204/1‑274 (MQ=255)
                                                                                                                                                        cATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGACAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGTAAACTCTTTTTCTCAGTATTTGCATTTGAGCGCGATATCATTGTCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAGGGCTGATAccg                                                                                                                          >  1:58904/1‑272 (MQ=255)
                                                                                                                                                          tatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACCGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCAACTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAACGCGATATCATTGGCCCTTTTTTGCACGGCCACCCTGGATGAAACCGTTTCGGCATGGcctatagcgtt                                                                                                                        >  4:349445/1‑264 (MQ=255)
                                                                                                                                                                   ttgcAACGCCCGGTACGAAAGGACCGCTAGGTGCCGGTGCATTTTCCGTCGCGATAGTTTTGTTCATGATTTTTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACCTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGGGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                                              <  4:31894/272‑1 (MQ=255)
                                                                                                                                                                    aTCAACGACCTTTACGTAAGGACCGTTAGCTTCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTGATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                                              <  2:486675/274‑1 (MQ=255)
                                                                                                                                                                         cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                                                                                                                                                                                                      >  8:629064/1‑118 (MQ=255)
                                                                                                                                                                         cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                                                                                                                                                                                                      <  7:629064/118‑1 (MQ=255)
                                                                                                                                                                              tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                                        >  6:240981/1‑111 (MQ=38)
                                                                                                                                                                              tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                                        <  5:240981/111‑1 (MQ=38)
                                                                                                                                                                                        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                                        >  7:325521/1‑101 (MQ=38)
                                                                                                                                                                                        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                                        <  8:325521/101‑1 (MQ=38)
                                                                                                                                                                                             aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGcc                                                                                                                                                                                                                                            >  3:82908/1‑123 (MQ=255)
                                                                                                                                                                                              tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCaattaa                                                                                                                                                                                                                                                  <  8:502544/116‑1 (MQ=255)
                                                                                                                                                                                                 cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGca                                                                                                                                                                                                                                        <  1:82908/123‑1 (MQ=255)
                                                                                                                                                                                                             tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGcc                                                                                                                                                                                                                                                                             <  5:490487/75‑1 (MQ=37)
                                                                                                                                                                                                             tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGcc                                                                                                                                                                                                                                                                             >  6:490487/1‑75 (MQ=37)
                                                                                                                                                                                                             tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGc                                                                                                                                                                                                                                                                              <  5:1104253/74‑1 (MQ=37)
                                                                                                                                                                                                                ccGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGcgcg                                                                   <  1:37675/273‑1 (MQ=255)
                                                                                                                                                                                                                  gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCCCATTGTCGATGCGCTCCGGCAGACTTGGGCGCGAttt                                                               >  2:358443/1‑275 (MQ=255)
                                                                                                                                                                                                                   tCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                                                                                                                                                                                                    <  5:977226/78‑1 (MQ=255)
                                                                                                                                                                                                                   tCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                                                                                                                                                                                                    >  6:977226/1‑78 (MQ=255)
                                                                                                                                                                                                                   tCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTCGCTCATTATAGGGAGCCAGa                                                                                                                                                                                                                                                                          >  7:645509/1‑72 (MQ=255)
                                                                                                                                                                                                                   tCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTCGCTCATTATAGGGAGCCAGa                                                                                                                                                                                                                                                                          <  8:645509/72‑1 (MQ=255)
                                                                                                                                                                                                                        tagtttttGCTCCTGATTTATCCTTTATTACATGGGCAAAAAGGCTTGATCATTTGAAGGGCCCAGACTTTTATTACNAA‑CCGAAATTAATTGGCCAGCCCCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                         <  3:123892/271‑1 (MQ=255)
                                                                                                                                                                                                                         tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                         <  2:349446/274‑1 (MQ=255)
                                                                                                                                                                                                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                                        >  7:1045738/1‑58 (MQ=255)
                                                                                                                                                                                                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                                        <  8:1045738/58‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACt                                                     <  2:196885/266‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGtt                                             >  1:67317/1‑274 (MQ=255)
                                                                                                                                                                                                                                      gATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGAACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTAc                                                             >  1:123892/1‑257 (MQ=255)
                                                                                                                                                                                                                                           ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                                                                                                                                                                                                  >  3:29857/1‑103 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                              <  1:29857/103‑1 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                                                                                                                                       <  8:518309/142‑1 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                                                                                                                                       >  7:518309/1‑142 (MQ=255)
                                                                                                                                                                                                                                                ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg                                                                                                                                         >  4:554140/1‑171 (MQ=255)
                                                                                                                                                                                                                                                    taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt                                                                                                                                     <  2:554141/171‑1 (MQ=255)
                                                                                                                                                                                                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTAc                                                             >  2:471028/1‑239 (MQ=255)
                                                                                                                                                                                                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATATTTGTCGATACGGTTAACCGTGGCtt                      >  4:219157/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                         <  4:471027/239‑1 (MQ=255)
                                                                                                                                                                                                                                                                aggTCTGGTTCATTATAGGGAGCCAGACGTTTATGACCAA‑CCGCAATGAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGTTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTgcggc                 <  2:89255/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaGTCTGGCTCATTATAGGGAGCCAGACTTTTAGTACCAA‑CCTAAATTAATTGGCCAGCACCACATTATGGGAAAACTCTGTTTCACAGTATTTGCATTTGAGCGGGATATCTTTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAGGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTgcggc                 <  1:6132/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 atgcTGGCTCATTATAGGGAGCCAGACTGTTATCAACAA‑ACGTGAGTAAGTGTCCAGCACCTCATGATGTGAAAACTCTGTGTAACAGAATTTGCATTTGAGCGTCCTACCATTGGGGCGTTTTCGGACGGAAAAGCTGGATGAAACCTGTCCGGCATGGCGGATACGGTTGCTGTTCGGCCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCTCCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTgcggcg                <  4:202246/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTAt                                                                                                                                                                                                                               >  7:807352/1‑65 (MQ=255)
                                                                                                                                                                                                                                                                    cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTAt                                                                                                                                                                                                                               <  8:807352/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCAACCACTCCATAGTTGTCGAAACGGTTAACCGTGGCTTGgcgcgccaa            >  4:483701/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                      ggCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                                                                                                                                                                                                  >  7:230650/1‑92 (MQ=255)
                                                                                                                                                                                                                                                                      ggCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATa            >  4:45092/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCaattaat                                                                                                                                                                                                                                                 >  6:342980/1‑43 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCaattaat                                                                                                                                                                                                                                                 <  5:342980/43‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                                                                                                                                       >  6:356216/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                                                                                                                                       <  5:356216/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGATCCAGCCTTTTATAATCAA‑CCACAATTAATTGGCCTGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGAGTTTTCGCACGGCTAAGCTGGATGAAACCGGTTCGGCa                                                                                                                                      <  5:282913/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                                                                                                                                                 <  7:779151/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                                                                                                                                                 >  8:779151/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                                                                                                                                                                                                    <  8:672848/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa                                                                                                                                      >  6:9969/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACATTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                                                                                                                                                                                                    >  7:672848/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                            ttATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat                                                                                                                                                                                       <  7:94214/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg                                                                                                                                    >  7:733519/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTc                                                                                                                                                                      >  7:731830/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTc                                                                                                                                                                      <  8:731830/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACg                                                                                                                                                                 >  4:47233/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                                             >  8:399845/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                                             <  7:399845/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg                                                                                                                                   >  5:1003661/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTTCGGCGCTTACAATGCCa    >  2:508713/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                               taGGGGGCCAAACTTTAATTACCAA‑CCGCAATTAATGTGCCCGCACCACATTATGGGAAAACTCTTTTTCACATTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGt  <  4:137957/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gggAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaa                                                                                                                                                             <  2:47234/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  ggAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTCCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGAt                                                                                                                              >  7:279512/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
AACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAA‑CCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGT  >  NC_000913/4470684‑4471229

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