Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,813,382 T→C L110L (TTA→TTG)  unknown ← D‑glycero‑beta‑D‑manno‑heptose 1‑phosphate adenylyltransferase (EC 2.7.7.70) / D‑glycero‑beta‑D‑manno‑heptose‑7‑phosphate kinase (EC 2.7.1.167)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,813,3820TC100.0% 166.4 / NA 57L110L (TTA→TTG) unknownD‑glycero‑beta‑D‑manno‑heptose 1‑phosphate adenylyltransferase (EC 2.7.7.70) / D‑glycero‑beta‑D‑manno‑heptose‑7‑phosphate kinase (EC 2.7.1.167)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/29);  total (28/29)

TAATAAGGTCGCACCGCGATAGCGTTCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCATGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGATAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACCACCGGCACCGGTGCTTCCGGTGAAATA  >  1/1813136‑1813619
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
taataaGGTCGCACCGCGATAGCGTTCAAAATCCGTCCCTTTCTGATCGATGATCACCGGCACATTGGCGTTACGGGCAATTTTTATCATTTTTTGTACTTCATTTATTGTGCCTTTGCCGTAATCGGACAAGTTCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGTCGGGACaaaa                                                                                                                                                                                                                                           >  1:550778/1‑251 (MQ=255)
             ccGCGATAGCGTTCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGt                                                                                                                                                                                                                              >  1:304307/1‑251 (MQ=255)
               gcgATAGCGTTCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTaa                                                                                                                                                                                                                            >  1:7667/1‑251 (MQ=255)
                cgATAGCGTTCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCAGTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTaa                                                                                                                                                                                                                            <  1:4181/250‑1 (MQ=255)
                       gTTCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCggg                                                                                                                                                                                                                    >  1:505133/1‑251 (MQ=255)
                        ttCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTaa                                                                                                                                                                                                                            >  1:160299/1‑242 (MQ=255)
                        ttCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCgggt                                                                                                                                                                                                                   >  1:12147/1‑251 (MQ=255)
                        ttCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCgggg                                                                                                                                                                                                                   >  1:648861/1‑250 (MQ=255)
                        ttCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCggg                                                                                                                                                                                                                    <  1:672535/250‑1 (MQ=255)
                                          cGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAGATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGa                                                                                                                                                                                                 <  1:582346/251‑1 (MQ=255)
                                          cGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGa                                                                                                                                                                                                 <  1:252109/251‑1 (MQ=255)
                                           ggATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGaa                                                                                                                                                                                                >  1:837157/1‑251 (MQ=255)
                                            gATCGATGAGCACCGGCTCATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGaa                                                                                                                                                                                                <  1:613294/250‑1 (MQ=255)
                                                           gCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTGCTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGAttttt                                                                                                                                                                                >  1:468758/1‑251 (MQ=255)
                                                                      ttACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGtaa                                                                                                                                                                     <  1:636692/251‑1 (MQ=255)
                                                                        aCGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGtaata                                                                                                                                                                   <  1:500709/251‑1 (MQ=255)
                                                                          gggCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTaataat                                                                                                                                                                 >  1:36702/1‑251 (MQ=255)
                                                                           ggCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGAACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATg                                                                                                                                                                >  1:454374/1‑251 (MQ=255)
                                                                               aTTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATg                                                                                                                                                                >  1:261652/1‑247 (MQ=255)
                                                                                     aTCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAAcgc                                                                                                                                                      <  1:258221/251‑1 (MQ=255)
                                                                                     aTCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAAcgc                                                                                                                                                      <  1:202062/251‑1 (MQ=255)
                                                                                         ttttttGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGtaataa                                                                                                                                                                  >  1:442178/1‑235 (MQ=255)
                                                                                                    tCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGaa                                                                                                                                                                                                >  1:881699/1‑194 (MQ=255)
                                                                                                     cATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTaa                                                                                                                                                                                                                            >  1:741857/1‑165 (MQ=255)
                                                                                                     cATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATcc                                                                                                                                      <  1:784493/251‑1 (MQ=255)
                                                                                                     cATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATcc                                                                                                                                      <  1:368941/251‑1 (MQ=255)
                                                                                                     cATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGC‑AATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAAcgc                                                                                                                                                      <  1:71097/234‑1 (MQ=255)
                                                                                                              tgCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGGCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGtt                                                                                                                             <  1:432493/251‑1 (MQ=255)
                                                                                                                   ttGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGcc                                                                                                                        <  1:100487/251‑1 (MQ=255)
                                                                                                                     gCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTC‑AAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGc                                                                                                                                                                                                        >  1:484712/1‑168 (MQ=255)
                                                                                                                     gCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTAACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCa                                                                                                                      >  1:367478/1‑251 (MQ=255)
                                                                                                                      ccGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAg                                                                                                                     >  1:706127/1‑251 (MQ=255)
                                                                                                                         tAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCaa                                                                                                                  <  1:803222/251‑1 (MQ=255)
                                                                                                                                aCAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCAc                                                                                                            <  1:77848/250‑1 (MQ=255)
                                                                                                                                   aGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGtaataa                                                                                                                                                                  <  1:604571/193‑1 (MQ=255)
                                                                                                                                                gTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGag                                                                                            >  1:492756/1‑249 (MQ=255)
                                                                                                                                                              cACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGc                                                                              <  1:421536/250‑1 (MQ=255)
                                                                                                                                                                  ttCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCtt                                                                                                                                           <  1:901496/185‑1 (MQ=255)
                                                                                                                                                                  ttCTAGTTTTTTCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATcc                                                                                                                                      >  1:159254/1‑190 (MQ=255)
                                                                                                                                                                           tttCAGCAAATCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCAc                                                                <  1:669001/251‑1 (MQ=255)
                                                                                                                                                                                   aaTCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACg                                                                                              >  1:287788/1‑213 (MQ=255)
                                                                                                                                                                                    aTCCGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACg                                                                          >  1:659094/1‑232 (MQ=255)
                                                                                                                                                                                      ccGCACTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACg                                                      >  1:72940/1‑250 (MQ=255)
                                                                                                                                                                                           cTTGTCACGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGccc                                                                                                                       <  1:536690/180‑1 (MQ=255)
                                                                                                                                                                                               tcacgTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCAtt                                             <  1:140216/246‑1 (MQ=255)
                                                                                                                                                                                               tcacgTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCAtt                                             <  1:294894/246‑1 (MQ=255)
                                                                                                                                                                                               tcacgTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGACGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAGCGCCGCGCCCGTTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCAtt                                             <  1:334323/246‑1 (MQ=255)
                                                                                                                                                                                                cacgTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCAtttt                                           <  1:505482/248‑1 (MQ=255)
                                                                                                                                                                                                cacgTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCAtttt                                           <  1:249556/248‑1 (MQ=255)
                                                                                                                                                                                                                     aaaaTCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACcaccggca                      <  1:746648/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCg                                                                                                                                              <  1:107114/130‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTCAAGGCGTAATAATTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACcaccggca                      >  1:727946/1‑250 (MQ=255)
                                                                                                                                                                                                                                 taataaTTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACCACCGGCACCGGTGCTTCCg          >  1:367300/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ataatTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGt                                                                                                                                                                                                                        <  1:35310/43‑1 (MQ=255)
                                                                                                                                                                                                                                      atTGTTGATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATg                                                                                                                                                                >  1:629007/1‑96 (MQ=255)
                                                                                                                                                                                                                                        tgttgATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACCACCGGCACCGGTGCTTCCGGTGAAAt   >  1:337892/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gttgATGGCGGGACAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACCACCGGCACCGGTGCTTCCGGTGAAATa  >  1:752293/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
TAATAAGGTCGCACCGCGATAGCGTTCAAAATCCGTCCCTTTCGGATCGATGAGCACCGGCACATTGGCGTTACGGGCAATTTGGATCATTTTTTGTACTTCATTTAATGTGCCTTTGCCGTAATCGGACAAGATCAAAGCGCCGTAATTTTTGACCGCACTTTCTAGTTTTTTCAGCAAATCCGCACTTGTCATGTTCTGGAAATCTTCTTCAAAATCAAGGCGTAATAATTGTTGATGGCGGGATAAAATACGCAATTTGGTAATTGTCGGGTGGGTAGCAAGCTGTACGAAGTCGCAAGCGATTTTTTGTTGTTGTAATAATGAGGTTAACGCCGCGCCCGCTTCATCCTGACCGGTTAAGCCCAGCAATTGCACCGGCACATTTAACGAAGCGATGTTCATTGCCACGTTTGCCGCACCGCCTGCACGCTCCTCATTTTGCTGAACACGCACCACCGGCACCGGTGCTTCCGGTGAAATA  >  1/1813136‑1813619

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