Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 2,058,578 G→A A220A (GCG→GCA)  unknown → Membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*12,058,5780GA100.0% 180.0 / NA 61A220A (GCG→GCA) unknownMembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/25);  total (36/25)

CATGACGGGCTATTATTCCCCGGTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCGGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCAGGCGTGCCTTTGGTGCCGATGGCGTCCGTGGCGTCCGACCGCGGTG  >  1/2058332‑2058822
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
cATGACGGGCTATTATTCCCCGGTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGAc                                                                                                                                                                                                                                                   <  1:18310/251‑1 (MQ=255)
       ggCTATTATTCCCCGGTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACgg                                                                                                                                                                                                                                            >  1:526752/1‑251 (MQ=255)
         cTATTATTCCCCGGTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACgg                                                                                                                                                                                                                                            >  1:677414/1‑249 (MQ=255)
         cTATTATTCCCCGGTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACgg                                                                                                                                                                                                                                            >  1:613632/1‑249 (MQ=255)
                      gtattggcaccggcgagtaCCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATgcgg                                                                                                                                                                                                                              >  1:761215/19‑250 (MQ=255)
                      gTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATgcgg                                                                                                                                                                                                                              <  1:134455/250‑1 (MQ=255)
                         aTTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATgcggc                                                                                                                                                                                                                             <  1:881439/248‑1 (MQ=255)
                          ttCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGcggt                                                                                                                                                                                                                          >  1:516593/1‑250 (MQ=255)
                                      gCACCCAACAAGGCAAATTTAACCAACCAATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATAGACAACTTTCTGCTGGGCGTGCAAGACAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTAATTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCctgc                                                                                                                                                                                                               <  1:517617/249‑1 (MQ=255)
                                           caacaaGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGa                                                                                                                                                                                                        <  1:118922/251‑1 (MQ=255)
                                               aaGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGAc                                                                                                                                                                                                    <  1:663235/251‑1 (MQ=255)
                                                               accGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCgaaag                                                                                                                                                                                    >  1:852585/1‑251 (MQ=255)
                                                               accGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTAAGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCgaaag                                                                                                                                                                                    >  1:356309/1‑251 (MQ=255)
                                                                  gATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGa                                                                                                                                                                                                        >  1:213230/1‑228 (MQ=255)
                                                                   aTTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCctgctg                                                                                                                                                                                                             <  1:795872/222‑1 (MQ=255)
                                                                    tttATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGaaaaa                                                                                                                                                                               >  1:837063/1‑251 (MQ=255)
                                                                       aTGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATg                                                                                                                                                                            >  1:595096/1‑251 (MQ=255)
                                                                       aTGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATg                                                                                                                                                                            >  1:199884/1‑251 (MQ=255)
                                                                                           cGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGa                                                                                                                                                        <  1:17889/251‑1 (MQ=255)
                                                                                                             aTTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAAcc                                                                                                                                       >  1:222858/1‑250 (MQ=255)
                                                                                                                      ggTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACg                                                                                                                                <  1:265314/248‑1 (MQ=255)
                                                                                                                           aTTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGAAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCaa                                                                                                                         >  1:425341/1‑250 (MQ=255)
                                                                                                                                gagGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGAt                                                                                                                                                       >  1:466997/1‑215 (MQ=255)
                                                                                                                                  ggTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGtt                                                                                                                 >  1:704303/1‑251 (MQ=255)
                                                                                                                                   gTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATGTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGtt                                                                                                                 <  1:435562/250‑1 (MQ=255)
                                                                                                                                        ggTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATg                                                                                                                                                                            >  1:280290/1‑186 (MQ=255)
                                                                                                                                        ggTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAgcg                                                                                                           >  1:440033/1‑251 (MQ=255)
                                                                                                                                              gAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAAcc                                                                                                      >  1:524991/1‑250 (MQ=255)
                                                                                                                                              gAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGAC‑AAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAAccc                                                                                                     >  1:216809/1‑250 (MQ=255)
                                                                                                                                                 ttGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCct                                                                                                   >  1:277204/1‑250 (MQ=255)
                                                                                                                                                   ggCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTag                                                                                                               <  1:561051/236‑1 (MQ=255)
                                                                                                                                                   ggCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCctc                                                                                                  >  1:749707/1‑249 (MQ=255)
                                                                                                                                                    gCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTa                                                                                               >  1:24244/1‑251 (MQ=255)
                                                                                                                                                     cGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTa                                                                                               <  1:40094/250‑1 (MQ=255)
                                                                                                                                                     cGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAACGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTa                                                                                               <  1:449434/250‑1 (MQ=255)
                                                                                                                                                       tACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATg                                                                                             <  1:767414/250‑1 (MQ=255)
                                                                                                                                                           gCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑t                                                                                           >  1:657127/1‑247 (MQ=255)
                                                                                                                                                              aTTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCttt                                                                                    >  1:583304/1‑251 (MQ=255)
                                                                                                                                                               ttCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTag                                                                                                               <  1:408454/224‑1 (MQ=255)
                                                                                                                                                                   aTGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGc                                                                                                                      >  1:787509/1‑213 (MQ=255)
                                                                                                                                                                           cAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATgc                                                                                                                                                                                                                                >  1:352510/1‑99 (MQ=255)
                                                                                                                                                                             aCTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATGTTTTTCTTTAAAAATGATCCGTa                                                                      >  1:28935/1‑251 (MQ=255)
                                                                                                                                                                             aCTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTa                                                                      >  1:764735/1‑250 (MQ=255)
                                                                                                                                                                             aCTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGCGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGcc                                                                                                                     >  1:680903/1‑204 (MQ=255)
                                                                                                                                                                              cTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATc                                                                          <  1:397003/245‑1 (MQ=255)
                                                                                                                                                                                  ctgctgGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGaaa                                                                <  1:615889/251‑1 (MQ=255)
                                                                                                                                                                                  ctgctgGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGaaa                                                                <  1:656396/251‑1 (MQ=255)
                                                                                                                                                                                           gTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTc                                                                                                                                                                     >  1:453338/1‑142 (MQ=255)
                                                                                                                                                                                           gTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACg                                                                                                                                <  1:771550/179‑1 (MQ=255)
                                                                                                                                                                                           gTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTagag                                                                                                             <  1:785635/198‑1 (MQ=255)
                                                                                                                                                                                            tGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAgg                                                       <  1:283136/250‑1 (MQ=255)
                                                                                                                                                                                            tGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAgg                                                       <  1:855808/250‑1 (MQ=255)
                                                                                                                                                                                                  gCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCa                                                >  1:335915/1‑251 (MQ=255)
                                                                                                                                                                                                     gTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAgg                                                       >  1:354590/1‑241 (MQ=255)
                                                                                                                                                                                                        gTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTc                                                                                                                                                                     >  1:296975/1‑129 (MQ=255)
                                                                                                                                                                                                        gTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTagag                                                                                                             >  1:28839/1‑185 (MQ=255)
                                                                                                                                                                                                        gTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCAGGCGTg                                          >  1:261850/1‑251 (MQ=255)
                                                                                                                                                                                                         ttATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCAGGCGTGc                                         <  1:831626/251‑1 (MQ=255)
                                                                                                                                                                                                          tATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGa                                                                                                                                                                                                        <  1:703156/92‑1 (MQ=255)
                                                                                                                                                                                                                                                ttACGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCAGGCGTGCCTTTGGTGCCGATGGCGTCCGTGGCGTCCGACCGCGgtg  <  1:626249/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCGCAGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCAGGCGTg                                          >  1:577228/1‑209 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CATGACGGGCTATTATTCCCCGGTGATTCATGCCCGTCGCACCCAACAAGGCAAATTTAACCAACCGATTTATGCACTTCCTTCGAAAAAACGTTTCTCCCGTGCCGAAATTTATGCCGGTGCATTGAGAGGTAAAGGTTTGGAATTGGCGTACAGCGATTCTATGATCGACAACTTTCTGCTGGGCGTGCAAGGCAGTGGTTATGTGGATTTTGGTGAAGGCAATTTGAATTATTTCGCTTACGCGGGACAAAACGGCTACCAATATGCGGCGGTCGGTCGCCTGCTGGTGGAAGACGGTGAAATTCCGAAAGAAAAAATGTCTATTCAGGCAATTCGCGATTGGGCGAATGCCAACCCTTCACGGGTGCAAAGCCTGTTAGAGCGCAACCCCTCTTATG‑TTTTCTTTAAAAATGATCCGTACGGAAAAGTAAAAGGTGCGGCAGGCGTGCCTTTGGTGCCGATGGCGTCCGTGGCGTCCGACCGCGGTG  >  1/2058332‑2058822

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