Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 730,123 A→C I112M (ATT→ATG)  unknown ← Type I restriction‑modification system, restriction subunit R (EC 3.1.21.3)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1730,1230AC100.0% 221.4 / NA 78I112M (ATT→ATG) unknownType I restriction‑modification system, restriction subunit R (EC 3.1.21.3)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (39/39);  total (39/39)

GTATAGCCTTCCATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAAAATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCATTAGCGGCAAAAGAAAAGATAGTGGCA  >  1/729875‑730369
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gTATAGCCTTCCATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGG‑AAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATa                                                                                                                                                                                                                                                     >  1:187668/1‑251 (MQ=255)
   taGCCTTCCATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATAcc                                                                                                                                                                                                                                                   <  1:461526/251‑1 (MQ=255)
         tCCATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGacaa                                                                                                                                                                                                                                              >  1:344338/1‑250 (MQ=255)
         tCCATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGacaa                                                                                                                                                                                                                                              >  1:165147/1‑250 (MQ=255)
          ccATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGacaaca                                                                                                                                                                                                                                            >  1:470113/1‑251 (MQ=255)
           cATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGAcaacaa                                                                                                                                                                                                                                           >  1:750207/1‑251 (MQ=255)
              atatTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATaa                                                                                                                                                                                                                                        >  1:624436/1‑251 (MQ=255)
                atTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:638747/251‑1 (MQ=255)
                atTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:898420/251‑1 (MQ=255)
                atTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:606589/251‑1 (MQ=255)
                atTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:777530/251‑1 (MQ=255)
                atTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:289271/251‑1 (MQ=255)
                atTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:212507/251‑1 (MQ=255)
                          gCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAAttt                                                                                                                                                                                                                            >  1:670890/1‑251 (MQ=255)
                           cTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAAttt                                                                                                                                                                                                                            >  1:800293/1‑250 (MQ=255)
                                ttCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAAt                                                                                                                                                                                                                                          >  1:260729/1‑231 (MQ=255)
                                ttCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      >  1:43549/1‑235 (MQ=255)
                                ttCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTAt                                                                                                                                                                                                                      >  1:375407/1‑251 (MQ=255)
                                     gTACTACGTTACGCGTATTGGCGTTGCCAAGAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTAtttttt                                                                                                                                                                                                                 >  1:262415/1‑251 (MQ=255)
                                     gTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:727524/230‑1 (MQ=255)
                                     gTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      >  1:526447/1‑230 (MQ=255)
                                     gTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:700298/230‑1 (MQ=255)
                                     gTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTAtttttt                                                                                                                                                                                                                 >  1:702629/1‑251 (MQ=255)
                                     gTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTAtttttt                                                                                                                                                                                                                 >  1:815056/1‑251 (MQ=255)
                                       actacGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:4616/228‑1 (MQ=255)
                                       actacGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTa                                                                                                                                                                                                                >  1:874124/1‑250 (MQ=255)
                                                   gTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCa                                                                                                                                                                                                                                     <  1:336188/217‑1 (MQ=255)
                                                   gTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCaa                                                                                                                                                                                                   >  1:190976/1‑251 (MQ=255)
                                                                      gCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTgg                                                                                                                                                                                >  1:275733/1‑251 (MQ=255)
                                                                                   aaTGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCAt                                                                                                                                                                   <  1:220205/251‑1 (MQ=255)
                                                                                   aaTGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCAt                                                                                                                                                                   <  1:321352/251‑1 (MQ=255)
                                                                                    aTGGTGTTTTGCTTCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:429874/183‑1 (MQ=255)
                                                                                    aTGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCATAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCAt                                                                                                                                                                   >  1:622265/1‑250 (MQ=255)
                                                                                    aTGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:106766/183‑1 (MQ=255)
                                                                                       gtgtTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      >  1:47826/1‑180 (MQ=255)
                                                                                         gtTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAAtc                                                                                                                                                                                                >  1:99770/1‑216 (MQ=255)
                                                                                          ttttGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTAcatc                                                                                                                                                            <  1:612249/251‑1 (MQ=255)
                                                                                               cTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAgt                                                                                                                                                       <  1:876547/251‑1 (MQ=255)
                                                                                                 ccAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAgttg                                                                                                                                                     <  1:676596/251‑1 (MQ=255)
                                                                                                      tCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGc                                                                                                                                                <  1:149064/251‑1 (MQ=255)
                                                                                                             aaGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATAc                                                                                                                                         <  1:597008/251‑1 (MQ=255)
                                                                                                             aaGTAACGATATTACCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATa                                                                                                                                          >  1:189028/1‑250 (MQ=255)
                                                                                                                    gATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTtgttgt                                                                                                                                  >  1:459118/1‑251 (MQ=255)
                                                                                                                     atatTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTa                                                                                                                                 >  1:368637/1‑251 (MQ=255)
                                                                                                                        ttGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGt                                                                                                                               >  1:701732/1‑250 (MQ=255)
                                                                                                                         tGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGtca                                                                                                                             <  1:574909/251‑1 (MQ=255)
                                                                                                                                ggTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:711763/139‑1 (MQ=255)
                                                                                                                                 gTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAg                                                                                                                      <  1:1011/250‑1 (MQ=255)
                                                                                                                                           cATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAg                                                                                                                                                        >  1:641009/1‑206 (MQ=255)
                                                                                                                                           cATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAg                                                                                                           >  1:544043/1‑251 (MQ=255)
                                                                                                                                             ttATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      <  1:9973/126‑1 (MQ=255)
                                                                                                                                             ttATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTAcc                                                                                                                                                                                                              <  1:37145/150‑1 (MQ=255)
                                                                                                                                             ttATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAgttg                                                                                                                                                     <  1:576531/207‑1 (MQ=255)
                                                                                                                                             ttATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACttgttg                                                                                                                                   >  1:708228/1‑225 (MQ=255)
                                                                                                                                             ttATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAg                                                                                                           <  1:589676/249‑1 (MQ=255)
                                                                                                                                             ttATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGaa                                                                                                         <  1:818573/251‑1 (MQ=255)
                                                                                                                                                      cGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGc                                                                                                                                                                                                                                      >  1:718312/1‑117 (MQ=255)
                                                                                                                                                                 ggAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAg                                                                                                           <  1:189004/229‑1 (MQ=255)
                                                                                                                                                                 ggAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTc                                                                                      <  1:681489/250‑1 (MQ=255)
                                                                                                                                                                           tGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCAc                                                                           <  1:178677/251‑1 (MQ=255)
                                                                                                                                                                            gAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCa                                                                            >  1:835937/1‑249 (MQ=255)
                                                                                                                                                                                  aCCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACttg                                                                                                                                      <  1:891013/185‑1 (MQ=255)
                                                                                                                                                                                  aCCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCaaaa                                                                    >  1:300344/1‑251 (MQ=255)
                                                                                                                                                                                  aCCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCaaaa                                                                    <  1:222511/251‑1 (MQ=255)
                                                                                                                                                                                  aCCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCaaaa                                                                    >  1:567076/1‑251 (MQ=255)
                                                                                                                                                                                                       tATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCt                                               <  1:144741/251‑1 (MQ=255)
                                                                                                                                                                                                          cGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAg                                                                                                           >  1:723071/1‑188 (MQ=255)
                                                                                                                                                                                                               aacaaGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTAcatc                                                                                                                                                            >  1:743468/1‑134 (MQ=255)
                                                                                                                                                                                                               aacaaGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGt                                        >  1:168562/1‑250 (MQ=255)
                                                                                                                                                                                                                acaaGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAgttggt                                                                                                                                                   >  1:319275/1‑142 (MQ=255)
                                                                                                                                                                                                                acaaGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGt                                        >  1:546957/1‑249 (MQ=255)
                                                                                                                                                                                                                     gtgtTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTgg                                                                                                 >  1:668793/1‑187 (MQ=255)
                                                                                                                                                                                                                     gtgtTTAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTCTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAg                                                                                                                                                        <  1:783005/132‑1 (MQ=255)
                                                                                                                                                                                                                        tttAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGttcttc                              <  1:144180/251‑1 (MQ=255)
                                                                                                                                                                                                                          tAGCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCAt                            <  1:13028/251‑1 (MQ=255)
                                                                                                                                                                                                                            gCGTTGGCGCATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCATTa                          >  1:226604/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cATCAAACCCGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCATTAGCGGCaaa                  <  1:884260/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGGTTAAAAACATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCATTAGCGGCAAAAGAAAAGAt         <  1:196312/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaaacaTACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGc                                                 >  1:10451/1‑204 (MQ=255)
                                                                                                                                                                                                                                                    aaaacaTACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCATTAGCGGCAAAAGAAAAGATAGTGGCa  >  1:505676/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GTATAGCCTTCCATATATTCTTGATAGCTTTTTTCCAGTACTACGTTACGCGTATTGGCGTTGCCAAAAAGCGTAATGGCATCAATGGTGTTTTGCTCCAAATCGCGGAAAGTAACGATATTGCCAAAGGTTTTGGTCGCATTATAAATTCGATTGGTGCGGGAAAAAGCTTGAATTAACCCGTGATAACGTAGATTTTTATCGACAAACAAGGTGTTTAGCGTTGGCGCATCAAACCCGGTTAAAAAAATACCGACAACAATAAGCAGGTCAATTTCTTTATTTTTTACCCGTTTTGCCAAATCTCGATAATAATTTTGGAAGGATTTTCCATCTACATCATAGTTGGTTGCAAAATACTTGTTGTAGTCATCAATAGCAGATTTTAAGAATTCTTTGGCGCTCGAATTCATGGCTTCCACTTCAAAACTTTCATCGGTAATATCGCCTACCGCGTCTTGTTCTTCATTAGCGGCAAAAGAAAAGATAGTGGCA  >  1/729875‑730369

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