Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,466,399 G→A G64S (GGT→AGT)  unknown → HlyD‑like membrane fusion protein YhiI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,466,3990GA98.6% 197.1 / ‑5.7 71G64S (GGT→AGT) unknownHlyD‑like membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (34/36);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTTTTACCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGCATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAGGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAACACCGCCAAGCCGCTTAT  >  1/1466152‑1466641
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
gTTTTTACCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTa                                                                                                                                                                                                                                                 >  1:47124/1‑251 (MQ=255)
  ttttACCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAAt                                                                                                                                                                                                                                               <  1:302839/251‑1 (MQ=255)
    ttACCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGa                                                                                                                                                                                                                                             >  1:7945/1‑251 (MQ=255)
      aCCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAg                                                                                                                                                                                                                                           <  1:146489/251‑1 (MQ=255)
      aCCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAg                                                                                                                                                                                                                                           <  1:774035/251‑1 (MQ=255)
      aCCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAg                                                                                                                                                                                                                                           <  1:226584/251‑1 (MQ=255)
          gATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTaaa                                                                                                                                                                                                                                       >  1:75865/1‑251 (MQ=255)
             ctctATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACg                                                                                                                                                                                                                                    >  1:103309/1‑251 (MQ=255)
             ctctATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACg                                                                                                                                                                                                                                    >  1:559666/1‑251 (MQ=255)
               ctATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGc                                                                                                                                                                                                                                   >  1:159844/1‑250 (MQ=255)
               ctATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCa                                                                                                                                                                                                                                  >  1:729389/1‑251 (MQ=255)
               ctATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCa                                                                                                                                                                                                                                  >  1:175796/1‑251 (MQ=255)
                      cATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGCATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCAACCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGCAAAACGCAATCAAc                                                                                                                                                                                                                            >  1:669512/1‑250 (MQ=255)
                          cGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTaa                                                                                                                                                                                                                                                <  1:763034/226‑1 (MQ=255)
                          cGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCa                                                                                                                                                                                                                          >  1:256207/1‑248 (MQ=255)
                                    gcgTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGGAATGAAGTAAAACGCAATCAACCACTGGCACGGCttt                                                                                                                                                                                                             >  1:564856/1‑251 (MQ=255)
                                             tAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGc                                                                                                                                                                                                                     <  1:3087/234‑1 (MQ=255)
                                                         gaTGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTaa                                                                                                                                                                                                                                                >  1:90373/1‑195 (MQ=255)
                                                              aaGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGtt                                                                                                                                                                                   <  1:814626/251‑1 (MQ=255)
                                                                       atTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGca                                                                                                                                                                          >  1:646707/1‑251 (MQ=255)
                                                                        tttCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGcac                                                                                                                                                                         <  1:82845/251‑1 (MQ=255)
                                                                                   gTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAg                                                                                                                                                              >  1:22785/1‑251 (MQ=255)
                                                                                    tttCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAg                                                                                                                                                              <  1:43110/250‑1 (MQ=255)
                                                                                      tCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCGCAGGCACAAAAGaaa                                                                                                                                                           <  1:32316/251‑1 (MQ=255)
                                                                                      tCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGaaa                                                                                                                                                           <  1:237017/251‑1 (MQ=255)
                                                                                         ttattaTCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAgg                                                                                                                                                                      >  1:365812/1‑237 (MQ=255)
                                                                                         ttattaTCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACgg                                                                                                                                                        <  1:620126/251‑1 (MQ=255)
                                                                                         ttattaTCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACgg                                                                                                                                                        >  1:610702/1‑251 (MQ=255)
                                                                                             tatCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAAcc                                                                                                                                                                                                                           <  1:139499/180‑1 (MQ=255)
                                                                                                          ggTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTcc                                                                                                                                                                                                        >  1:744872/1‑186 (MQ=255)
                                                                                                           gTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTccgc                                                                                                                                                                              >  1:612061/1‑211 (MQ=255)
                                                                                                                 ccATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGcc                                                                                                                                <  1:573440/251‑1 (MQ=255)
                                                                                                                                          ggCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCcaa                                                                                                        >  1:365610/1‑250 (MQ=255)
                                                                                                                                          ggCCAGTGAACTGAACGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGAGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTttcttc                                                                                                                                                                                                         >  1:797683/1‑153 (MQ=255)
                                                                                                                                                tgaactgaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCcaacaac                                                                                                    <  1:567793/248‑1 (MQ=255)
                                                                                                                                                 gaactgaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGt                                                                                                >  1:560991/1‑251 (MQ=255)
                                                                                                                                                  aactgaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTATAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTc                                                                                               <  1:652558/251‑1 (MQ=255)
                                                                                                                                                  aactgaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTc                                                                                               <  1:40427/251‑1 (MQ=255)
                                                                                                                                                  aactgaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTc                                                                                               <  1:661309/251‑1 (MQ=255)
                                                                                                                                                  aactgaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTACTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTc                                                                                               <  1:112555/251‑1 (MQ=255)
                                                                                                                                                      gaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTccgc                                                                                                                                                                              >  1:339351/1‑168 (MQ=255)
                                                                                                                                                      gaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGcac                                                                                                                                                                   >  1:767207/1‑179 (MQ=255)
                                                                                                                                                      gaacGGTATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACg                                                                                           >  1:463910/1‑251 (MQ=255)
                                                                                                                                                           gtaTTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTccgc                                                                                                                                                                              <  1:68985/161‑1 (MQ=255)
                                                                                                                                                           gtaTTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACg                                                                                      >  1:580098/3‑251 (MQ=255)
                                                                                                                                                           gtaTTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACg                                                                                      >  1:171996/3‑251 (MQ=255)
                                                                                                                                                           gtaTTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACg                                                                                      >  1:80052/3‑251 (MQ=255)
                                                                                                                                                           gtaTTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAAGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAg                                                                                                 >  1:436390/3‑240 (MQ=255)
                                                                                                                                                                          aCGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAAt                                                                       <  1:483530/251‑1 (MQ=255)
                                                                                                                                                                          aCGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAAt                                                                       <  1:799093/251‑1 (MQ=255)
                                                                                                                                                                                     gAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCAttt                                                                                                                                                                                                 >  1:525401/1‑118 (MQ=255)
                                                                                                                                                                                              cGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAGGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTgcgc                                                     >  1:508832/1‑249 (MQ=255)
                                                                                                                                                                                                        ttGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAAt                                            >  1:334544/1‑248 (MQ=255)
                                                                                                                                                                                                               ccTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAAt                                                                            <  1:204511/209‑1 (MQ=255)
                                                                                                                                                                                                               ccTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCaa                                  >  1:610154/1‑251 (MQ=255)
                                                                                                                                                                                                                cTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGCTGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAAt                                                                       >  1:424099/1‑213 (MQ=255)
                                                                                                                                                                                                                           ccGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAg                      <  1:556335/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGTGTAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGa                     <  1:487025/251‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAg                                                                                                                                                              <  1:772595/113‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTg                                                        <  1:353790/215‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTg                                                        <  1:685262/215‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGaa                    >  1:601843/1‑251 (MQ=255)
                                                                                                                                                                                                                              tgtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAac                   <  1:746004/251‑1 (MQ=255)
                                                                                                                                                                                                                              tgtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAac                   <  1:616452/251‑1 (MQ=255)
                                                                                                                                                                                                                               gtAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAac                   <  1:619325/250‑1 (MQ=255)
                                                                                                                                                                                                                                tAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGTGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAacac                 <  1:427609/251‑1 (MQ=255)
                                                                                                                                                                                                                                tAGAAGAAATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAacac                 <  1:415632/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gaagaaATTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGACAATCTGATTTCAAGCAGCGAATTAGAacac                 <  1:281967/249‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTTTGGTGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAACACCGCCAAg          >  1:448065/1‑250 (MQ=255)
                                                                                                                                                                                                                                               tGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAACACCGCCAAGCCGCTtat  <  1:417264/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tGCAGGAAAGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAACACCGCCAAGCCGCTtat  <  1:843468/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GTTTTTACCCGATCTCTATAATCATTCGCATAAACCGCGTTAAAGTAAGGATACAGAGATGAAAGCAAAATATTTCGTGTTCGGTTTCATTATTATCGGTGCCGTCGGTTTCGCCATTTGGAAAAACAGGCAATTACAGGCCAGTGAACTGAACGGCATTGCGGCAGTTAACGGACGTTTGGAATTAAAACGGCTGGACATTGCCACCCTTTATCCCGGCCGTGTAGAAGAAATTTTGGTGCAGGAAGGTAATGAAGTAAAACGCAATCAACCACTGGCACGGCTTTCTTCCAGCATTTCGCAAGCGCAGGTTTCCGCCGCACAGGCACAAAAGAAACGGGCGGAAGAATCCGTTGCCCGCGCCCTATCAGAAATTGATGCCCGCCAACAACAAGTCAACGTAACGAAATTGGAATTGAATAACGCACAAAAATTGCGCCTCGATAATCTGATTTCAAGCAGCGAATTAGAACACCGCCAAGCCGCTTAT  >  1/1466152‑1466641

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