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.

AACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGT  >  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)
                                aaGGCCATCGGGTTCGACGATCGCTTTTACGTTAGCCAGCGACGGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaa                                                                                                                                                                                                                                                  <  2:149912/273‑1 (MQ=255)
                                    ccAGCGGCTTCTACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                                             <  2:50908/275‑1 (MQ=255)
                                       ggggtctCGACGATAGTTTTTGCGTTATCCAACGACGTACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAATGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                                                                                                                                                          <  4:469786/269‑1 (MQ=255)
                                       gCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTTCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCNACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATATGTTTTCTCATGCTTTCTNCTTTATTACAGCGGTAAACAGTCTTGCTACTTATATGGAGCCAGANTNTTTTATCAATCCGCaattaatttaccat                                                                                                                                                                                                                                          >  3:37675/1‑269 (MQ=255)
                                       gCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACATCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                                                                                                                                                          >  2:95349/1‑275 (MQ=255)
                                          gCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                       <  2:46375/275‑1 (MQ=255)
                                          gCTTCGACGAGCGATTTTACGTTATCCCGCGACGGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                       <  1:17172/275‑1 (MQ=255)
                                                       cTTTTACCCTATCCAGCGACTGCCCTTGCTGTGCGTCTACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCTTGATCTGACCTGAAGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGTTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                                                                                                                                                                                                          <  1:157627/275‑1 (MQ=255)
                                                          ttACGTTATTCAGCGACTGAATGGCTCGGGCAGCCACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                                                                                                                                                                                                        <  2:294815/274‑1 (MQ=255)
                                                                  tCCGGCGACTGCAGGGGCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGGTTTCGGATTTACCGGGATCTGACCGGAGGGGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                                                                                                                                                                                                <  4:502622/274‑1 (MQ=255)
                                                                  tCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGCTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCGGTCGCGATAGTTTTGCGCAGGCTTTCTCCTTTCTTCCAGCGTTCAACAGTCTAGCTCCTTATCTGGAGCCAGACTTTTATTAACAACCGAACTTAATTGGCCGACACCCCACTAGGGGCAAACATTTTTTc                                                                                                                                                                                                               >  4:89254/1‑275 (MQ=255)
                                                                  tCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAATTCTGGCTCATTATAGGGATCCAGACTTTTATTACCAACCGCAATTAATTGGACAGGCCCACATTATGGGAAAACTCTTTTTc                                                                                                                                                                                                               >  4:72303/1‑275 (MQ=255)
                                                                   ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                              <  1:92508/275‑1 (MQ=255)
                                                                   ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGGTTTCGGATTTACCGGGATCTTACCGAAGGTGATGATCATATTGCCTAGATCAACGCCCTGTACGTAAGGACCGATATCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                              <  2:16371/275‑1 (MQ=255)
                                                                            gACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                                                                                                                                                                                                           >  4:337841/1‑237 (MQ=255)
                                                                                tGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                       <  2:337842/237‑1 (MQ=255)
                                                                                      tgCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                                                                                                                                                                                             <  2:311046/272‑1 (MQ=255)
                                                                                        ctgCGACGTCGGCGGGTACGTCGCCCGGTTTCGGATTTACCGGGTTCTGACCGGAGGGGATGATCATATTGCCCAGAGCAACGCCCTGTACGTAAGGACCGATAGCGGCCGGTGCATTTTCCGTAGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                                                                                                                                                          <  4:430710/272‑1 (MQ=255)
                                                                                                tCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTCCAGCGGTAAAAAGTCTGGCTCATTCACGGGGGTCTGACTTTTATTCCCAACCGGACTTAATTGGCCAGCCCCAATATAAGGGAAAAATTTTTTTCCAAGTATTTGCATTTGCGCCTGATATcct                                                                                                                                                                                  >  4:486674/1‑272 (MQ=255)
                                                                                                                                  ggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                                                                                                                                                                                                               <  5:543853/148‑1 (MQ=255)
                                                                                                                                          cggggggTGATGTTCAGATTTCCCTGTTCACTGCCCTGTCCGTAAGAGCCGATGTGTGCCGGTTCATTTTCCGTCCGGATAGTTTTGCTCATGATTTCTCCTTTATTACAGAGGTAAAAAGTCTGGCTAATTATAGGGAGCCAGGCTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGc                                                                                                                                       <  4:455164/270‑1 (MQ=255)
                                                                                                                                            ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt                                                                                                                                     <  2:217497/275‑1 (MQ=255)
                                                                                                                                             gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                                                                                                                                                              >  4:126981/1‑169 (MQ=255)
                                                                                                                                             gAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa                                                                                                                                      >  2:137958/1‑273 (MQ=255)
                                                                                                                                                 tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                                                                                                                                                          <  2:126982/169‑1 (MQ=255)
                                                                                                                                                   atgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTACGTCGCGCTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCGTTATAGGGGGCCAGACTTTTATTACCAAACGCCATTAATTCGCCAACACCACATTATGGTAAAACTCTTTTTCACAGTATTTGCCTTTGAACGGGATACCCTTGGGCCGTTTTCACACGGACAAGCTGGATAAAACCGGTTCGGCCTGGCCGac                                                                                                                              >  4:447204/1‑274 (MQ=255)
                                                                                                                                                        cATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGACAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGTAAACTCTTTTTCTCAGTATTTGCATTTGAGCGCGATATCATTGTCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAGGGCTGATAccg                                                                                                                          >  1:58904/1‑272 (MQ=255)
                                                                                                                                                          tatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACCGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAACCGCAATTAATTGGCCAGCACCAACTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAACGCGATATCATTGGCCCTTTTTTGCACGGCCACCCTGGATGAAACCGTTTCGGCATGGcctatagcgtt                                                                                                                        >  4:349445/1‑264 (MQ=255)
                                                                                                                                                                   ttgcAACGCCCGGTACGAAAGGACCGCTAGGTGCCGGTGCATTTTCCGTCGCGATAGTTTTGTTCATGATTTTTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACCTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGGGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                                              <  4:31894/272‑1 (MQ=255)
                                                                                                                                                                    aTCAACGACCTTTACGTAAGGACCGTTAGCTTCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTGATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                                              <  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)
                                                                                                                                                                                             aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                                                                                                                                                                                            >  3:82908/1‑123 (MQ=255)
                                                                                                                                                                                              tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaa                                                                                                                                                                                                                                                  <  8:502544/116‑1 (MQ=255)
                                                                                                                                                                                                 cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                                                                                                                                                                                                        <  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)
                                                                                                                                                                                                                ccGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGcgcg                                                                   <  1:37675/273‑1 (MQ=255)
                                                                                                                                                                                                                  gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCCCATTGTCGATGCGCTCCGGCAGACTTGGGCGCGAttt                                                               >  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)
                                                                                                                                                                                                                        tagtttttGCTCCTGATTTATCCTTTATTACATGGGCAAAAAGGCTTGATCATTTGAAGGGCCCAGACTTTTATTACNAACCGAAATTAATTGGCCAGCCCCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                         <  3:123892/271‑1 (MQ=255)
                                                                                                                                                                                                                         tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                         <  2:349446/274‑1 (MQ=255)
                                                                                                                                                                                                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                                        >  7:1045738/1‑58 (MQ=255)
                                                                                                                                                                                                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                                        <  8:1045738/58‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACt                                                     <  2:196885/266‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGtt                                             >  1:67317/1‑274 (MQ=255)
                                                                                                                                                                                                                                      gATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGAACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTAc                                                             >  1:123892/1‑257 (MQ=255)
                                                                                                                                                                                                                                           ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                                                                                                                                                                                                  >  3:29857/1‑103 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                              <  1:29857/103‑1 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                                                                                                                                       <  8:518309/142‑1 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                                                                                                                                       >  7:518309/1‑142 (MQ=255)
                                                                                                                                                                                                                                                ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg                                                                                                                                         >  4:554140/1‑171 (MQ=255)
                                                                                                                                                                                                                                                    taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAt                                                                                                                                     <  2:554141/171‑1 (MQ=255)
                                                                                                                                                                                                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTAc                                                             >  2:471028/1‑239 (MQ=255)
                                                                                                                                                                                                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATATTTGTCGATACGGTTAACCGTGGCtt                      >  4:219157/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                         <  4:471027/239‑1 (MQ=255)
                                                                                                                                                                                                                                                                aggTCTGGTTCATTATAGGGAGCCAGACGTTTATGACCAACCGCAATGAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGTTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTgcggc                 <  2:89255/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaGTCTGGCTCATTATAGGGAGCCAGACTTTTAGTACCAACCTAAATTAATTGGCCAGCACCACATTATGGGAAAACTCTGTTTCACAGTATTTGCATTTGAGCGGGATATCTTTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCAGGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTgcggc                 <  1:6132/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 atgcTGGCTCATTATAGGGAGCCAGACTGTTATCAACAAACGTGAGTAAGTGTCCAGCACCTCATGATGTGAAAACTCTGTGTAACAGAATTTGCATTTGAGCGTCCTACCATTGGGGCGTTTTCGGACGGAAAAGCTGGATGAAACCTGTCCGGCATGGCGGATACGGTTGCTGTTCGGCCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCTCCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTgcggcg                <  4:202246/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                                                                                                                                                                                                               >  7:807352/1‑65 (MQ=255)
                                                                                                                                                                                                                                                                    cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                                                                                                                                                                                                               <  8:807352/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCAACCACTCCATAGTTGTCGAAACGGTTAACCGTGGCTTGgcgcgccaa            >  4:483701/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                      ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                                                                                                                                                                                                  >  7:230650/1‑92 (MQ=255)
                                                                                                                                                                                                                                                                      ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATa            >  4:45092/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                                                 >  6:342980/1‑43 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                                                 <  5:342980/43‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                                                                                                                                       >  6:356216/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                                                                                                                                       <  5:356216/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGATCCAGCCTTTTATAATCAACCACAATTAATTGGCCTGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGAGTTTTCGCACGGCTAAGCTGGATGAAACCGGTTCGGCa                                                                                                                                      <  5:282913/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                                                                                                                                                 <  7:779151/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                                                                                                                                                 >  8:779151/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                                                                                                                                                                                                    <  8:672848/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCa                                                                                                                                      >  6:9969/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          cATTATAGGGAGCCAGACATTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                                                                                                                                                                                                    >  7:672848/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                            ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat                                                                                                                                                                                       <  7:94214/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATg                                                                                                                                    >  7:733519/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTc                                                                                                                                                                      >  7:731830/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTc                                                                                                                                                                      <  8:731830/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACg                                                                                                                                                                 >  4:47233/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                                             >  8:399845/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                                             <  7:399845/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg                                                                                                                                   >  5:1003661/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             tataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTTCGGCGCTTACAATGCCa    >  2:508713/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                               taGGGGGCCAAACTTTAATTACCAACCGCAATTAATGTGCCCGCACCACATTATGGGAAAACTCTTTTTCACATTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGt  <  4:137957/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaa                                                                                                                                                             <  2:47234/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  ggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTCCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGAt                                                                                                                              >  7:279512/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
AACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGT  >  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: