Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,951 A→C 8.6% intergenic (‑38/+35) 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,9510AC 8.6% 289.9 / 20.9 102intergenic (‑38/+35)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (40/51);  new base C (5/5);  total (46/56)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-02

Only 100 of 103 total aligned reads displayed.

TTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATA  >  NC_000913/4470679‑4471219
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                             
ttACAAACACGGTAGTATTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGACGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCAGACGGTGCATTTTCCGGCGCGATAGTTTTTCTCATGATTTGTCCTTTATTACCGCGGTAAAAAGTCTGGCACAtt                                                                                                                                                                                                                                                                             >  1:114855/1‑275 (MQ=255)
     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)
                        tGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGGAAGGACCGATAGCTGCCGGTGCCTTTTCCGTCCCGATAGTTTTGCTCTTTGTTTCCCCTTTATTACCGCGTTAAAAAGTCTGGCTCATTATGGGGTGCCAGTCTTttatttac                                                                                                                                                                                                                                                     >  4:453613/1‑272 (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)
                                                                tACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGTTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAct                                                                                                                                                                                                            <  2:102366/275‑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)
                                                                            cGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGt                                                                                                                                                                                                <  2:340819/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)
                                                                                                                                                                         aTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                                    <  2:72304/274‑1 (MQ=255)
                                                                                                                                                                         aTCAACGACCTTTACGTAAGGACCGTTAGCTTCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTGATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCggg                                                                                                    <  2:486675/274‑1 (MQ=255)
                                                                                                                                                                             aCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                                                                                                                                                                                             <  6:251514/102‑1 (MQ=255)
                                                                                                                                                                             aCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                                                                                                                                                                                             >  5:251514/1‑102 (MQ=255)
                                                                                                                                                                              cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                                                                                                                                                                                            >  8:629064/1‑118 (MQ=255)
                                                                                                                                                                              cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                                                                                                                                                                                            <  7:629064/118‑1 (MQ=255)
                                                                                                                                                                                   tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                              <  5:240981/111‑1 (MQ=38)
                                                                                                                                                                                   tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                              >  6:240981/1‑111 (MQ=38)
                                                                                                                                                                                             gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                              <  8:325521/101‑1 (MQ=38)
                                                                                                                                                                                             gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                                                                                                                                                                                              >  7:325521/1‑101 (MQ=38)
                                                                                                                                                                                                  aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                                                                                                                                                                                  >  3:82908/1‑123 (MQ=255)
                                                                                                                                                                                                   tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaa                                                                                                                                                                                                                                        >  7:502544/1‑116 (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)
                                                                                                                                                                                                                  tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                                                                                                                                                                                             >  5:1047634/1‑65 (MQ=255)
                                                                                                                                                                                                                  tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                                                                                                                                                                                             <  6:1047634/65‑1 (MQ=255)
                                                                                                                                                                                                                     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                                                                                                                                                                                                                                                                <  8:645509/72‑1 (MQ=255)
                                                                                                                                                                                                                        tCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTCGCTCATTATAGGGAGCCAGa                                                                                                                                                                                                                                                                >  7:645509/1‑72 (MQ=255)
                                                                                                                                                                                                                           cggtAGTTTTGTTCATGGTTTCTCCTTTATTACACGGGTAAAAAGGCTTGTTAAGTTTAGGGAGCCAGACGTTTATTACCAACCGCAAGTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTAcc                                                  <  4:95348/271‑1 (MQ=255)
                                                                                                                                                                                                                              tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                               <  2:349446/274‑1 (MQ=255)
                                                                                                                                                                                                                                      cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                               >  4:196884/1‑266 (MQ=255)
                                                                                                                                                                                                                                        cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                              <  8:1045738/58‑1 (MQ=255)
                                                                                                                                                                                                                                        cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                              >  7:1045738/1‑58 (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                                                                                                                                                             >  7:518309/1‑142 (MQ=255)
                                                                                                                                                                                                                                                    tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtttt                                                                                                                                                             <  8:518309/142‑1 (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)
                                                                                                                                                                                                                                                                          tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCAACCACTCCATAGTTGTCGAAACGGTTAACCGTGGCTTGgcgcgccaa  >  4:483701/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                           ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                                                                                                                                                                                        >  7:230650/1‑92 (MQ=255)
                                                                                                                                                                                                                                                                           ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAt                                                                                                                                                                                        <  8:230650/92‑1 (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)
                                                                                                                                                                                                                                                                                attattGGGAGCCAGACATTTATTACCAACCGCAATTAATTGGCCAACACCACATTATGGAAAAACTCTTTTGCACAGTATTTGCATTTGAGCGCGATATCATTTGCACGTTTTCGCACGGCAAAGCTGTATGATACCGGTTCGGCAt                                                                                                                           >  6:282913/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                             
TTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATA  >  NC_000913/4470679‑4471219

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