Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 516,770 A→G V26A (GTA→GCA)  unknown ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1516,7700AG100.0% 224.8 / NA 79V26A (GTA→GCA) unknownhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (43/36);  total (43/36)

TAATAAGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATACGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCTCATCAATTTTACTGTTATCGGCGTGTGGTTCGTTATCTGAAATGAAA  >  1/516524‑517014
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
taataaGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGa                                                                                                                                                                                                                                                  >  1:745245/1‑251 (MQ=255)
  ataaGGCACAGTGTCGAAGGTGAGATTCAATTTTGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATg                                                                                                                                                                                                                                                >  1:570890/1‑251 (MQ=255)
  ataaGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATg                                                                                                                                                                                                                                                >  1:868749/1‑251 (MQ=255)
  ataaGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATg                                                                                                                                                                                                                                                >  1:479641/1‑251 (MQ=255)
  ataaGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATg                                                                                                                                                                                                                                                >  1:28869/1‑251 (MQ=255)
    aaGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGtt                                                                                                                                                                                                                                              >  1:621354/1‑251 (MQ=255)
    aaGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGtt                                                                                                                                                                                                                                              >  1:286967/1‑251 (MQ=255)
            gtgtCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGAt                                                                                                                                                                                                                                           <  1:228219/246‑1 (MQ=255)
            gtgtCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGt                                                                                                                                                                                                                                      >  1:468445/1‑251 (MQ=255)
             tgtCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGt                                                                                                                                                                                                                                      <  1:840626/250‑1 (MQ=255)
             tgtCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTAGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGGCCGATTCTGGCACTGAAATCGTGCAGGCTGTTATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATg                                                                                                                                                                                                                                       <  1:809726/249‑1 (MQ=255)
              gtCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTg                                                                                                                                                                                                                                    >  1:150705/1‑251 (MQ=255)
                  aaGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAgg                                                                                                                                                                                                                                 <  1:838961/250‑1 (MQ=255)
                          aTTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCaaa                                                                                                                                                                                                                         <  1:768169/250‑1 (MQ=255)
                                tttGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGcac                                                                                                                                                                                                                  <  1:72902/251‑1 (MQ=255)
                                          aaCTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGt                                                                                                                                                                                                        <  1:739186/251‑1 (MQ=255)
                                                tCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGACTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACaaa                                                                                                                                                                                                  >  1:321996/1‑251 (MQ=255)
                                                        cAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGt                                                                                                                                                                                           <  1:340080/250‑1 (MQ=255)
                                                              gAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAAc                                                                                                                                                                                     >  1:637573/1‑250 (MQ=255)
                                                                    gCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGt                                                                                                                                                                                                        <  1:878129/225‑1 (MQ=255)
                                                                       gTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACa                                                                                                                                                                           >  1:339132/1‑251 (MQ=255)
                                                                        tCGATAATGCCGCTCAATTCCAAGCGTTTTCTCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACa                                                                                                                                                                           <  1:536238/250‑1 (MQ=255)
                                                                        tCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACa                                                                                                                                                                           <  1:420258/250‑1 (MQ=255)
                                                                                    cTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTcc                                                                                                                                                               <  1:793714/250‑1 (MQ=255)
                                                                                    cTCAATTCCAAGCGTTTTCCCGCGGAGTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTcc                                                                                                                                                               <  1:485142/250‑1 (MQ=255)
                                                                                        aTTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAAt                                                                                                                                                           >  1:570171/1‑250 (MQ=255)
                                                                                         ttCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCg                                                                                                                                                         >  1:778828/1‑251 (MQ=255)
                                                                                                cGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAgg                                                                                                                                                    <  1:763731/249‑1 (MQ=255)
                                                                                                cGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGt                                                                                                                                                   >  1:794853/1‑250 (MQ=255)
                                                                                                      cccGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTaaaaa                                                                                                                                            <  1:675684/251‑1 (MQ=255)
                                                                                                      cccGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAACAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTaaaaa                                                                                                                                            >  1:453029/1‑251 (MQ=255)
                                                                                                           ggATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTAt                                                                                                                                                                    <  1:472700/222‑1 (MQ=255)
                                                                                                                                 aaaTCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTcc                                                                                                                                                               >  1:618820/1‑205 (MQ=255)
                                                                                                                                             cTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAgg                                                                                                                                                    >  1:747716/1‑204 (MQ=255)
                                                                                                                                             cTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAAc                                                                                                                        >  1:607230/1‑232 (MQ=255)
                                                                                                                                             cTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCaa                                                                                                     >  1:302122/1‑251 (MQ=255)
                                                                                                                                             cTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCaa                                                                                                     >  1:570470/1‑251 (MQ=255)
                                                                                                                                               gTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTGACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAAtt                                                                                                   >  1:44909/1‑251 (MQ=255)
                                                                                                                                               gTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAAtt                                                                                                   >  1:497180/1‑251 (MQ=255)
                                                                                                                                               gTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGATGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGTAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATTGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAAtt                                                                                                   >  1:144135/1‑251 (MQ=255)
                                                                                                                                                  aTCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCaa                                                                                                >  1:748949/1‑251 (MQ=255)
                                                                                                                                                      aCGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAg                                                                                            >  1:333308/1‑251 (MQ=255)
                                                                                                                                                           ggCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGt                                                                                                                                                                                                             >  1:848632/1‑133 (MQ=255)
                                                                                                                                                           ggCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTcc                                                                                                                                                               <  1:249655/179‑1 (MQ=255)
                                                                                                                                                                gACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATc                                                                                                                                                          >  1:678702/1‑179 (MQ=255)
                                                                                                                                                                gACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACt                                                                                                                       >  1:108383/1‑214 (MQ=255)
                                                                                                                                                                gACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGa                                                                                             >  1:522854/1‑240 (MQ=255)
                                                                                                                                                                         tGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGttt                                                                                                                                                                                         >  1:856401/1‑139 (MQ=255)
                                                                                                                                                                           ccATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATaa                                                                                                                                                                                                                                         >  1:880043/1‑89 (MQ=255)
                                                                                                                                                                           ccATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTcc                                                                                                                                                               <  1:260467/163‑1 (MQ=255)
                                                                                                                                                                           ccATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATc                                                                                                                                                          <  1:280204/168‑1 (MQ=255)
                                                                                                                                                                           ccATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAatat                                                                       <  1:86212/251‑1 (MQ=255)
                                                                                                                                                                           ccATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTCTTATCAATTCAAGGAGCTTATTTTCATGCAGAAatat                                                                       >  1:30344/1‑251 (MQ=255)
                                                                                                                                                                            cATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAatat                                                                       >  1:727888/1‑250 (MQ=255)
                                                                                                                                                                                    tttttATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAAc                                                                                                                        <  1:568187/193‑1 (MQ=255)
                                                                                                                                                                                     ttttATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAggg                                                             <  1:877071/251‑1 (MQ=255)
                                                                                                                                                                                            attattGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCAc                                                                  >  1:718178/1‑239 (MQ=255)
                                                                                                                                                                                             ttattGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATtctc                                                                                                                                                                <  1:282488/144‑1 (MQ=255)
                                                                                                                                                                                             ttattGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCAc                                                                  <  1:594095/238‑1 (MQ=255)
                                                                                                                                                                                             ttattGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTa                                                     >  1:873011/1‑251 (MQ=255)
                                                                                                                                                                                                ttGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  <  1:898590/251‑1 (MQ=255)
                                                                                                                                                                                                ttGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  <  1:686143/251‑1 (MQ=255)
                                                                                                                                                                                                ttGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  <  1:330494/251‑1 (MQ=255)
                                                                                                                                                                                                ttGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  <  1:770804/251‑1 (MQ=255)
                                                                                                                                                                                                ttGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGACACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  <  1:559901/251‑1 (MQ=255)
                                                                                                                                                                                                  gtgCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  >  1:100058/1‑249 (MQ=255)
                                                                                                                                                                                                  gtgCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCt                                                 >  1:896788/1‑250 (MQ=255)
                                                                                                                                                                                                        aGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCAcc                                                                 <  1:900973/228‑1 (MQ=255)
                                                                                                                                                                                                         gCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTc                                                        <  1:619287/236‑1 (MQ=255)
                                                                                                                                                                                                             gCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCTCATCAATTTTa                                      >  1:155790/1‑250 (MQ=255)
                                                                                                                                                                                                               aaTTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCTCATCAATTTTACTg                                   <  1:205510/251‑1 (MQ=255)
                                                                                                                                                                                                                       gATTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  >  1:276228/1‑228 (MQ=255)
                                                                                                                                                                                                                        aTTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTa                                                     <  1:70392/224‑1 (MQ=255)
                                                                                                                                                                                                                        aTTTCATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  >  1:151573/1‑227 (MQ=255)
                                                                                                                                                                                                                            cATTTTCGCCGATGGCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCAc                                                                                                              <  1:737225/163‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCAAGGGCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCTCATCAATTTTACTGTTATCGGCGTGTGGTTCGTTATc            >  1:122802/1‑248 (MQ=255)
                                                                                                                                                                                                                                              gggCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  >  1:891949/1‑205 (MQ=255)
                                                                                                                                                                                                                                               ggCAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGc                                                  <  1:450775/204‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAAATGCGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCTCATCAATTTTACTGTTATCGGCGTGTGGTTCGTTATCtgaaatgaaa  <  1:791682/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
TAATAAGGCACAGTGTCGAAGGTGAGATTCAATTTGGCGGAAAACTGATCGGTAGGCAGTTTGAATAAGCCGTCGATAATGCCGCTCAATTCCAAGCGTTTTCCCGCGGATTGTCCGATTCTGGCACTGAAATCGTGCAGACTGTAATCCACGCTGGCAAGACCGGGGATGCCATAATTATTTTTATAATTATTGTGCCGAGCCAGCAATTGGTTGATTTCATTTTCGCCGATGGCAAGGGCAAATACGGATGTTGATAATGTTGAGGTCAGCAAAAGTGCACAAAGTAATGTCACAAATTTTGGTTTTAACCAAGCAAACATAGGTATTCTCCTAATCGTAAGGTTTAAAAAAAATAAAAAAACGGTTAAACTAGACCGCACTTTTATCAATTCAAGGAGCTTATTTTCATGCAGAAATATATCACCAGGGTGGTCTTATGCTCATCAATTTTACTGTTATCGGCGTGTGGTTCGTTATCTGAAATGAAA  >  1/516524‑517014

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