Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 12.9% intergenic (‑82/+36) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8550CA12.9% 179.0 / 17.5 69intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base A (3/6);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

AGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGC  >  NC_000913/4366563‑4367150
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                       
aGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGACACGTCATCTAAAGTGc                                                                                                                                                                                                                                                                                                 >  2:534926/1‑301 (MQ=255)
   ccAGACGGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                                                                                                                                                                              <  1:621946/301‑1 (MQ=255)
   ccAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                                                                                                                                                                              <  2:55557/301‑1 (MQ=255)
     aGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATGACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtttt                                                                                                                                                                                                                                                                                            <  2:466344/301‑1 (MQ=255)
        ccGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                                                                                                                                                                                         <  2:70935/301‑1 (MQ=255)
           ccAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTGTTAAGTAACTGCTAAAAAGAAAAAAGGAACGTAATCTGACGTGCCTTTTTTAttt                                                                                                                                                                                                                                                                                      >  2:129680/1‑301 (MQ=255)
           ccAGCAACCTGAATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCttg                                                                                                                                                                                                                                                                                      <  1:593687/301‑2 (MQ=255)
             aGCAACCTGGATGGCAGAAATGGCGTCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAGTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGAGCCTTTTTTATTTGt                                                                                                                                                                                                                                                                                    <  1:142758/301‑1 (MQ=255)
              ttctccTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                                                                                   <  2:274391/297‑1 (MQ=255)
              gCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                                                                                   <  2:68368/301‑1 (MQ=255)
                  ccTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGCTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtaagat                                                                                                                                                                                                                                                                               >  1:274088/1‑297 (MQ=255)
                    tgtaTGGCAGAAAGAGCGGCGATAACCGCCATGATAATGGAGATGACTGCGAACGGGATGTTACCGAGTTTAACGACAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTCTTCCCCCCCATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAccc                                                                                                                                                                                                                                                                             <  1:129556/298‑1 (MQ=255)
                               aaTAGCGGCTATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTGACCGGGTTTAACGCCAATATCTGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATGGTTATTTTTAAGTAACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTAc                                                                                                                                                                                                                                                                  <  2:606085/301‑1 (MQ=255)
                               aaTAGCGGCGATAACCGGCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTAc                                                                                                                                                                                                                                                                  <  1:526641/301‑1 (MQ=255)
                                         ttAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                                                                                                                                                                        <  2:281999/300‑1 (MQ=255)
                                           aaCCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCtttat                                                                                                                                                                                                                                                      >  1:281696/1‑299 (MQ=255)
                                            aCCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTAGTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACACGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAAAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGAACTACCCTTTACGATTACTTTTCGCTTtcat                                                                                                                                                                                                                                                     >  1:399218/1‑301 (MQ=255)
                                               gCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatcat                                                                                                                                                                                                                                                  >  1:476304/1‑300 (MQ=255)
                                                 ttttaTAATGGAGTTGACTTCGAACGGGATGTTTCCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGta                                                                                                                                                                                                                                                <  1:47167/297‑1 (MQ=255)
                                                   tGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtaaa                                                                                                                                                                                                                                              >  2:47220/1‑299 (MQ=255)
                                                       aaTGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCATTCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCACGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTCCTGTTCGCTTTCATCAGTATAGCGt                                                                                                                                                                                                                                          >  1:585039/1‑301 (MQ=255)
                                                             gaTTACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCGATTCCCACCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTGGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                                                                                                                                                                                    <  1:77798/301‑1 (MQ=255)
                                                                           ggATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                                                                                                                                                                      >  1:154784/1‑301 (MQ=255)
                                                                            gATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATATCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTTTAAGCCGGGTGc                                                                                                                                                                                                                     >  1:116054/1‑301 (MQ=255)
                                                                                     gggTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                                                                                                                                                                            <  2:175750/301‑1 (MQ=255)
                                                                                              gCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTTCATCCAATTCTGTACGgag                                                                                                                                                                                                   >  1:11239/1‑301 (MQ=255)
                                                                                                  aTAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaaaa                                                                                                                                                                                               <  2:399732/301‑1 (MQ=255)
                                                                                                   aaGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCATTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTTATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                                                                                                                                                                              <  1:534253/300‑1 (MQ=255)
                                                                                                       ggCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCg                                                                                                                                                                                          <  2:313995/301‑1 (MQ=255)
                                                                                                         cAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                                                                                                                                                                        >  2:584313/1‑301 (MQ=255)
                                                                                                                 aGCACCCCCAAGCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                                                                                                                                                                                <  1:212677/301‑1 (MQ=255)
                                                                                                                   cACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTACCACGTTGATATGAGATTGTTATTTTTAAGTTACTGCTCACACGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTCCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTACATCAGAGTCTGTACGGAGCAAATAACGTACAGTTCCGcc                                                                                                                                                                              >  1:570614/1‑300 (MQ=255)
                                                                                                                     ctcccAATTCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                                                                                                                                                                            <  2:11239/299‑1 (MQ=255)
                                                                                                                     cccccAATCCGCCTTCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTGGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTACATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                                                                                                                                                                            <  2:571328/301‑1 (MQ=255)
                                                                                                                        ccAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGACGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTCTTTAACAAGCTGATATTCGATTGTTATTTTTAAGTTAATGCTCACAAGCAAACACGCACGACAACAGACGTGCCTTTTTTATTTGTACTACCCTGTACGCTTACTGTTCGCTTTCATCAGTCTAGCGTTTTGCTTTGTAAGCCGGGCGAAGCCGATTCTGTACGGCCAAAATACCGTCACGTTccccacccat                                                                                                                                                                         >  1:272080/1‑293 (MQ=255)
                                                                                                                          aaTCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCa                                                                                                                                                                       <  2:154935/301‑1 (MQ=255)
                                                                                                                            tCCGCCTGCAAAACCAATATCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                                                                                                                                                                     <  2:116178/301‑1 (MQ=255)
                                                                                                                              ctttccctacatgacgctcttccgatctcAATCTGGCGCCCAAGAAGTTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg                                                                                                                                                                   <  1:466738/273‑1 (MQ=255)
                                                                                                                                  tGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAccc                                                                                                                                                               >  2:503199/1‑300 (MQ=255)
                                                                                                                                    acgtgtgctcttccgatctccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg                                                                                                                                                             <  2:465539/282‑1 (MQ=255)
                                                                                                                                      aaaaCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGtt                                                                                                                                                           <  1:583589/301‑1 (MQ=255)
                                                                                                                                                     ccccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGCCCGCGTTCCAGAACGACTTccc                                                                                                                                            >  1:230945/1‑299 (MQ=255)
                                                                                                                                                       ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcgagatcggaagagcgtcgtg                                                                                                                                          >  1:464975/1‑282 (MQ=255)
                                                                                                                                                          cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGagatcggaagagcacaagtctgaactac                                                                                                                                       >  2:467302/1‑274 (MQ=255)
                                                                                                                                                              ctggcgatcttcttgggcGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAGGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGCTTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACAccct                                                                                                                                   >  1:530598/18‑301 (MQ=255)
                                                                                                                                                              cTGGTGCCCAAGTAGATAGCAAGCAATACTATGATGAGTTCGAAAACTAGCATATTAGCCTTCCTTGTTGTTTAACAAGTTGATAGTAGATTGTTATTTTTAAGGTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTGTTTTATTAGTACTACCCTGTACGATTACTGTGCGCTGTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACActcg                                                                                                                                   <  1:502527/301‑2 (MQ=255)
                                                                                                                                                                  cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTAcc                                                                                                                               >  1:286699/1‑301 (MQ=255)
                                                                                                                                                                           cgtgttctcttccgatcttgatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGc                                                                                                                      <  2:555815/283‑1 (MQ=255)
                                                                                                                                                                           agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGc                                                                                                                      >  2:174215/1‑301 (MQ=255)
                                                                                                                                                                           agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGc                                                                                                                      <  2:476905/301‑1 (MQ=255)
                                                                                                                                                                             aTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGcaa                                                                                                                    >  2:386812/1‑300 (MQ=255)
                                                                                                                                                                                gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCCGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCAAGAACGACATCACGTACACTCTTACCGGATTCGACACAgc                                                                                                                 >  1:225235/1‑300 (MQ=255)
                                                                                                                                                                                             tgatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGCACGACTTCACGTACACTCATACCGGTTTCGGcagaacggaagagcgacgt                                                                                                    >  1:555101/1‑284 (MQ=255)
                                                                                                                                                                                                            ttAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCAccacaacggtg                                                                                     <  2:163889/300‑10 (MQ=255)
                                                                                                                                                                                                                  tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGAt                                                                               <  2:287101/301‑1 (MQ=255)
                                                                                                                                                                                                                    ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGGTATGGTGAACGATGa                                                                             >  2:335922/1‑301 (MQ=255)
                                                                                                                                                                                                                               tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACAGGTTCAgg                                                                  >  2:580824/1‑301 (MQ=255)
                                                                                                                                                                                                                                   ttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTAACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTCACCGATGAACGGGCTCAGgtaaa                                                              >  1:146504/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 ttgctcttccgatctaaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAAc                                <  2:46810/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCg                           <  2:231188/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCATTGTAAGCCGGGTGCATCCGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGcac                        <  2:531271/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAAcagatcggaaacgcgt                 >  1:46757/5‑286 (MQ=255)
                                                                                                                                                                                                                                                                                aaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATAACGATAGAGTTGTAAACGTAACCTTCGAACActtctt                 >  2:14723/5‑301 (MQ=255)
                                                                                                                                                                                                                                                                                   gatctaaGGCACGTCAGATGGCGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGt              <  1:52450/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   gAAATAAGTCACTTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGt              <  2:129894/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTAGTTAGCa         >  1:141289/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        aaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGtagatc         >  2:52518/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                               gTCATCTGACGTTCCTTTTTTATTTGTTCTACCCTTTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGc  <  1:174018/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAAAGTTGTAAACGTAACCTTCGCACACTTCCTTGTTAGCAGTGATGc  >  1:461324/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                       
AGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGC  >  NC_000913/4366563‑4367150

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